Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,245 G→T 23.2% 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,2450GT23.2% 196.0 / 11.3 99intergenic (‑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 (32/39);  new base T (12/14);  total (44/55)
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

ACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATA  >  NC_000913/2177101‑2177393
                                                                                                                                                |                                                                                                                                                    
accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAatt                                                                                                                                                 >  2:1123693/1‑150 (MQ=255)
  caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattat                                                                                                                                               <  2:1327220/150‑1 (MQ=255)
  caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattat                                                                                                                                               <  2:461614/150‑1 (MQ=255)
   accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                              <  1:1008710/150‑1 (MQ=255)
   accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                                              <  1:1008698/150‑1 (MQ=255)
       cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAATAATGAAATtttttttt                                                                                                                                          >  1:1606988/1‑150 (MQ=255)
          gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAAAAAATGAAATTATTTGttttt                                                                                                                                       >  1:1551704/1‑150 (MQ=255)
           cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTa                                                                                                                                      <  2:50407/150‑1 (MQ=255)
             tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTCAAAATATAAATAAATTATTTTTTTTTaaa                                                                                                                                    >  1:346283/1‑150 (MQ=255)
             tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaaa                                                                                                                                    <  2:592514/150‑1 (MQ=255)
                tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                                 <  2:23203/150‑1 (MQ=255)
                tttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGAATTCAAAATATAATTAAATTATTTTTTTTTAAaaat                                                                                                                                 >  1:187270/1‑150 (MQ=255)
                tgtCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                                 <  2:1089636/148‑1 (MQ=255)
                  tCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCg                                                                                                                               <  2:451716/150‑1 (MQ=255)
                     agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                                            <  2:765525/150‑1 (MQ=255)
                        ttGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATaa                                                                                                                         <  1:1474752/150‑1 (MQ=255)
                          gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAATATAATTAAATTATTTTTTTTTTAATATAGAGATaaag                                                                                                                       >  1:619516/1‑148 (MQ=255)
                                 ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaa                                                                                                                <  1:826338/150‑1 (MQ=255)
                                   gAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaa                                                                                                              <  1:971117/150‑1 (MQ=255)
                                    aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                             <  1:1192998/150‑1 (MQ=255)
                                       gCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGa                                                                                                          <  2:443922/150‑1 (MQ=255)
                                       gCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGa                                                                                                          <  2:443930/150‑1 (MQ=255)
                                       gCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCATGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATGGAGATAACGATCACAAAAAcca                                                                                                          >  1:1630182/1‑148 (MQ=255)
                                        ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATATTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGAc                                                                                                         <  2:860350/150‑1 (MQ=255)
                                        ccGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGAc                                                                                                         <  2:860356/150‑1 (MQ=255)
                                          ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATATAAATAAATTAATTTTTTTTAAATATAGAGATAAAGATCACAAAAAAGACaa                                                                                                       >  1:909212/1‑150 (MQ=255)
                                          ggAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAAATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACaa                                                                                                       >  2:1510639/1‑150 (MQ=255)
                                          ggAACCGCATAACCGCCGCGCTGGGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTGTTCAATATCGCGGTAACGATCACAAAAACGACac                                                                                                       >  1:1292152/1‑149 (MQ=255)
                                                gCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTAAAAAAAAAATGAAATTATTTGTTTTTAAATATAGAGATAAAGATAACAAAAACCACAAAATGaa                                                                                                 >  2:1026924/1‑150 (MQ=255)
                                                 cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGaaa                                                                                                <  2:1623522/150‑1 (MQ=255)
                                                 cATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTGAAAATATAATTTAATTATTTTTTTTTAAATTTAGATATAAAGATAACAAAAAAGACAatatttaa                                                                                                >  2:207806/1‑146 (MQ=255)
                                                  aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTATAAATATAAAGAAATAATATTTTTTTAAAATACAATAGAACGAACAAAAAAACGACAACATGtaaa                                                                                               >  1:1444083/1‑150 (MQ=255)
                                                  aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGTTTTAAAAATAAAATGTAATTATTTTTTTTTTATTATGGGGATAACGATTAAACAAAAGACCAAATGgaaa                                                                                               >  1:1228855/1‑150 (MQ=255)
                                                  aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTGAAAAAATAATGTAAATATTTTTTTTTTAATATAGAGATAAAGATCACAAAAAACACAAAATGaaaa                                                                                               >  1:585351/1‑150 (MQ=255)
                                                   tAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAAAAATTAAAATATTTTTTTTTTAAAATAGGGATAACGATAAAAAAAAAGACAAAAAGAAAat                                                                                              >  1:1489471/1‑150 (MQ=255)
                                                        gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAAAGACAAAATGAAAATTATTc                                                                                         >  2:1549720/1‑150 (MQ=255)
                                                        gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTc                                                                                         >  2:1549692/1‑150 (MQ=255)
                                                        gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATTAAATTATTTTTTTTTTAATATAGAGATAAAGATAACAAAAAAGACAATATGAAAattattg                                                                                         >  1:1116864/1‑149 (MQ=255)
                                                         ccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAAGAAATTATTTGTTTTTAAATATCGAGATAACGATAACAAAAACGACAAAATGTAAATTATTCg                                                                                        >  2:565354/1‑150 (MQ=255)
                                                           gcgcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTCGATTTCAAAATAAAATGAAAATATTTTTTTTTTAATATAGAGATAAAGATAACAAAAAACACAATATGAAAAATATTagag                                                                                      >  1:1098099/1‑150 (MQ=255)
                                                             gcTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCgagga                                                                                    <  2:1427878/150‑1 (MQ=255)
                                                                  gCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATTTCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGaa                                                                               <  2:646737/150‑1 (MQ=255)
                                                                  gCGTTGTTCAGCACCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATCAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGaa                                                                               <  2:1362681/150‑1 (MQ=255)
                                                                   cttttttCAGAATCTTCTTTTTCGATACCACTTAATTTTTAATATCGTTTCTTTTTTTTTCTATTTAAAAATATAATTAAATTTTTTTTTTTTAAATATCTATATAACGAACAAAAAAACGACAATATGAAAATTATTCGAGGAGTGaaa                                                                              <  1:1614370/145‑1 (MQ=255)
                                                                   cGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAATTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACCACCAAATGAAAATTATTAGAGGAGTGaaa                                                                              >  2:1670917/1‑150 (MQ=255)
                                                                   cGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAATAATTAAAATAATTTTTTTTAAATATAGAGATAACCATAAAAAAAACGACAAAAAGAAAATTATTCGAGGAGTGaaa                                                                              >  2:1670909/1‑150 (MQ=255)
                                                                          cAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTGGATTTAAAAAAATAAATAAATTTTTTTTTTTGATATATTGTGATAACGATAACCAAAAACACAAAATTTAAATTATTCGAGGGGTGGAAGGcccaa                                                                       >  1:1719841/1‑146 (MQ=255)
                                                                          cAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAATAAAAATAAATTATTTGTTTTTTAATATAGAAATAAAGATAAAAAAAAACACAATAATTAAATTATTAGGGGGGTGGAAAGCaaaa                                                                       >  1:1290443/1‑150 (MQ=255)
                                                                           aGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAAAATGAAAATTATTCGAGGAGTGAAAGGCaaaaa                                                                      >  2:1726721/1‑150 (MQ=255)
                                                                           aGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATAAAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAAAGACAAAATGAAAATTATTCGAGGAGTGAAAGGCaaaaa                                                                      >  2:532886/1‑150 (MQ=255)
                                                                           aGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGATTTAAAAATATAAATAAATTTTTTTTTTTTTAATTTATAGATAAACGTAACAAAAAACACAATATTAAAAATATTAGCTGAGTGTTAGGGaaaaa                                                                      >  1:346621/1‑150 (MQ=255)
                                                                             aaTCTTCTTTGTCGATACCACTTAAATTTCCATATCCTGTCTTTTTTTTTCTTTTTAAAAATATAATGAAATTTTTTTTTTTTAAATATCGATATAACGAAAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaaaaa                                                                    <  1:95719/149‑1 (MQ=255)
                                                                             aaTCGGCTTTTTCTATACCACTTTCTTTTTCATTTCCTGTCTTTTTTTTTCGATTTCAAAATATTTTTAATTTTTTTTTTTTTAAATATCTATATAACTATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaaaaa                                                                    <  1:1464699/149‑1 (MQ=255)
                                                                              aTCTGGTTTGTAGATCCCACGTAAATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATTAATTTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                                   <  1:1163842/150‑1 (MQ=255)
                                                                              aTCTGCTTTGTCTACACCAGGTAATTTTAAATTTCCTGTTTTTTGTTTTCTTTTAAAAAATATAATGAAATTTTTTGTTTTAAAATATCTAGATAACGATCAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                                   <  1:752294/150‑1 (MQ=255)
                                                                              aTCTGCTTTGTCGATACCAATTAATTTTTCATATCCTGTCGTTTTTTTTATTTTTAAAAATATATTTAAATTTTTTTTTTTTAAATATCGATATAACTATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                                   <  1:752439/150‑1 (MQ=255)
                                                                               aatGCTTTGTCGTTACCACGTACTTTTTCATATCCTTTGGTTTTTTTTCTATTTCAAAATTTAATGAAATTTTTTTTTTTTAAATATCTTTATAACTATAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACg                                                                  <  2:592639/148‑1 (MQ=255)
                                                                                        tCGATACCACTTACATTTTCATATCCTTTCTTTTGTTTTCGATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGa                                                         <  2:1083745/150‑1 (MQ=255)
                                                                                        tCGATACCACGTACATTTTCATATCCTGTCTTTTGTTTTCTATTTCAAAATATAATTAAATTATTTTTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGa                                                         <  2:755048/150‑1 (MQ=255)
                                                                                        tCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGa                                                         <  2:1083835/150‑1 (MQ=255)
                                                                                            tACCACTTAAATTTTCATATCCTGTCTTTTGTTTTCGATTAAAAAATATAATTAAATTTTTTTTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgg                                                     <  1:942536/150‑1 (MQ=255)
                                                                                              ccACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGTAAGGCAAAAAAAAGGCCTCCCGATAGGGa                                                   >  2:1484591/1‑150 (MQ=255)
                                                                                              ccACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAAATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAAAATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAAAGGGa                                                   >  2:672114/1‑150 (MQ=255)
                                                                                               cACTTACATTTTATTATCCTGTGTTTTTTTTTCGATTTCAAAATATAATTAATTTATTTTTTTTTAAATATCTAGATAACGATAAAAAAAACTACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                                  <  1:1144804/150‑1 (MQ=255)
                                                                                               cACGTACATTTTCATATCCTGTCTTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                                  <  2:1370567/150‑1 (MQ=255)
                                                                                               cACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                                  <  2:578467/150‑1 (MQ=255)
                                                                                               cACGTAAATTTTCATATCCTGTCTTTTGTTTTCTATTAAAAAATATATTGAAATTTTTTTTTTTTAAATATCTATATAACTATAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                                  <  1:1423029/150‑1 (MQ=255)
                                                                                               cACGTAAATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATATTGAATTTATTTTTTTTTAAATATCTAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGaa                                                  <  1:1144788/150‑1 (MQ=255)
                                                                                                aCGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAAGGCCTCCCGATAGGGAAg                                                 >  2:255889/1‑150 (MQ=255)
                                                                                                 cttACTTTTACTTATCTTGTCTTTTTTTTTCTTTTTAAAAATATATTTAATTTTTTTTTTTTTAAATATCTATATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                                <  1:1404960/148‑1 (MQ=255)
                                                                                                 cGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATCATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGc                                                <  2:1725395/150‑1 (MQ=255)
                                                                                                     attttaCATATCCTTTGTTTTTTTTTCTTTTAAAAAATATATTAAATTTTTTTTTTTTAAAATATCTATATCCCCATCACAACACCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                            <  1:640769/148‑1 (MQ=255)
                                                                                                     aaTTTTACTATCCTGGCGTTTTTTTTCTTTTTCAAAATATAATTAAATTATTTTTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                            <  1:640795/149‑1 (MQ=255)
                                                                                                      aTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                           <  2:1489471/150‑1 (MQ=255)
                                                                                                                ccTGTCTTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                                 <  1:694452/150‑1 (MQ=255)
                                                                                                                 cTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTATCAAAGTGcgcg                                >  2:252188/1‑150 (MQ=255)
                                                                                                                 cTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGAAAATATGAAAATTATTAGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                                >  2:252187/1‑150 (MQ=255)
                                                                                                                  tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCGAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                               >  2:1258264/1‑150 (MQ=255)
                                                                                                                  tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                               >  2:1183232/1‑150 (MQ=255)
                                                                                                                    tCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtgt                             >  2:1684502/1‑150 (MQ=255)
                                                                                                                          tttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                       <  1:639028/149‑1 (MQ=255)
                                                                                                                          tttttCGATTTAAAAATATAATTAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                       <  1:868325/149‑1 (MQ=255)
                                                                                                                            tttCTATTTCAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                     <  1:1559146/150‑1 (MQ=255)
                                                                                                                            tttCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                     <  1:1559143/150‑1 (MQ=255)
                                                                                                                                gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGAGAGGGAAGCCGTAGCAAAGTGCGCGTGTGTTTATGCCGGa                 >  2:1538582/1‑150 (MQ=255)
                                                                                                                                  tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg               >  1:164516/1‑150 (MQ=255)
                                                                                                                                  tttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATg               >  1:1323189/1‑150 (MQ=255)
                                                                                                                                   ttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc              >  1:606169/1‑150 (MQ=255)
                                                                                                                                   ttCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc              >  1:606173/1‑150 (MQ=255)
                                                                                                                                     cAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg            <  1:1505106/150‑1 (MQ=255)
                                                                                                                                      aaaaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgt           >  1:807832/1‑150 (MQ=255)
                                                                                                                                        aaTATAATGAAATTATTTGTTTTTAAATATCGAGATAACGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtgt         <  1:1123693/150‑1 (MQ=255)
                                                                                                                                          tataATTAATTTATTTGTTTTAAAATATAGATATAACGATAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa       <  1:629877/150‑1 (MQ=255)
                                                                                                                                          tataATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa       <  2:1098099/150‑1 (MQ=255)
                                                                                                                                          tataATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa       <  1:1236738/150‑1 (MQ=255)
                                                                                                                                          tataATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaa       <  1:1236790/150‑1 (MQ=255)
                                                                                                                                               tGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa  <  2:187270/150‑1 (MQ=255)
                                                                                                                                               tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCACAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa  >  2:1508766/1‑150 (MQ=255)
                                                                                                                                               tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa  >  2:55407/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATA  >  NC_000913/2177101‑2177393

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