Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,432 G→A 6.4% intergenic (‑255/+142) 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,4320GA 6.4% 234.6 / 4.8 81intergenic (‑255/+142)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base G (38/36);  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 = 6.38e-01

AAACGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAG  >  NC_000913/699283‑699577
                                                                                                                                                     |                                                                                                                                                 
aaaCGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACTGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATg                                                                                                                                                   >  2:1266439/1‑150 (MQ=255)
aaaCGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATg                                                                                                                                                   >  2:1312391/1‑150 (MQ=255)
aaaCGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATg                                                                                                                                                   <  1:1696601/150‑1 (MQ=255)
 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:1244984/1‑150 (MQ=255)
                       ggTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                            >  2:634447/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:76895/1‑150 (MQ=255)
                                                                 cGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAg                                                                                  <  1:1597122/150‑1 (MQ=255)
                                                                  gATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAgg                                                                                 >  2:369112/1‑150 (MQ=255)
                                                                  gATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAgg                                                                                 >  1:1006593/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                                                        >  2:1729309/1‑150 (MQ=255)
                                                                                           aTGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                        <  1:620795/150‑1 (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)
                                                                                                                                                     |                                                                                                                                                 
AAACGTAATGGTGATTGTCGATTGGTGAAAAAAGGTCTGGTGTGAGGATATATTTATTGATGAATCGATAATTTTTAGCGGGTTTTATTGAATGTTTATATTTTACGGGGGCCAAATTGCTGACAAAGTGCGATTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAG  >  NC_000913/699283‑699577

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: