Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,440 A→C 5.9% intergenic (‑263/+134) asnB ← / ← umpH asparagine synthetase B/UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,4400AC 5.9% 245.8 / 6.1 86intergenic (‑263/+134)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base C (3/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

GGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCG  >  NC_000913/699292‑699586
                                                                                                                                                    |                                                                                                                                                  
ggTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                                                   >  2:1353859/1‑150 (MQ=255)
  tGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACg                                                                                                                                                 <  2:406771/150‑1 (MQ=255)
  tGATTGTCGAATGGTGAAAAAAGGGCTGGTGGGAGGATATAATTATTGATGAATCGATAATTTTTAGCGGGTATTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCACGCCGCATCCGGCATGAAca                                                                                                                                                 <  2:1142654/150‑2 (MQ=255)
     ttGTCGATTGTGGAAATAAGTTGTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTTAATGTTTATATTTTACGGGTGGCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCTGCGTGAACGCCt                                                                                                                                              <  2:581166/150‑1 (MQ=255)
      ggTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                             <  1:1080389/149‑1 (MQ=255)
       gTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                            >  2:805059/1‑150 (MQ=255)
         cGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                          >  2:1329619/1‑150 (MQ=255)
         cGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                          >  2:471809/1‑150 (MQ=255)
              ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGTTGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     >  2:1647669/1‑150 (MQ=255)
              ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     >  2:634447/1‑150 (MQ=255)
              ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     >  2:1585731/1‑150 (MQ=255)
              ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     >  2:1244984/1‑150 (MQ=255)
               gTGAAAAAAGGTCTGGTGTGAGGATATATTTATTTATGAATCGATAATTTTTAGCTGTTTTTATTGAATTTTTATATTTTACGTGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGACGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                                                    >  2:383605/1‑150 (MQ=38)
               gTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                    <  1:833059/150‑1 (MQ=255)
               gTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                    >  1:1754241/1‑150 (MQ=255)
                tGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGGTTTATTGAATGTTTATATTTTACGGGGGCCAACTTGCTGACAAAGTGCGATGTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcatg                                                                                                                                   <  1:78719/150‑4 (MQ=38)
                ggAAAAAACGGCGGGTGGGAGGATATATTTATTGAAGAATCCCTAATTTTTACCGGGTTTTATTGAAAGTTTATAATTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcatg                                                                                                                                   <  1:1728806/149‑4 (MQ=25)
                 gAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                  >  2:1576459/1‑150 (MQ=255)
                   aaaaaGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaa                                                                                                                                >  2:426319/1‑150 (MQ=255)
                   aaaaaGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGATTATATTTTACGGGGGCCAAATAGCTGACAAAGTGCGATATGTTCATGCCGGATGCGGCGTGAACGCATTATCCGGCCTACaa                                                                                                                                >  1:1476072/1‑150 (MQ=255)
                        ggTCTGGTGTGAGGATATATATATTGATGTATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCa                                                                                                                           <  1:354886/150‑1 (MQ=255)
                          tCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATg                                                                                                                         <  2:346175/150‑1 (MQ=255)
                                  gAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCa                                                                                                                 >  2:764644/1‑150 (MQ=255)
                                        atatTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGTATTCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAatatat                                                                                                           <  1:1574149/150‑1 (MQ=255)
                                           tttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGc                                                                                                        <  1:1472999/150‑1 (MQ=255)
                                            ttattGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCa                                                                                                       <  1:667706/150‑1 (MQ=255)
                                                 gatgaATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGc                                                                                                  >  1:310216/1‑150 (MQ=255)
                                                  atgaATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCAAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCt                                                                                                 <  1:226808/150‑1 (MQ=255)
                                                        cGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAg                                                                                           <  1:1597122/150‑1 (MQ=255)
                                                        cGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAg                                                                                           >  1:76895/1‑150 (MQ=255)
                                                         gATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAgg                                                                                          >  1:1006593/1‑150 (MQ=255)
                                                         gATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAgg                                                                                          >  2:369112/1‑150 (MQ=255)
                                                               ttttAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATACGGCCTACAAAAGCATGAAAATTCAATATATTGAAGGAGCTGCGTAGGCCTGAt                                                                                    >  1:1554822/1‑150 (MQ=255)
                                                                tttAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATa                                                                                   <  1:708610/150‑1 (MQ=255)
                                                                  tAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAg                                                                                 <  2:985950/150‑1 (MQ=255)
                                                                        gTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAgc                                                                           <  1:893687/150‑1 (MQ=255)
                                                                           ttATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCa                                                                        <  1:902421/150‑1 (MQ=255)
                                                                             aTTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATc                                                                      <  2:950085/150‑1 (MQ=255)
                                                                               tGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAg                                                                    <  2:1764996/150‑1 (MQ=255)
                                                                                gAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                   >  1:1108007/1‑150 (MQ=255)
                                                                                gAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                   >  1:1639835/1‑150 (MQ=255)
                                                                                gAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                   >  2:160615/1‑150 (MQ=255)
                                                                                gAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                   >  2:1139158/1‑150 (MQ=255)
                                                                                 aaTGATTATATTTTACGGGGGCCAATTTGCTGACAATGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAATAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                  <  1:70698/150‑1 (MQ=255)
                                                                                  aTGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 <  1:620795/150‑1 (MQ=255)
                                                                                  aTGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:1729309/1‑150 (MQ=255)
                                                                                   tGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAg                                                                >  1:1625546/1‑150 (MQ=255)
                                                                                     tttATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAAATCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                                                              >  2:1109169/1‑150 (MQ=255)
                                                                                      ttATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                                                             <  1:274233/150‑1 (MQ=255)
                                                                                      ttATATTTTACGGGGGCCAAATTGCGGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                                                             <  1:1169959/150‑1 (MQ=255)
                                                                                       tataTTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTg                                                            >  1:687211/1‑150 (MQ=255)
                                                                                          atTTTACGGGCGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCACATCAGGCAGTTTGGCg                                                         <  2:793831/150‑1 (MQ=255)
                                                                                              tACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTg                                                     <  2:1265989/150‑1 (MQ=255)
                                                                                              tACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTg                                                     >  1:1266964/1‑150 (MQ=255)
                                                                                                cGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtc                                                   >  1:1201662/1‑150 (MQ=255)
                                                                                                cGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtc                                                   >  1:672851/1‑150 (MQ=255)
                                                                                                  ggggCCAAATTCCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtcat                                                 >  2:976580/1‑150 (MQ=255)
                                                                                                   gggCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtcatc                                                <  1:634447/150‑1 (MQ=255)
                                                                                                     gCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCag                                              >  2:1582034/1‑150 (MQ=255)
                                                                                                             gCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAAcc                                      <  1:107448/150‑1 (MQ=255)
                                                                                                               tGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCCTTTGACATCAGCGCCAACCAc                                    >  1:681327/1‑150 (MQ=255)
                                                                                                                gACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACg                                   <  1:764644/150‑1 (MQ=255)
                                                                                                                    aaGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCg                               <  1:694897/150‑1 (MQ=255)
                                                                                                                       tGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAg                            <  2:1537488/150‑1 (MQ=255)
                                                                                                                          gATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGACAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                         >  2:1127238/1‑150 (MQ=255)
                                                                                                                              tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGt                     >  1:1679774/1‑150 (MQ=255)
                                                                                                                              tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGt                     >  1:1091497/1‑150 (MQ=255)
                                                                                                                              tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACCTCCGCAGACGTGGt                     >  1:1672306/1‑150 (MQ=255)
                                                                                                                               gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGtt                    <  1:1602090/150‑1 (MQ=255)
                                                                                                                                ttCATGCCGGATGCGTCGTGAACGCCTTATCCTGCCTACAAAAGCATGAAAATTAAATATATTGCAGGAGCTGCGTAGTACTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGATCCAACCACGTCCGCAGACGttgttg                   <  2:1695725/150‑1 (MQ=21)
                                                                                                                                ttCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTg                   <  1:1266439/150‑1 (MQ=255)
                                                                                                                                 tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCGGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGc                  <  2:334352/150‑1 (MQ=255)
                                                                                                                                 tCAAGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATATATTGCAGCAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTgg                  <  1:903321/150‑2 (MQ=21)
                                                                                                                                   aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTCGCGTTTGTCATCAGAGCCAACCACGTCAGCAGACGTGGTTGCTa                >  1:208162/3‑150 (MQ=21)
                                                                                                                                     gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCAGCAGAGCCAACCACGTCCGCAGACGTGGTTGCTAtt              >  2:1604694/1‑150 (MQ=25)
                                                                                                                                      ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc             <  1:130462/150‑1 (MQ=255)
                                                                                                                                      ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc             >  1:784892/1‑150 (MQ=25)
                                                                                                                                       cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa            <  2:1006593/150‑1 (MQ=255)
                                                                                                                                       cGGAAGCGGCGTGAACGCCTTATCCGGACTACAAAAGCATGAAAATTCAATATATTGCATGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa            <  1:1037791/150‑1 (MQ=255)
                                                                                                                                        ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGTTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAg           <  1:978991/150‑1 (MQ=255)
                                                                                                                                             cggcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAAc      >  1:1685383/1‑150 (MQ=255)
                                                                                                                                              ggcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGTCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACg     <  2:1073279/150‑1 (MQ=255)
                                                                                                                                              ggcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACg     >  2:1883799/1‑150 (MQ=255)
                                                                                                                                              ggcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACg     >  2:1800211/1‑150 (MQ=255)
                                                                                                                                              ggcgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACg     >  2:1346722/1‑150 (MQ=255)
                                                                                                                                                cgTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTc   <  1:1353859/150‑1 (MQ=255)
                                                                                                                                                 gTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCTTCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAtcgtcg  <  2:385985/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCG  >  NC_000913/699292‑699586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: