Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,428 C→A 5.2% intergenic (‑224/+84) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,177,4280CA 5.2% 389.8 / 7.4 110intergenic (‑224/+84)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (55/49);  new base A (3/3);  total (58/52)
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.00e+00

Only 100 of 111 total aligned reads displayed.

AAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGT  >  NC_000913/2177282‑2177576
                                                                                                                                                  |                                                                                                                                                    
aaagaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg                                                                                                                                                    >  3:296501/1‑149 (MQ=255)
aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg                                                                                                                                                    >  3:296502/1‑149 (MQ=255)
aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg                                                                                                                                                    >  2:933820/1‑149 (MQ=255)
aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAg                                                                                                                                                    >  3:593875/1‑149 (MQ=255)
       aCAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTGAAGCGTctta                                                                                                                                             <  2:569731/149‑4 (MQ=255)
       aCAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTGAAGCGTctta                                                                                                                                             <  2:400472/149‑4 (MQ=255)
          atatGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCa                                                                                                                                          <  2:1039999/149‑1 (MQ=255)
             tGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGc                                                                                                                                       >  2:928300/1‑149 (MQ=255)
              gAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa                                                                                                                                      >  3:382614/1‑149 (MQ=255)
              gAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa                                                                                                                                      >  3:673392/1‑149 (MQ=255)
               aaaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                                                                                                                                     >  3:440077/1‑149 (MQ=255)
                aaaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATa                                                                                                                                    >  2:210489/1‑149 (MQ=255)
                 aaTTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAg                                                                                                                                   >  3:203886/1‑149 (MQ=255)
                  attattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAgg                                                                                                                                  >  2:871221/1‑149 (MQ=255)
                   ttattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGt                                                                                                                                 >  3:52699/1‑149 (MQ=255)
                    tattCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGtt                                                                                                                                <  1:546160/149‑1 (MQ=255)
                       tCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGcc                                                                                                                             <  4:522458/149‑1 (MQ=255)
                         gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                                                                                                                           >  2:368265/1‑149 (MQ=255)
                         gaggagTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                                                                                                                           >  1:379536/1‑149 (MQ=255)
                              gTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAACTGCGCGTGATTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCa                                                                                                                      <  3:904938/149‑1 (MQ=255)
                                     gCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTaaa                                                                                                               <  3:418243/149‑1 (MQ=255)
                                      cAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGGGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc                                                                                                              <  3:675441/149‑1 (MQ=255)
                                       aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                                                                                                                         <  2:1011189/137‑1 (MQ=255)
                                       aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGAt                                                                                                                         >  1:1011189/1‑137 (MQ=255)
                                        aaaaaaCGGCCTCCCGATAGGGAAGCCGTATCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa                                                                                                                                      >  1:784456/1‑123 (MQ=255)
                                        aaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCa                                                                                                                                      <  2:784456/123‑1 (MQ=255)
                                         aaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGcc                                                                                                           <  2:379536/149‑1 (MQ=255)
                                           aaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGc                                                                                                            <  2:233760/146‑1 (MQ=255)
                                           aaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGc                                                                                                            >  1:233760/1‑146 (MQ=255)
                                            taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTa                                                                                                        >  3:953660/2‑149 (MQ=255)
                                                  cTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTc                                                                                                  >  4:596651/1‑149 (MQ=255)
                                                          aGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACgggg                                                                                          >  3:719002/1‑149 (MQ=255)
                                                           gggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggt                                                                                         <  3:708541/149‑1 (MQ=255)
                                                            ggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGcc                                                                                                                             >  4:756148/1‑112 (MQ=255)
                                                            ggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGcc                                                                                                                             <  3:756148/112‑1 (MQ=255)
                                                            ggAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtg                                                                                        >  1:514323/1‑149 (MQ=255)
                                                               aGCCGTAGCAAAGTTCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGgtggtg                                                                                     >  1:662189/1‑149 (MQ=255)
                                                                  cGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGtgttg                                                                                  >  1:482709/1‑149 (MQ=255)
                                                                    tAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTgttgtt                                                                                <  4:440077/149‑1 (MQ=255)
                                                                    tAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTgttgtt                                                                                <  2:957345/149‑1 (MQ=255)
                                                                     aGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTc                                                                               <  1:871221/149‑1 (MQ=255)
                                                                     aGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTc                                                                               <  2:837391/149‑1 (MQ=255)
                                                                         aaGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGtat                                                                           <  1:253469/149‑1 (MQ=255)
                                                                            tGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAg                                                                        >  3:396217/1‑149 (MQ=255)
                                                                               gcgTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCg                                                                     >  3:481498/1‑149 (MQ=255)
                                                                                cgTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGt                                                                    <  1:725408/149‑1 (MQ=255)
                                                                                  tgtTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTAt                                                                  <  2:196437/149‑1 (MQ=255)
                                                                                   gtTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATc                                                                 >  4:368853/1‑149 (MQ=255)
                                                                                     ttttATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAg                                                               <  4:593875/149‑1 (MQ=255)
                                                                                       ttATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGc                                                             <  4:242931/149‑1 (MQ=255)
                                                                                         aTGCCGGATGCTGTGTAAACATACGTGGCAGTGCTATGTCGGCATTATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTCGTGTTGTTCGTATAAGGCGTATCAGGCGa                                                           <  4:961896/149‑1 (MQ=255)
                                                                                          tGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtg                                                                                        <  2:209868/119‑1 (MQ=255)
                                                                                          tGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtg                                                                                        >  1:209868/1‑119 (MQ=255)
                                                                                          tGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGAt                                                          <  3:547435/149‑1 (MQ=255)
                                                                                           gCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATa                                                         <  4:52699/149‑1 (MQ=255)
                                                                                              ggATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtg                                                                                        <  2:411104/115‑1 (MQ=255)
                                                                                              ggATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtg                                                                                        >  1:411104/1‑115 (MQ=255)
                                                                                               gATGTGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTaa                                                     >  3:197408/1‑149 (MQ=255)
                                                                                                 tGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAAtt                                                   >  1:94038/1‑149 (MQ=255)
                                                                                                      tgtAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTc                                                                               >  3:273879/1‑116 (MQ=255)
                                                                                                      tgtAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTc                                                                               <  4:273879/116‑1 (MQ=255)
                                                                                                      tgtAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAg                                                                        >  4:719790/1‑123 (MQ=255)
                                                                                                      tgtAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAg                                                                        <  3:719790/123‑1 (MQ=255)
                                                                                                      tgtAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCt                                              >  4:698172/1‑149 (MQ=255)
                                                                                                        tAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTAt                                            >  4:761420/1‑149 (MQ=255)
                                                                                                        tAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTAt                                            >  4:469515/1‑149 (MQ=255)
                                                                                                         aaaCATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATc                                           <  1:657209/149‑1 (MQ=255)
                                                                                                         aaaCATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATc                                           <  3:374790/149‑1 (MQ=255)
                                                                                                           aCATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCga                                         <  4:696103/149‑1 (MQ=255)
                                                                                                           aCATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCga                                         <  4:900181/149‑1 (MQ=255)
                                                                                                           aCATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCga                                         <  1:651666/149‑1 (MQ=255)
                                                                                                              tACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTc                                                                               <  3:536500/108‑1 (MQ=255)
                                                                                                              tACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTc                                                                               >  4:536500/1‑108 (MQ=255)
                                                                                                                      aGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTcc                              <  4:197408/149‑1 (MQ=255)
                                                                                                                            aTGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATa                                                         >  2:912048/1‑116 (MQ=255)
                                                                                                                             tGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtggt                                                                                      >  3:320403/1‑86 (MQ=255)
                                                                                                                             tGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGggtggt                                                                                      <  4:320403/86‑1 (MQ=255)
                                                                                                                             tGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACgg                       >  1:702574/1‑149 (MQ=255)
                                                                                                                                aGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGt                    >  1:1019573/1‑149 (MQ=255)
                                                                                                                                   cATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTaa                 <  3:469515/149‑1 (MQ=255)
                                                                                                                                    ctgatgcgACGCTGAAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAAt                >  3:652385/8‑149 (MQ=255)
                                                                                                                                     tGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATc               <  1:28351/149‑1 (MQ=255)
                                                                                                                                      gatgcgACGCTGAAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGt                    >  3:947236/6‑143 (MQ=255)
                                                                                                                                      gatgcgACGCTGAAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGt                    <  4:947236/138‑1 (MQ=255)
                                                                                                                                      gATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCa              >  2:768779/1‑149 (MQ=255)
                                                                                                                                       aTAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAg             >  4:79104/1‑149 (MQ=255)
                                                                                                                                        tAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGt            <  2:451323/149‑1 (MQ=255)
                                                                                                                                         aaGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGtt           <  2:482709/149‑1 (MQ=255)
                                                                                                                                           gACGCTTCAGCGTCGCATCAGGCGTAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTc         >  4:760916/1‑149 (MQ=255)
                                                                                                                                           gACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAGACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTc         >  3:733087/1‑149 (MQ=255)
                                                                                                                                           gACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTc         >  2:1030561/1‑149 (MQ=255)
                                                                                                                                             cGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATACACGTCCTGCACGGCGTTAATCAGTTTCAc       <  4:296502/149‑1 (MQ=255)
                                                                                                                                             cGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCAc       <  4:296501/149‑1 (MQ=255)
                                                                                                                                              gCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACg      >  2:488421/1‑149 (MQ=255)
                                                                                                                                              gCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACg      >  4:42707/1‑149 (MQ=255)
                                                                                                                                               cTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAg                                                        >  1:555287/1‑98 (MQ=255)
                                                                                                                                               cTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAg                                                        <  2:555287/98‑1 (MQ=255)
                                                                                                                                                ttCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGcc    >  2:64889/1‑149 (MQ=255)
                                                                                                                                                 tCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCg   >  3:475034/1‑149 (MQ=255)
                                                                                                                                                  cAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGt  >  2:408519/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGT  >  NC_000913/2177282‑2177576

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: