Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,177,245 G→T 14.5% 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,2450GT14.5% 1014.5 / 31.9 455intergenic (‑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 (169/199);  new base T (39/42);  total (211/244)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 462 total aligned reads displayed.

TTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAC  >  NC_000913/2177096‑2177394
                                                                                                                                                      |                                                                                                                                                     
tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGGTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAATATATAT‑Tg                                                                                                                                                       >  2:1471928/1‑150 (MQ=255)
 ttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGa                                                                                                                                                      <  2:100872/150‑1 (MQ=255)
  tCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGCTTTCAAAATATAA‑TGaa                                                                                                                                                     >  1:1442796/1‑150 (MQ=255)
     accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAATAAAA‑TGaaaat                                                                                                                                                  >  2:1739035/1‑148 (MQ=255)
      ccaccaCTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAatta                                                                                                                                                 <  1:3539371/150‑1 (MQ=255)
       caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAA‑TGAAattat                                                                                                                                                >  1:5267501/1‑150 (MQ=255)
        accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAattatt                                                                                                                                               <  2:7199253/150‑1 (MQ=255)
        accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAattatt                                                                                                                                               <  2:6111572/150‑1 (MQ=255)
        accacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAattatt                                                                                                                                               <  2:6499187/150‑1 (MQ=255)
         ccacGTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTAttt                                                                                                                                              <  2:2346470/150‑1 (MQ=255)
          cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTAAAAATATAA‑AGAAATTAtttt                                                                                                                                             >  1:6303717/1‑149 (MQ=255)
          cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAAATAtttt                                                                                                                                             >  2:4172333/1‑149 (MQ=255)
            cTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑AGAAAATATTTGtt                                                                                                                                           >  1:2771957/1‑150 (MQ=255)
              tGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGtttt                                                                                                                                         <  1:7296249/150‑1 (MQ=255)
               gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATtttttttt                                                                                                                                        >  2:1339059/1‑150 (MQ=255)
                cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTa                                                                                                                                       <  2:5433470/150‑1 (MQ=255)
                cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTa                                                                                                                                       <  1:1806367/150‑1 (MQ=255)
                  tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTaaa                                                                                                                                     <  1:4665131/150‑1 (MQ=255)
                  tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATAAAA‑TTAAATTATTTGTTTTTaaa                                                                                                                                     >  2:8594088/1‑150 (MQ=255)
                  tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAA‑TGAAATTATTTTTTTTTaaa                                                                                                                                     >  2:8500405/1‑150 (MQ=255)
                      ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATc                                                                                                                                 <  2:302474/150‑1 (MQ=255)
                      ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAA‑TGAAATTATTTTTTTTTAAATATc                                                                                                                                 >  2:4722453/1‑150 (MQ=255)
                       tCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCg                                                                                                                                <  2:7181289/150‑1 (MQ=255)
                        cGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TTAAATTATTTTTTTTTAAATATAga                                                                                                                               >  2:5909286/1‑150 (MQ=255)
                         gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TTAAATTATTTGTTTTTAAATATagag                                                                                                                              >  1:6799664/1‑150 (MQ=255)
                                 gAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TTAAATTATTTTTTTTTTAATATAGAGATAAAGAt                                                                                                                      >  1:2369788/1‑150 (MQ=255)
                                 gAATATTGAATGCCGGAACCGCATAACCGCCGCGATGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGAt                                                                                                                      <  2:2640036/150‑1 (MQ=255)
                                  aaTATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATc                                                                                                                     <  1:3364959/150‑1 (MQ=255)
                                   atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                                    <  1:8346966/150‑1 (MQ=255)
                                    tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAgataacgataac                                                                                                                   >  2:2926509/1‑150 (MQ=255)
                                    tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAAATATTTGTTTTTAAATATCGAGATAACGATcac                                                                                                                   >  2:1192819/1‑150 (MQ=255)
                                        gAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaa                                                                                                               <  1:4163323/150‑1 (MQ=255)
                                         aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACtaaaa                                                                                                              <  1:3843057/150‑1 (MQ=255)
                                         aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                              <  2:3369938/150‑1 (MQ=255)
                                         aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                                              <  2:5362285/150‑1 (MQ=255)
                                         aaTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAA‑TGAAATTATTTTTTTTTAAATATAGAGATAACGATCACcaaaa                                                                                                              >  1:5341475/1‑150 (MQ=255)
                                          aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAc                                                                                                             <  2:7095413/150‑1 (MQ=255)
                                            gCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGa                                                                                                           <  2:3897442/150‑1 (MQ=255)
                                              cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTGAAAATATAA‑TTAAAATATTTTTTTTTTAATATAGTGATAAAGATCAaaaaaaacaca                                                                                                         >  1:4817566/1‑145 (MQ=255)
                                              cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCCAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACa                                                                                                         >  2:8131621/1‑150 (MQ=255)
                                                       aTAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAAAATAA‑AGAAAATATTTTTTTTTAAATATCGGGCTAAACATAACCAAAACGCCCAAAAGaaaa                                                                                                >  1:7659654/1‑150 (MQ=255)
                                                          aCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAA‑TGAAAATATTTTTTTTTAAAAATAGAGATAACGATCACAAAAACGACCATATGAAAatta                                                                                             >  2:3164322/1‑150 (MQ=255)
                                                             gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTTTTTTTAAATATCGAGATAACGATAACAAAAAAGACAATATGAAAattatta                                                                                          >  2:5409631/1‑149 (MQ=255)
                                                               cgcgCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCga                                                                                        <  2:6750826/150‑1 (MQ=255)
                                                               cgcgCTGTGCGTTGTTCAGAATCTGTTTTGTCTACACCACTTAATTTTTAATATCCTGTTGTTTGTTTTCTATTAAAAAATATAA‑TTAAATTTTTTTTTTTAAATTTTCGATAAAACTAACAAAAAAACGACAATATGAAAATTATTCga                                                                                        <  1:6731250/150‑1 (MQ=255)
                                                                             ttCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAA‑TTAAATTATTTTTTTTTAAATATAGAGATAACGATAACAAAAACGACAATAAGAAAATTATTCGAGGAGTGAAAGGCaa                                                                          >  2:1770180/1‑150 (MQ=255)
                                                                               cAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGATTTCAAAATATAA‑TGAAAATATTTTTTTTTAAATATAGAGATAACGGTAACAAAAACCGCAAAATTAAAATTATTAGATGAGTGGAAGGCaaaa                                                                        >  2:4613183/1‑150 (MQ=255)
                                                                               cAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAA‑TGAAATTATTTTTTTTTAAATATAGAGATAACGATAACAAAAACGACAATATTAAAATTATTCGAGGAGTGAAAGGCaaaa                                                                        >  2:4674722/1‑150 (MQ=255)
                                                                                 gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTGAAAATATAA‑TGAAAATAATTTTTTTTAAAAATAGAGATAAAGATAACAAAAAACACAAAATGAAAATTAATGAAGGAGTGAAAGGCaaaaaa                                                                      >  2:6581436/1‑150 (MQ=255)
                                                                                 gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTGGTTTTGAAAATATAA‑TTAAAATATTTTTTTTTAAATATAGATATAAAGATAACAAAAAAGAAAAAATTAAAATTATTAGATGAGTGTAATTCaaaaaa                                                                      >  2:8060813/1‑150 (MQ=255)
                                                                                   aTCTGCTTTGTCGATACCACGTAATTTTTAATTTCCTGTCGTTTTTTTTCTATTTAAAAATTTAT‑TTAATTTTTTTTTTTTTAAATATGTAGATAACTATCAAAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAc                                                                    <  1:4003844/150‑1 (MQ=255)
                                                                                    tCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTGAAAATATAA‑TTAAAATAATTTTTTTTTAATATAGAGATAAAGAACAAAAAAAACACCATATGAAAAATATTAGATGGGGGGAAGGGaaaaaaaag                                                                   >  1:5002909/1‑148 (MQ=255)
                                                                                     ctgctgTGTCTATACCACGTACATTTTCATACCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACgg                                                                  <  2:8361373/150‑1 (MQ=255)
                                                                                               gATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                                        <  2:5608592/150‑1 (MQ=255)
                                                                                               gATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATa                                                        <  2:1663260/150‑1 (MQ=255)
                                                                                               gATACCACGTACATGTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTGGTTAATATCGAGATAACGATCACAAAAACGACAATCTGAAAAGTATTCGCGGTGGGAAAGGCGAAATAACGGCCTCCCGATa                                                        >  2:4957230/1‑150 (MQ=255)
                                                                                                  acatcttACTTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATTTAT‑TTAACTTTTTTTTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAggg                                                     <  2:2454899/144‑1 (MQ=255)
                                                                                                   ccACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGCGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAgggc                                                    >  2:2174447/1‑149 (MQ=255)
                                                                                                    cccGTACATTTTCATATCCTGTGGTTTGTTTTCGTTTTCAAAATATAA‑TGAATTTATTTGTTTTTAATTATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAACGGCAAAAAAACGGCCTCCCGATAGGGaa                                                   <  1:5132431/148‑1 (MQ=255)
                                                                                                    cACGTACATTTTCATATCCTGTCGTTTTTTTTGGATTTGAAAAAAAAA‑TTAAAATATTTTTTTTTAAATATAGAGATAAAGATCACAAAAAACACAATATTAAAATTATTAGATGATGGAAATGCAAAAAAAAAGCCAACCGATATGGaa                                                   >  2:5185790/1‑150 (MQ=255)
                                                                                                    cACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAAAATAA‑AGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAAAATGAAAATTATTAGGGGGGTGTAAGGCAAAAAAACGGCCTCCCGGGAGGGaa                                                   >  2:6785948/1‑150 (MQ=255)
                                                                                                      cGTACATTTTCATATCCTGTAGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCGAAAAAACGGCCTCCCGATAGGGAAGc                                                 <  2:580053/150‑1 (MQ=255)
                                                                                                       ttACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGcc                                                <  2:4456800/149‑1 (MQ=255)
                                                                                                        tACATTTTCATTCCGTGTGGTTTTTTTTCGATTACAAAATATAT‑TTAATTTTTTTGTTTTTAATTATCGAGATACCGACCACAAAAACGACAATCTGAAAATTATTCGAGGAGTGAACGGCAAAAAAACGGCCTCCCGATAGGGAAGCCg                                               <  1:1986666/150‑1 (MQ=255)
                                                                                                          cATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAACCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                             <  2:3690975/150‑1 (MQ=255)
                                                                                                          cATTTTCATATCCTGTCGTTTGTTTACGATTTCAAAATATAA‑TTAACTTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                             <  1:5561894/150‑1 (MQ=255)
                                                                                                           aTTTTCATATCCTGTCTTTTGTTTTCTATTAAAAAATATAA‑TTAAATTTTTTTTTTTTAAATATCGATATAACGATAAAAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCTTCCCTATCGGGAGGCCGttt                                            <  1:3071130/150‑3 (MQ=255)
                                                                                                           aTTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAA‑TCAAATTATTTTTTTTTAAATATCGATATAACTATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                            <  1:8594088/150‑1 (MQ=255)
                                                                                                            ttttCATATCCTGTCGTGTTTTTTCGATTTCAAAATATAA‑TGAAATTATTTTTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGc                                           >  2:5481044/1‑150 (MQ=255)
                                                                                                             tttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGTTAACGATCACAAAAAAGACAATATGAAAATTATTCGTGGAGTGAAAGGGAAAAAAAAGGCCTCCCGATAGGGAAGCCGTAGCa                                          >  2:514123/1‑150 (MQ=255)
                                                                                                              ttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TTAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaa                                         <  2:8319398/150‑1 (MQ=255)
                                                                                                              ttCATATCCTGTCGTGTGTTTTCGAATTAAAAAAATAA‑TGAAAATATTTGTTTTTAAATAACGAGATAACGAACACAAAAACGACAATATGAAAATTAATAGAGGAGGGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaa                                         >  2:6150727/1‑150 (MQ=255)
                                                                                                               taatatCCTTTGGTTTTTTTTCTTTTAAAAAATATAT‑TTAAATTTTTTTTTTTTAAATATCGATATACCTATCAAAAAAACGACAATATTAAAATTATTCGCGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCaaa                                        <  1:2969913/148‑1 (MQ=255)
                                                                                                                 atatCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGt                                      <  1:1223578/150‑1 (MQ=255)
                                                                                                                     ccTGTCGTTTGTTTTCTATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                                  <  2:7750715/150‑1 (MQ=255)
                                                                                                                     ccTGTCGTTTGTTTTCG‑TTTTATAATTTAACTGAATTTATTTTTTTTTAAATATCTATATACCGATCACAAAAACGACAATTTGAAATTTATTCGAGGAGTGAATGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgcgc                                  <  1:8678855/150‑1 (MQ=255)
                                                                                                                       tGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGTAAGGCACAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgt                                >  2:6648100/1‑150 (MQ=255)
                                                                                                                         tcttttgtttTCTATTTCAAAATATAA‑TTAAATTTTTTTTTTTTAAATATCTAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtgt                              <  1:3685037/147‑1 (MQ=255)
                                                                                                                           gtttgtttTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt                            >  1:1137771/1‑150 (MQ=255)
                                                                                                                            tttgtttTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtttt                           >  1:252682/1‑150 (MQ=255)
                                                                                                                             ttgtttTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttttt                          <  2:4347361/150‑1 (MQ=255)
                                                                                                                               gtttTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                        >  1:6092823/1‑150 (MQ=255)
                                                                                                                                ttttCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCAC‑AAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                      >  1:2160938/1‑150 (MQ=255)
                                                                                                                                      aTTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGAt                 >  2:3524463/1‑150 (MQ=255)
                                                                                                                                        ttCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc               >  1:5738466/1‑150 (MQ=255)
                                                                                                                                         taaaaaTATAA‑TTAAATTATTTGTTTTTAAATATCGAGATACCGAACACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg              <  1:8050702/148‑1 (MQ=255)
                                                                                                                                         tCAAAATATAA‑TGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg              <  1:78209/150‑1 (MQ=255)
                                                                                                                                         tCAAAATATAA‑TGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg              <  2:7617126/150‑1 (MQ=255)
                                                                                                                                         tCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg              <  2:155113/150‑1 (MQ=255)
                                                                                                                                         tCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg              <  2:1851678/150‑1 (MQ=255)
                                                                                                                                         tCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg              <  2:1987746/150‑1 (MQ=255)
                                                                                                                                         acAAAATATAA‑TTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCg              <  1:2446653/149‑1 (MQ=255)
                                                                                                                                          aaaaaTATAA‑TGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCgg             <  1:6648788/149‑1 (MQ=255)
                                                                                                                                            aaaTATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGgtg           >  2:2795182/1‑150 (MQ=255)
                                                                                                                                            aaaTATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATGATTCGGGGAGTGAAAGGCAAAAAAACGGCCTCCCGAGAGGGAAGCCGTAGCAAAGTGCGCGGGGGTTTATGCCGGATGCGgtg           >  2:6327659/1‑150 (MQ=255)
                                                                                                                                              atatAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTa         >  1:4398304/1‑150 (MQ=255)
                                                                                                                                                ataA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTaaa       <  2:2771951/150‑1 (MQ=255)
                                                                                                                                                  aa‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACa     >  1:8133879/1‑150 (MQ=255)
                                                                                                                                                      gAAATTATTTTTTTTTAAATATCTATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc  <  2:7936606/150‑1 (MQ=255)
                                                                                                                                                      gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc  <  1:585082/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAA‑TGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAC  >  NC_000913/2177096‑2177394

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: ‑