Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,524,463 Δ1 bp 5.5% intergenic (+95/+18) yncE → / ← ansP ATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,524,4630A. 5.5% 563.2 / 23.3 163intergenic (+95/+18)yncE/ansPATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter
Reads supporting (aligned to +/- strand):  ref base A (83/71);  new base . (3/6);  total (86/77)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 164 total aligned reads displayed.

CTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGA  >  NC_000913/1524316‑1524611
                                                                                                                                                    |                                                                                                                                                    
cTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGACTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAg                                                                                                                                                     >  4:1403421/1‑149 (MQ=255)
 tAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGc                                                                                                                                                    >  1:826746/1‑149 (MQ=255)
   aaCAGCAGTAAGCTATCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCg                                                                                                                                                  <  2:390412/149‑1 (MQ=255)
   aaCAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCg                                                                                                                                                  <  4:1000462/149‑1 (MQ=255)
         agGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCaaa                                                                                                                                            <  2:1103500/149‑2 (MQ=255)
            aaGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCAAattc                                                                                                                                         <  4:889216/149‑2 (MQ=255)
             aGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCAAattca                                                                                                                                        <  2:384378/149‑3 (MQ=255)
             aGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAAtttt                                                                                                                                         <  4:469465/148‑1 (MQ=255)
              gCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCaaattcag                                                                                                                                       <  1:1341974/149‑7 (MQ=255)
              gCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGc                                                                                                                                       <  4:832084/149‑1 (MQ=255)
                tACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCaa                                                                                                                                     <  1:85130/149‑1 (MQ=255)
                  cccAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAAtt                                                                                                                                   <  3:1066994/149‑1 (MQ=255)
                  cccAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAAtt                                                                                                                                   <  4:1178505/149‑1 (MQ=255)
                  cccAGCCAGACGATGTGATTCGTATTCCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAAtt                                                                                                                                   >  3:1244237/1‑149 (MQ=255)
                      gccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTaa                                                                                                                                 <  4:297988/147‑1 (MQ=255)
                      gccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACt                                                                                                                               >  2:1374811/1‑149 (MQ=255)
                       ccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTc                                                                                                                              <  2:150705/149‑1 (MQ=255)
                          gACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGcc                                                                                                                           <  4:16516/149‑1 (MQ=255)
                             gATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGtt                                                                                                                        >  3:8352/1‑149 (MQ=255)
                             gATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGtt                                                                                                                        >  2:1065166/1‑149 (MQ=255)
                              aTGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGttt                                                                                                                       >  1:1064706/1‑149 (MQ=255)
                               tgtgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTc                                                                                                                      >  2:464120/1‑149 (MQ=255)
                               tgtgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTc                                                                                                                      <  3:1454449/149‑1 (MQ=255)
                               tgtgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTc                                                                                                                      <  1:393458/149‑1 (MQ=255)
                                 tgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAg                                                                                                                    <  2:1016344/149‑1 (MQ=255)
                                  gATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAgg                                                                                                                   <  2:561766/149‑1 (MQ=255)
                                        tATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTaaa                                                                                                             <  1:688276/149‑1 (MQ=255)
                                             cgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                                                                        <  3:248449/149‑1 (MQ=255)
                                             cgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                                                                        >  1:1315755/1‑149 (MQ=255)
                                               cTGTAATAGATATGTCGAGGGTGAAAACTGTGATACTTGCCTGCCGTTTAGCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCAGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAtt                                                                                                      >  2:1031214/1‑149 (MQ=255)
                                                tGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAttt                                                                                                     >  4:769129/1‑149 (MQ=255)
                                                 gTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTc                                                                                                    <  3:852144/149‑1 (MQ=255)
                                                  tAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                                                                   >  3:472327/1‑149 (MQ=255)
                                                       gATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGttt                                                                                              >  3:64543/1‑149 (MQ=255)
                                                        atatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGtttt                                                                                             >  4:567812/1‑149 (MQ=255)
                                                         tatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGttttt                                                                                            <  1:488137/149‑1 (MQ=255)
                                                          atGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttc                                                                                           <  2:1308152/149‑1 (MQ=255)
                                                            ggCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                                                                         >  1:1131208/1‑149 (MQ=255)
                                                             gCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttcttc                                                                                        <  1:566190/149‑1 (MQ=255)
                                                               gatggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                          <  1:681067/142‑1 (MQ=255)
                                                                 gggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCt                                                                                    >  1:583765/1‑149 (MQ=255)
                                                                   gTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTc                                                                                  >  4:749122/1‑149 (MQ=255)
                                                                    tGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTcc                                                                                 >  4:1291022/1‑149 (MQ=255)
                                                                     gAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCt                                                                                >  4:1455795/1‑149 (MQ=255)
                                                                       aaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcg                                                                              >  1:1190896/1‑149 (MQ=255)
                                                                       aaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcg                                                                              >  4:1130217/1‑149 (MQ=255)
                                                                        aaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcga                                                                             >  3:328629/1‑149 (MQ=255)
                                                                               agaCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTgg                                                                      <  4:1426783/149‑1 (MQ=255)
                                                                                 aCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcg                                                                    >  1:111280/1‑149 (MQ=255)
                                                                                       ctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAt                                                                                                       <  1:1217009/108‑1 (MQ=255)
                                                                                       ctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAATAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCt                                                                                                  <  2:651987/113‑1 (MQ=255)
                                                                                        tgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGT‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCt                                                             >  2:571217/1‑149 (MQ=255)
                                                                                        tgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGTAGTGCt                                                             >  3:987826/1‑149 (MQ=255)
                                                                                        tgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCt                                                             >  2:243565/1‑149 (MQ=255)
                                                                                          cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCAt                                                                                      <  4:66957/122‑1 (MQ=255)
                                                                                          cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgt                                                           >  1:680430/1‑149 (MQ=255)
                                                                                           ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                                                                                 >  3:764011/1‑110 (MQ=255)
                                                                                           ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg                                                          <  3:948213/149‑1 (MQ=255)
                                                                                            ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTcc                                                                                 >  1:157657/1‑125 (MQ=255)
                                                                                            ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTcc                                                                                 <  2:157657/125‑1 (MQ=255)
                                                                                            ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGa                                                         >  3:1215774/1‑149 (MQ=255)
                                                                                            ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGa                                                         >  4:605070/1‑149 (MQ=255)
                                                                                              ttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAt                                                       <  4:1364853/149‑1 (MQ=255)
                                                                                                 tCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTc                                                    >  3:86209/1‑149 (MQ=255)
                                                                                                  cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                                   <  3:1212472/149‑1 (MQ=255)
                                                                                                  cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                                   <  3:749122/149‑1 (MQ=255)
                                                                                                  cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                                   <  2:962375/149‑1 (MQ=255)
                                                                                                  cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                                   <  1:228988/149‑1 (MQ=255)
                                                                                                   aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAg                                                  >  4:1318248/1‑149 (MQ=255)
                                                                                                      gCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCaa                                               >  2:1369483/1‑149 (MQ=255)
                                                                                                         tATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAAcgc                                            >  3:498996/1‑149 (MQ=255)
                                                                                                         tATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAAcgc                                            >  1:698896/1‑149 (MQ=255)
                                                                                                             cGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTGCATCTGCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGttt                                        <  4:939247/149‑1 (MQ=255)
                                                                                                              ggCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTg                                       >  4:575021/1‑149 (MQ=255)
                                                                                                               gCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTgc                                      <  3:401047/149‑1 (MQ=255)
                                                                                                                  cTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCACTTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTAGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGca                                   <  4:911691/149‑1 (MQ=255)
                                                                                                                   tttCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGcac                                  >  4:350773/1‑149 (MQ=255)
                                                                                                                         cAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCaaa                            <  4:817326/149‑1 (MQ=255)
                                                                                                                           gTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGttt                                        >  1:1292852/1‑135 (MQ=255)
                                                                                                                              aTCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                       <  2:725210/149‑1 (MQ=255)
                                                                                                                              aTCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                       >  1:725210/1‑149 (MQ=255)
                                                                                                                               tCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTGCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc                      >  2:672948/1‑149 (MQ=255)
                                                                                                                               tCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc                      >  2:1227504/1‑149 (MQ=255)
                                                                                                                                 tGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCTGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCta                    >  2:1291490/1‑149 (MQ=255)
                                                                                                                                 tGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCta                    >  1:1406240/1‑149 (MQ=255)
                                                                                                                                    aTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATaa                 <  4:86209/149‑1 (MQ=255)
                                                                                                                                     tttGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                       >  1:394339/1‑142 (MQ=255)
                                                                                                                                      ttGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGcac                                  <  4:185582/130‑1 (MQ=255)
                                                                                                                                       tGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                          <  2:683536/73‑1 (MQ=38)
                                                                                                                                       tGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                          >  1:683536/1‑73 (MQ=38)
                                                                                                                                         cccGGC‑ACGTAGCCGTGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCGTCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGc            <  1:1117022/149‑1 (MQ=255)
                                                                                                                                          ccGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCa           >  3:1370720/1‑149 (MQ=255)
                                                                                                                                          ccGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCa           >  3:1406100/1‑149 (MQ=255)
                                                                                                                                            ggC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCaaa         >  2:1093976/1‑149 (MQ=255)
                                                                                                                                                aCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                          >  3:1048917/1‑65 (MQ=255)
                                                                                                                                                aCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                          <  4:1048917/65‑1 (MQ=255)
                                                                                                                                                aCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATa      >  2:1172246/1‑149 (MQ=255)
                                                                                                                                                aCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATa      <  1:64027/149‑1 (MQ=255)
                                                                                                                                                 cGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATAc     <  2:637188/149‑1 (MQ=255)
                                                                                                                                                    aGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGa  <  1:1076703/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGA  >  NC_000913/1524316‑1524611

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑