Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,279 C→A 5.2% intergenic (‑75/+233) 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,2790CA 5.2% 797.5 / 6.4 248intergenic (‑75/+233)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (105/124);  new base A (16/3);  total (121/127)
Fisher's exact test for biased strand distribution p-value = 1.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

ATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTT  >  NC_000913/2177131‑2177427
                                                                                                                                                     |                                                                                                                                                    
atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                                                                     <  1:3149648/150‑1 (MQ=255)
 tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcac                                                                                                                                                    <  1:4695714/150‑1 (MQ=255)
 tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATcac                                                                                                                                                    <  1:1793502/150‑1 (MQ=255)
   ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaa                                                                                                                                                  >  2:5767042/1‑150 (MQ=255)
    tGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGTTTTCAAAATATAATGAAATTTTTTGTTTTTTAATTTGGAGTTAACGTTCACaat                                                                                                                                                 >  2:4841165/1‑149 (MQ=255)
     gAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaa                                                                                                                                                >  2:1558276/1‑150 (MQ=255)
       aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAc                                                                                                                                              <  2:1260833/150‑1 (MQ=255)
          ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTAAAAATATAATTAAATTATTTTTTTTTAAATTTAGAGATTAAGATAAAAAAAAAGaa                                                                                                                                           >  2:4074468/1‑149 (MQ=255)
          ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAATAATGTAATTATTTGTTTTTTAATATCGAGATAAAGATCACAAAAACGAc                                                                                                                                           >  2:1072624/1‑150 (MQ=255)
           cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACa                                                                                                                                          >  2:3683484/1‑150 (MQ=255)
            ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAAGAAATTATTTGTTTTTAAATATAGAGATAAAGCTCAAAAAAAAGaaaa                                                                                                                                         >  1:1923488/1‑150 (MQ=255)
              aaCCTCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGATTTGTCGATACAACTTAAATTTTCATATCCTGTCGTTTTTTTTCTTTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAata                                                                                                                                       <  2:587842/150‑1 (MQ=255)
                  gCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATTAAATTATTTTTTTTTTAATATAGAGATAAAGATAACAAAAACGAAAAAATTaa                                                                                                                                   >  1:2178445/1‑150 (MQ=255)
                   cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATATTGAAATTATTTGTTTTTAAATATCGAGATAACGGTCACAAAAACGACAATATGaaa                                                                                                                                  <  2:3923748/150‑1 (MQ=255)
                    aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACTTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGaaaa                                                                                                                                 <  2:5114329/150‑1 (MQ=255)
                                 tgtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGCTTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAAGACAAAATGAAAAATAATCgaggagg                                                                                                                    >  2:2057411/1‑149 (MQ=255)
                                   tgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAAATAAATTAATTTTTTTTAAAAAAACAGATAAACATAAACAAAAACAAAAAATTAAAAATAATAGAGGAGTGa                                                                                                                  >  2:1197414/1‑150 (MQ=255)
                                          ttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATAGAGATAAACATCACAAAAAAGACAAAATGAAAATTATTATATGAGTTAAAGGCaa                                                                                                           >  2:1557661/1‑150 (MQ=255)
                                           tCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAATAATTTTTTTTTTAAAAAAGAGTTAAAGATAACAAAAAAAACAAAAAGAAAATTATTCGTGGTGTGATATGCaaa                                                                                                          >  1:3125180/1‑150 (MQ=255)
                                           tCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGGTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACCATCACAAAAACGACAATAAGAAAATTATTCGAGGAGTGAAAGGCaaa                                                                                                          >  2:3330074/1‑150 (MQ=255)
                                            cAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAAAATGAAAATTTTTAGAGGAGTGAAAGGCaaaa                                                                                                         >  2:310247/1‑150 (MQ=255)
                                              gCATCTGCTTTGTCGATACCACGTACTTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaaaa                                                                                                       <  2:4256202/150‑1 (MQ=255)
                                              gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTTTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaaaa                                                                                                       <  2:1581101/150‑1 (MQ=255)
                                                aTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTGAAAATATAATGAAATTATTTGTTTTTAAATATAGAGATAACGATAACAAAAACGACAATATTAAAATTATTCGAGGGGGGAAAGGCCAAAAAAc                                                                                                     >  2:3645457/1‑150 (MQ=255)
                                                   tGCTTTGTCTATACCACGTAAATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGc                                                                                                  <  2:5164233/150‑1 (MQ=255)
                                                   tGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATTATGAAATTATTTGCTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAATGCAAAAAAACGGc                                                                                                  >  2:4634595/1‑150 (MQ=255)
                                                    gCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTGGATTTAAAAATATAATTAAAATATTTTTTTTTAAAAATAGATATAAAGATAAAAAAAAAAACAAAATGAAAAATATTAGTGGAGTGAAATGCaaaaaaaacgcc                                                                                                 >  1:3212240/1‑145 (MQ=255)
                                                    gCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGTTTTAAAAATATAATTAAATTATTTTTTTTTTAAAATCTAGTAAAAGATAACAAAAAAAACAAAAATTAAATTATTAGGGGAGTGTAAGGaaaaaaaaaagca                                                                                                 >  1:3811585/1‑145 (MQ=255)
                                                            gATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCTAtc                                                                                         <  2:4335125/150‑2 (MQ=255)
                                                                 cACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                                                                    <  2:923424/150‑1 (MQ=255)
                                                                  aCGTACTTTTTCATATCCTGTCGTTTTTTTTCTTTTTAAAAATATAATCAAATTATTTTTTTTTAAATATCGAGATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAg                                                                                   <  1:4189884/150‑1 (MQ=255)
                                                                    gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGcc                                                                                 <  2:3201043/150‑1 (MQ=255)
                                                                     accATTTTCATAGCGTGTCTTTTGTTTTCGATTACAAAATTAAATGAAATTATTGGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCg                                                                                <  1:3290545/148‑1 (MQ=255)
                                                                      tatttttCATATCCTTTCGTTTTTTTTCTATTTCAAAATATATTGAATTTATTTGTTTTTAAATATCTAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGt                                                                               <  2:237638/147‑1 (MQ=255)
                                                                      aCATTTTCATATCCTGTCTTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGt                                                                               <  2:3808041/150‑1 (MQ=255)
                                                                       cATTTTCATATCCTGTCGTTTGTTTTCTATTTAAAAATATAATCAAATTTTTTTTTTTTAAATATCGAGATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                              <  1:607100/150‑1 (MQ=255)
                                                                       cATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtt                                                                              <  2:1752361/150‑2 (MQ=255)
                                                                        aTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                             <  2:2847707/150‑1 (MQ=255)
                                                                           ttCATCTCCAGACGATTGAGTTCCATGACAAAATATCACGAAATTATTGGATTTTAAATATCTAGATAACGACCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGtttttt                                                                          <  2:5742815/150‑6 (MQ=255)
                                                                            tcttatCCTGTGTTTTTTTTTCGATTTCAAAATATAATTAATTTATTTTTTTTTAAATATCGATATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                                                                         <  1:2166427/147‑1 (MQ=255)
                                                                            tCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                                                                         <  2:5069999/150‑1 (MQ=255)
                                                                            tCATAGCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTTTTTGTTTTTAAATATCGAGATACCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                                                                         <  2:2480363/150‑1 (MQ=255)
                                                                                 tCCTGTCGTTTGTTTTCTATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcg                                                                    <  2:5065702/150‑1 (MQ=255)
                                                                                  cgtGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGATTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                                                                   <  2:3519561/148‑1 (MQ=255)
                                                                                  ccTGTCGTTTGTTTTCGTTTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                                                                   <  1:3545113/150‑1 (MQ=255)
                                                                                     gTCTTTTTTTTTCTTTTACAAAATATAATGAAATTTTTTTTTTTAAATTATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg                                                                <  1:1147606/150‑1 (MQ=255)
                                                                                       cGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                                                              <  2:474619/150‑1 (MQ=255)
                                                                                        ttttgtttTCTATTTCAAAATATAATTAAATTATTTTTTTTTAAATATCGATATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                                                             <  1:1553184/149‑1 (MQ=255)
                                                                                        gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                                                             <  2:4091778/150‑1 (MQ=255)
                                                                                           tgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                                                          >  2:2590825/1‑150 (MQ=255)
                                                                                            gtttTCTATTACAACATATACTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                                                         <  1:2244183/150‑1 (MQ=255)
                                                                                             ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                        >  1:4538430/1‑150 (MQ=255)
                                                                                             ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                                                        <  2:5787168/150‑1 (MQ=255)
                                                                                              tttCTATTTCAAAATATAATTAAATTTTTTGTTTTTAAATATCGAGATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                                                       <  1:4893387/150‑1 (MQ=255)
                                                                                               ttCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCTTGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc                                                      <  2:2243484/150‑1 (MQ=255)
                                                                                                tCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCg                                                     <  1:750124/150‑1 (MQ=255)
                                                                                                 cGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg                                                    >  1:2826936/1‑150 (MQ=255)
                                                                                                  gATTTCAAAATATAATGAAATAATTTGATTTTAAATAACGAGATAACGATCACAAAAAAGACAAAATGAAAAGTATTCGAGGATTGAAAGTCAAAAAAACGGCCTCCCGACAGGGAAGCCGAAGCAAAGTGCGCGTGTTTTTATGCCGGa                                                   >  2:2482826/1‑150 (MQ=255)
                                                                                                    tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg                                                 >  2:1659807/1‑150 (MQ=255)
                                                                                                        aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgt                                             >  2:62156/1‑150 (MQ=255)
                                                                                                        aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgt                                             >  2:3497672/1‑150 (MQ=255)
                                                                                                         aaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtg                                            <  2:3887540/150‑1 (MQ=255)
                                                                                                           atatAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGTAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGGGTa                                          >  2:1927608/1‑150 (MQ=255)
                                                                                                             ataATCAAATTATTTGTTTTAAAATATCGATAAAACTAAAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa                                        <  1:503700/150‑1 (MQ=255)
                                                                                                               aaTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa                                      >  1:5752588/1‑150 (MQ=255)
                                                                                                                  gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGGTTTTATGCCGGATGCGGTGTAAACATAc                                   >  2:782848/1‑150 (MQ=255)
                                                                                                                    aaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGt                                 <  2:1138685/150‑1 (MQ=255)
                                                                                                                        attTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                             <  2:3281119/150‑1 (MQ=255)
                                                                                                                         tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAGAGGGAAGCCGTAGCAAAGTGCGCGGGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAg                            >  2:5335570/1‑150 (MQ=255)
                                                                                                                          ttGTTTTTAAATATCGAGATAACGATCACAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTg                          <  1:1012283/150‑1 (MQ=255)
                                                                                                                            gTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGc                         >  2:4716435/1‑150 (MQ=255)
                                                                                                                              ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTa                       >  1:5570256/1‑150 (MQ=255)
                                                                                                                               tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTAt                      <  1:5460935/150‑1 (MQ=255)
                                                                                                                               tttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTGTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTAt                      >  2:3574080/1‑150 (MQ=255)
                                                                                                                                ttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATg                     >  2:5099284/1‑150 (MQ=255)
                                                                                                                                 tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGt                    >  2:2959184/1‑150 (MQ=255)
                                                                                                                                 tAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGt                    >  1:817090/1‑150 (MQ=255)
                                                                                                                                  aaaTATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                   <  2:499989/150‑1 (MQ=255)
                                                                                                                                  aaaTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                   >  1:5372188/1‑150 (MQ=255)
                                                                                                                                  aaaTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTa                   >  1:3738778/1‑150 (MQ=255)
                                                                                                                                   aaTATCGATATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAg                  <  1:2409585/150‑1 (MQ=255)
                                                                                                                                   aaTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAg                  >  1:5801308/1‑150 (MQ=255)
                                                                                                                                    atatCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                 <  2:2723884/150‑1 (MQ=255)
                                                                                                                                    atatCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAgg                 <  1:3018762/150‑1 (MQ=255)
                                                                                                                                      atCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATCGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCa               >  1:67794/1‑150 (MQ=255)
                                                                                                                                      atCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCa               >  2:4182427/1‑150 (MQ=255)
                                                                                                                                       tCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCAt              <  2:1574699/150‑1 (MQ=255)
                                                                                                                                         gagaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGa            <  2:5592332/150‑1 (MQ=255)
                                                                                                                                         gagaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGa            <  2:771266/150‑1 (MQ=255)
                                                                                                                                          agaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGAt           >  2:2036821/1‑150 (MQ=255)
                                                                                                                                          agaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGAt           >  2:4418061/1‑150 (MQ=255)
                                                                                                                                          agaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGAt           >  2:334331/1‑150 (MQ=255)
                                                                                                                                           gaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTCTGTAGGCATGATa          >  2:5436099/1‑150 (MQ=255)
                                                                                                                                           gaTAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATa          <  1:1678884/150‑1 (MQ=255)
                                                                                                                                             tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAg        <  1:4799294/150‑1 (MQ=255)
                                                                                                                                             tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAg        <  1:443895/150‑1 (MQ=255)
                                                                                                                                              aaCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa       >  1:4061303/1‑150 (MQ=255)
                                                                                                                                                 gATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGc    >  1:403944/1‑150 (MQ=255)
                                                                                                                                                  aTCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGCGGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCAGGATAAGACGCt   >  2:282768/1‑150 (MQ=255)
                                                                                                                                                   tCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtt  >  1:4201168/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTT  >  NC_000913/2177131‑2177427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: