Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,433 C→A 5.8% intergenic (‑256/+141) 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,4330CA 5.8% 214.4 / 3.3 81intergenic (‑256/+141)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base A (3/2);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

AACGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAAC  >  NC_000913/699284‑699582
                                                                                                                                                     |                                                                                                                                                     
aaCGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATgc                                                                                                                                                       <  2:1663243/150‑1 (MQ=255)
  cctAAAGGTGATTGTCGATTGGTGAATAAAGGTCTTGTGTGAGGATATATTTATTGATGAATCGATTATTTTTAGCGGGGTTTTTTTACTGTTTATATTTTACGGGGTCCAAATTGCTGACATAGTGGGATTTGTTCATGCCGGATgcgg                                                                                                                                                     <  2:174794/148‑1 (MQ=255)
      aTGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTCCTATTTGTTCATGCCGGATGCGGCGTg                                                                                                                                                 >  1:1658164/1‑150 (MQ=255)
        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:471809/1‑150 (MQ=255)
                 cGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                      >  2:1329619/1‑150 (MQ=255)
                      ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGTTGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 >  2:1647669/1‑150 (MQ=255)
                      ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 >  2:1585731/1‑150 (MQ=255)
                      ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 >  2:634447/1‑150 (MQ=255)
                      ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 >  2:1244984/1‑150 (MQ=255)
                       gTGAAAAAAGGTCTGGTGTGAGGATATATTTATTTATGAATCGATAATTTTTAGCTGTTTTTATTGAATTTTTATATTTTACGTGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGACGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                                                >  2:383605/1‑150 (MQ=38)
                       gTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                >  1:1754241/1‑150 (MQ=255)
                       gTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                <  1:833059/150‑1 (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                                                 >  1:1266964/1‑150 (MQ=255)
                                                                                                      tACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTg                                                 <  2:1265989/150‑1 (MQ=255)
                                                                                                        cGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtc                                               >  1:672851/1‑150 (MQ=255)
                                                                                                        cGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGtc                                               >  1:1201662/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)
                                                                                                                                                     |                                                                                                                                                     
AACGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAAC  >  NC_000913/699284‑699582

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: