Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,079,616 G→T 5.2% intergenic (‑286/+28) cmtB ← / ← tktA putative mannitol‑specific enzyme IIA component of PTS/transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,079,6160GT 5.2% 428.4 / 6.9 147intergenic (‑286/+28)cmtB/tktAputative mannitol‑specific enzyme IIA component of PTS/transketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (63/76);  new base T (2/6);  total (65/82)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 150 total aligned reads displayed.

TAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGC  >  NC_000913/3079470‑3079765
                                                                                                                                                    |                                                                                                                                                     
taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAGTTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAg                                                                                                                                                    <  1:841945/150‑1 (MQ=255)
 aGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGc                                                                                                                                                   <  2:62551/150‑1 (MQ=255)
   atatAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGa                                                                                                                                                 >  1:763424/1‑150 (MQ=255)
    tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGAc                                                                                                                                                <  1:1424437/150‑1 (MQ=255)
    tataACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGAc                                                                                                                                                >  1:1385013/1‑150 (MQ=255)
      taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCt                                                                                                                                              >  2:1705477/1‑150 (MQ=255)
        actacaaATGTCAATACTCATATCGTGATGCATATTTCATAATTGATATTTTT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCTTCGGCGACCttt                                                                                                                                            <  2:235733/147‑1 (MQ=255)
        acaacaAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCttt                                                                                                                                            >  2:1716890/1‑150 (MQ=255)
         caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCTTCGGCGACCtttt                                                                                                                                           <  1:401569/150‑1 (MQ=255)
           acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCTTCGGCGACCtttttt                                                                                                                                         <  2:1641810/150‑2 (MQ=255)
              aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTTATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTAccc                                                                                                                                      <  1:1112898/150‑1 (MQ=255)
              aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTAccc                                                                                                                                      <  1:726624/150‑1 (MQ=255)
              aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACAAGCATTATCTATTGCTTGTCAGAAAACGTACATATCAAGGTAA‑TAAAAAAGGTCGCAGAAGAGACCTATTTAccc                                                                                                                                      >  1:225657/1‑150 (MQ=255)
               aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCg                                                                                                                                     >  2:1095625/1‑150 (MQ=255)
               aTGTCAATAAGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCg                                                                                                                                     <  1:143130/150‑1 (MQ=255)
                 gTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTAAAACACGCCTTATCTATTGCTTGTCCGAAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGaa                                                                                                                                   >  1:113291/1‑150 (MQ=255)
                   cAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                                 <  1:435108/150‑1 (MQ=255)
                   cAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGGCGCCGAAGCGACCTTTTTACCCGAAAt                                                                                                                                 <  2:1032812/150‑1 (MQ=255)
                    aaTACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATg                                                                                                                                >  2:925171/1‑150 (MQ=255)
                      tACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCt                                                                                                                              >  2:1210251/1‑150 (MQ=255)
                         gCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAt                                                                                                                           <  2:1276682/150‑1 (MQ=255)
                          aatatCGTGATGAATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAAtt                                                                                                                          <  1:1086811/149‑1 (MQ=255)
                           atatCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTa                                                                                                                         <  1:1403125/150‑1 (MQ=255)
                             atCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCAGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAca                                                                                                                       >  2:1592239/1‑150 (MQ=255)
                              tCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTTTTATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcag                                                                                                                      <  2:187300/150‑1 (MQ=255)
                               cGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagc                                                                                                                     <  1:1020387/150‑1 (MQ=255)
                               cGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTAATGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagc                                                                                                                     <  2:244545/150‑1 (MQ=255)
                                gTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagca                                                                                                                    >  2:600008/1‑150 (MQ=255)
                                  gATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGt                                                                                                                  >  1:1664070/1‑150 (MQ=255)
                                  gATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGt                                                                                                                  <  2:763424/150‑1 (MQ=255)
                                     gCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGAAGTTCt                                                                                                               >  1:1685528/1‑150 (MQ=255)
                                      cATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCtt                                                                                                              <  1:709345/150‑1 (MQ=255)
                                        tatTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGACGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCtttt                                                                                                            >  1:322267/1‑150 (MQ=255)
                                         atTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCTGTTCTTTTg                                                                                                           <  1:1178031/150‑1 (MQ=255)
                                         atTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTg                                                                                                           <  2:1196370/150‑1 (MQ=255)
                                         atTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTg                                                                                                           <  2:385027/150‑1 (MQ=255)
                                          ttACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCAGCAAACGGACATATCAAGGTAATTAAAAAAGGTCGCAGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTg                                                                                                           <  2:1315378/150‑1 (MQ=255)
                                           tACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCt                                                                                                         <  1:18788/150‑1 (MQ=255)
                                            aCATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACCCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCtt                                                                                                        >  2:98278/1‑150 (MQ=255)
                                            aCATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCtt                                                                                                        <  2:1361382/150‑1 (MQ=255)
                                               tAAATGATATTTAT‑TGCATATAAAAATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGCTCGCCGAAGCGACCTTTTTACCCGAAATGCAAATTACAGCAGTTCTGTTGCTTTCg                                                                                                     <  1:1308019/150‑1 (MQ=255)
                                                aaTTGATTTTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACTGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACATTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGc                                                                                                    <  1:1612316/150‑1 (MQ=255)
                                                aaTTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGc                                                                                                    <  1:1542522/150‑1 (MQ=255)
                                                aaTTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGc                                                                                                    <  1:827281/150‑1 (MQ=255)
                                                 aTTGATATTTAT‑TGCATATAAATATTTGTGATCTACCACACGCCTTATCTATTGCTTGTCAGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTAGCCGAAGCTAACTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGca                                                                                                   >  1:111056/1‑150 (MQ=255)
                                                  ttGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGcaa                                                                                                  >  1:1603560/1‑150 (MQ=255)
                                                  caaatatTTATATGCA‑ATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGcaa                                                                                                  >  2:613143/4‑150 (MQ=255)
                                                      tatttat‑tGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGcatcaag                                                                                              >  2:246932/1‑149 (MQ=255)
                                                      tatttat‑tGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCAGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCaacaac                                                                                              >  1:603964/1‑150 (MQ=255)
                                                       atttat‑tGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACg                                                                                             >  2:489032/1‑150 (MQ=255)
                                                         ttat‑tGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGtt                                                                                           <  2:1124274/150‑1 (MQ=255)
                                                          tat‑tGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTa                                                                                          >  2:146001/1‑150 (MQ=255)
                                                            t‑tGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATc                                                                                        <  1:974043/150‑1 (MQ=255)
                                                                cATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAA‑GTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACa                                                                                    >  2:190547/1‑150 (MQ=255)
                                                                    taAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCCAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                                 >  1:1312176/1‑150 (MQ=255)
                                                                             gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                        >  2:352868/1‑150 (MQ=255)
                                                                             gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                        <  1:1070514/150‑1 (MQ=255)
                                                                             gtgATATACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                        <  1:218928/150‑1 (MQ=255)
                                                                               gATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActc                                                                      >  2:591516/1‑150 (MQ=255)
                                                                                  cTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                   <  2:1370182/150‑1 (MQ=255)
                                                                                  cTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                   >  2:98486/1‑150 (MQ=255)
                                                                                  cTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                   >  1:971687/1‑150 (MQ=255)
                                                                                  cTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                   >  1:964504/1‑150 (MQ=255)
                                                                                  cTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                   >  1:758140/1‑150 (MQ=255)
                                                                                   tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCa                                                                  >  2:907773/1‑150 (MQ=255)
                                                                                   tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCa                                                                  >  1:881904/1‑150 (MQ=255)
                                                                                    acaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaa                                                                 >  2:867372/1‑150 (MQ=255)
                                                                                      aacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAc                                                               <  1:1308131/150‑1 (MQ=255)
                                                                                      aacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAc                                                               >  1:920805/1‑150 (MQ=255)
                                                                                       acacGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAAAGTTATCAACAGTGAAGCCGAACTCTTCAAAca                                                              >  2:701880/1‑150 (MQ=255)
                                                                                        cacGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcag                                                             <  2:760199/150‑1 (MQ=255)
                                                                                         acGACTTATCTATTGCTTGTACTCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcagc                                                            <  1:147267/150‑1 (MQ=255)
                                                                                            ccTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACagcagc                                                         <  1:894921/150‑1 (MQ=255)
                                                                                                   tATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCg                                                  <  1:687446/150‑1 (MQ=255)
                                                                                                   tATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCg                                                  <  1:1552430/150‑1 (MQ=255)
                                                                                                    aTTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCgg                                                 <  2:743530/150‑1 (MQ=255)
                                                                                                       gCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGc                                              <  2:1385013/150‑1 (MQ=255)
                                                                                                       gCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGc                                              >  2:1010157/1‑150 (MQ=255)
                                                                                                        cTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                             <  1:591516/150‑1 (MQ=255)
                                                                                                        cTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGAGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                             >  2:1602606/1‑150 (MQ=255)
                                                                                                        cTTGTCCGCAAACGGACATATCAACGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                             <  1:321232/150‑1 (MQ=255)
                                                                                                              cGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCa                                       >  1:1238201/1‑150 (MQ=255)
                                                                                                               gCAAACGGACATATCAAGGTAA‑TAAAAAAGGTAGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACATTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCAc                                      <  2:1234475/150‑1 (MQ=255)
                                                                                                                cAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCAcc                                     <  1:925517/150‑1 (MQ=255)
                                                                                                                  aaCGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGa                                   >  1:919191/1‑150 (MQ=255)
                                                                                                                          tatCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATc                                                                                        >  1:650127/1‑89 (MQ=255)
                                                                                                                          tatCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATc                                                                                        <  2:650127/89‑1 (MQ=255)
                                                                                                                          tatCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTc                           >  1:1018995/1‑150 (MQ=255)
                                                                                                                          tatCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTc                           >  1:673460/1‑150 (MQ=255)
                                                                                                                                  taa‑taaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcgac                   <  2:881904/150‑1 (MQ=255)
                                                                                                                                   aa‑taaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcgacg                  <  1:122053/150‑1 (MQ=255)
                                                                                                                                      taaAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACAGACTAt                >  2:1655121/1‑150 (MQ=255)
                                                                                                                                        aaaaaGGTCGCTTCGGCGCCCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGACCTCTTCAAACAGCCGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATCCCGACGATAg              >  2:278047/1‑150 (MQ=255)
                                                                                                                                        aaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAg              >  2:43815/1‑150 (MQ=255)
                                                                                                                                           aaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCAc           <  2:1238201/150‑1 (MQ=255)
                                                                                                                                             ggTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCg         <  2:1310452/150‑1 (MQ=255)
                                                                                                                                              gTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGt        <  1:1716890/150‑1 (MQ=255)
                                                                                                                                              gTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGt        >  2:653446/1‑150 (MQ=255)
                                                                                                                                                    gAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGc  <  1:613143/150‑1 (MQ=255)
                                                                                                                                                    gAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGc  <  2:111056/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAT‑TGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAA‑TAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAGGC  >  NC_000913/3079470‑3079765

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑