Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,245 G→T 11.4% intergenic (‑41/+267) 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,2450GT11.4% 419.0 / 6.1 154intergenic (‑41/+267)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base G (63/65);  new base T (13/11);  total (76/78)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

TACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATA  >  NC_000913/2177100‑2177393
                                                                                                                                                 |                                                                                                                                                     
tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                                   >  1:2072059/1‑150 (MQ=255)
  ccaccaCTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatta                                                                                                                                                 >  2:1322381/1‑150 (MQ=255)
   caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATGAAattat                                                                                                                                               >  1:1648128/1‑150 (MQ=255)
   caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattat                                                                                                                                               <  2:1164513/150‑1 (MQ=255)
   caccacTTGCATCGTTTAGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattat                                                                                                                                               <  2:54250/150‑1 (MQ=255)
    accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                              <  1:2474009/150‑1 (MQ=255)
    accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                              <  1:521864/150‑1 (MQ=255)
    accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                              <  1:2380152/150‑1 (MQ=255)
            cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTa                                                                                                                                      <  2:14401/150‑1 (MQ=255)
              tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAAATAATGAAATTATTTGTTTTTaaa                                                                                                                                    >  1:1302303/1‑150 (MQ=255)
              tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAAAGAAATTATTTTTTTTTtaa                                                                                                                                    >  1:1302401/1‑150 (MQ=255)
                 tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                                 <  1:2000847/150‑1 (MQ=255)
                 tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                                 <  1:2514092/150‑1 (MQ=255)
                 tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                                 <  1:1236347/150‑1 (MQ=255)
                 tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                                 >  1:307823/1‑150 (MQ=255)
                 tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAATAAAGAAAATATTTGTTTTTAAatat                                                                                                                                 >  1:1230748/1‑150 (MQ=255)
                  ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTTGATTTCAAAAAATAAAGAAATTAATTGTTTTTaaaaaaa                                                                                                                                >  2:339432/1‑146 (MQ=255)
                    cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCga                                                                                                                              <  2:2532137/150‑1 (MQ=255)
                      agaTTGTGAATATAGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAAAGTAATTATTTGTTTTTAAATATCgaga                                                                                                                            >  1:2223570/1‑150 (MQ=255)
                       gaTTGTGAATATTGAATGTCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                           <  2:189863/150‑1 (MQ=255)
                       gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                           <  2:2491963/150‑1 (MQ=255)
                       gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                           <  2:766999/150‑1 (MQ=255)
                       gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                           <  2:2057854/150‑1 (MQ=255)
                       gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                                           <  2:2057830/150‑1 (MQ=255)
                          tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                                        <  1:1687859/150‑1 (MQ=255)
                          tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATTAAATTATTTGTTTTTAAATATGGAGATaaa                                                                                                                        >  1:2229058/1‑149 (MQ=255)
                          tgtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAATTATTTTTTTTTAAATATAGAGATAAc                                                                                                                        >  1:1247910/1‑150 (MQ=255)
                           gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACg                                                                                                                       >  1:1801275/1‑150 (MQ=255)
                               atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                                   >  1:430614/1‑150 (MQ=255)
                                     aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                             <  2:1298225/150‑1 (MQ=255)
                                     aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                             <  2:80663/150‑1 (MQ=255)
                                     aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                             <  2:125123/150‑1 (MQ=255)
                                     aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                             <  1:790132/150‑1 (MQ=255)
                                     aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                             <  1:150775/150‑1 (MQ=255)
                                         ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATCGAGGTAACGATCACAAAAACGAc                                                                                                         >  1:287461/1‑150 (MQ=255)
                                          cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAAGATAACAAAAACGACa                                                                                                        >  1:2111012/1‑150 (MQ=255)
                                          cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTTAATATAGAGATTACGATAACAAAAACGACa                                                                                                        >  1:974646/1‑150 (MQ=255)
                                           ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTAAAAATATAATGAAAATATTTTTTTTTAAATATAGAGATAACGGTGACAAAAAAGACaa                                                                                                       >  1:2444047/1‑150 (MQ=255)
                                            gAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAATGAAATTATTTGTTTTTTAATATGGAGATAAAGATAAAAAAATCGACAat                                                                                                      >  1:1642158/1‑150 (MQ=255)
                                                  cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAAATAATTAAATTATTTGTTTTTAAATATAGAGATAACGATAACAAAAAAGACAATATGaaa                                                                                                >  1:2342929/1‑150 (MQ=255)
                                                  cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAAATATTTTTTTTTAAATATAGAGATAACGATCAAACAAAAGACAAAATGaaa                                                                                                >  1:1523273/1‑150 (MQ=255)
                                                   aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAATTAAATTATTTTTTTTTAAATATCGAGATAAAGAAAACAAAAAAGACAAAATTaaaa                                                                                               >  1:2521612/1‑150 (MQ=255)
                                                           cgcgCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCga                                                                                       <  2:1802950/150‑1 (MQ=255)
                                                           cgcgCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCga                                                                                       <  2:234702/150‑1 (MQ=255)
                                                              gcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCgagga                                                                                    <  2:176547/150‑1 (MQ=255)
                                                              gcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCgagga                                                                                    <  2:2274243/150‑1 (MQ=255)
                                                                tgtgCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGt                                                                                  <  2:1621647/150‑1 (MQ=255)
                                                                   gCGTGGTTCAGCATCTTCTTTTTCGATACCCCTTAATTTTAAATATCCTGTGGTTTTTTTTCTTTTACAAAATATATTTAAATTATTTTTTTTTAAATATCTATATACCTAACACAAACACGACAATATGAAAATTATTCGAGGAGTGaa                                                                               <  1:1044177/150‑1 (MQ=255)
                                                                       tttacatcATCGGCTTTGTAGATACAACTTAAATTTACATACCCTTCCTTTTTTTTTCTATTAAAAAATATAATTAAATTTTTTTTTTTAAAATATCGACATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGc                                                                           <  1:715051/143‑1 (MQ=255)
                                                                         ttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGCTTTAAAAATATAATTAAATTATTTTTTTTTAAATATAGAGAAAACGATCACAAAAAAGAAAATATGAAAATTATTAGAGGTGTGGAAGGGaa                                                                         >  2:2400895/1‑150 (MQ=255)
                                                                         ttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaa                                                                         <  2:2488495/150‑1 (MQ=255)
                                                                         ttCAGCATCTGCTTTGTCGATACCACGTAAATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaa                                                                         <  2:2488595/150‑1 (MQ=255)
                                                                              cATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATTAAAATTATTCGAGGAGTGAAAGGCaaaaaaa                                                                    >  2:1942036/1‑150 (MQ=255)
                                                                                  tGCTTTGTCGATACCACGTAATTTTTCATATCCTGTCTTTTGTTTTCTATTTCAAAATATAATTAAATTATTTTTTTTTAATTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGc                                                                <  2:559974/150‑1 (MQ=255)
                                                                                   gCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGTAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGcc                                                               >  2:569869/1‑150 (MQ=255)
                                                                                   gCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAAAATAATTAAAATATTTTTTTTTTAATATGGAGATAACGATCACAAAAAAGACAAAAATAAAAATATTAGTGGAGTGTAAGGCGAAAAAAAGGcc                                                               >  1:364522/1‑150 (MQ=255)
                                                                                      ttGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAAATAATGAAATTATTTTTTTTTAAATATAGAGATAAAGATAACAAAAAAGAAAATATTAAAAATATTAGATGAGTGTAAGGCAAAAAAAAGGCCTcc                                                            >  2:1982997/1‑150 (MQ=255)
                                                                                            aTACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAg                                                      <  2:1292661/150‑1 (MQ=255)
                                                                                            aTACCACGTAAATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATTAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAg                                                      <  2:883765/150‑1 (MQ=255)
                                                                                                   gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                               <  2:1648128/150‑1 (MQ=255)
                                                                                                     aCATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGt                                             <  2:287461/150‑1 (MQ=255)
                                                                                                          ttCATAGCCTTTCGTTTTTTTTCTATTTAAAATTATAATTAAATTTTTTTTTTTTAAATATCTATATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaa                                        <  1:2072067/150‑1 (MQ=255)
                                                                                                                 ccTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                                 <  1:1021449/150‑1 (MQ=255)
                                                                                                                   tGTCGTTTGTTTTCGATTTCAAAAAATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                               >  2:973637/1‑150 (MQ=255)
                                                                                                                      gttttgtttTCGATTTCAAAATATATTGAATTTATTTGTTTTTAAATTTCGAGATACCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                            <  1:1118639/148‑1 (MQ=255)
                                                                                                                      cGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                            <  1:1904838/150‑1 (MQ=255)
                                                                                                                        tttgtttTCGATTTCAAAATATAATGAAATTATTTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                           >  1:15081/1‑150 (MQ=255)
                                                                                                                        tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                          >  2:84845/1‑150 (MQ=255)
                                                                                                                        tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                          >  1:362474/1‑150 (MQ=255)
                                                                                                                        tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                          >  1:1327764/1‑150 (MQ=255)
                                                                                                                         ttgtttTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt                         <  2:1366090/150‑1 (MQ=255)
                                                                                                                         ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt                         >  2:1369773/1‑150 (MQ=255)
                                                                                                                         ttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt                         >  1:1971681/1‑150 (MQ=255)
                                                                                                                          ttttttCTATTTCAAAATATAATTAATTTATTTTTATTTAAATATCGATATAACTATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTa                        <  1:2383899/148‑1 (MQ=255)
                                                                                                                           gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                       <  2:416347/150‑1 (MQ=255)
                                                                                                                           gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                       <  2:1480765/150‑1 (MQ=255)
                                                                                                                            ttttCGTTTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                      <  2:1112705/150‑1 (MQ=255)
                                                                                                                            ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATg                      <  2:1801273/150‑1 (MQ=255)
                                                                                                                            tgtattttttCAAAATTTTTTTAATTTTTTTGTTTTTATATATCTTTTTAACTTTCAAAATAACGACATTATTAATTTTATTCGAGTTGTGAAAGGAAAAAAAACGGCCTCCCGATAGGTAAGCAGTAGCAAAGTGCGCGTGTTTTTATg                      <  2:2393070/143‑1 (MQ=255)
                                                                                                                             tttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                     >  2:69890/1‑150 (MQ=255)
                                                                                                                             tttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTAATCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                     >  1:1339914/1‑150 (MQ=255)
                                                                                                                             tttCGAATTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGTAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGAGTGTTTTTATGc                     >  1:459359/1‑150 (MQ=255)
                                                                                                                              ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc                    >  2:529530/1‑150 (MQ=255)
                                                                                                                              ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGcc                    >  2:529568/1‑150 (MQ=255)
                                                                                                                                 gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCACAAAAACGGCCTCCCGATAGGGCAGCCGTAGCAAAGGGCGCGTGTGTTTATGCCGGa                 >  2:603829/1‑150 (MQ=255)
                                                                                                                                 gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGa                 >  2:2332577/1‑150 (MQ=255)
                                                                                                                                     tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg             >  1:1632219/1‑150 (MQ=255)
                                                                                                                                     tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg             >  2:764419/1‑150 (MQ=255)
                                                                                                                                       aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAGGTGCGCGTGTTTTTATGCCGGATGCGgt           >  2:1317285/1‑150 (MQ=255)
                                                                                                                                       aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgt           >  1:2186861/1‑150 (MQ=255)
                                                                                                                                          atatAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTa        <  1:2043032/150‑1 (MQ=255)
                                                                                                                                          atatAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTa        <  1:2043041/150‑1 (MQ=255)
                                                                                                                                           tataATGAAATTATTTGTTTTTAAATATCGAGATCACGATCACACAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAGAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa       <  1:318100/150‑1 (MQ=255)
                                                                                                                                           tataATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa       <  1:923040/150‑1 (MQ=255)
                                                                                                                                           tataATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa       <  1:559948/150‑1 (MQ=255)
                                                                                                                                            ataATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa      <  2:7794/150‑1 (MQ=255)
                                                                                                                                            ataATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCATAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa      <  2:1613046/150‑1 (MQ=255)
                                                                                                                                             taATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc     <  2:2450331/150‑1 (MQ=255)
                                                                                                                                                tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACAtg  >  2:2182306/1‑149 (MQ=255)
                                                                                                                                                tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa  >  2:2182311/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATA  >  NC_000913/2177100‑2177393

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: