Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,451,286 T→G 16.7% L145R (CTG→CGG)  thiE → thiamine‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,451,2860TG16.7% 375.5 / 18.5 217L145R (CTG→CGG) thiEthiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (109/64);  new base G (28/14);  total (139/78)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

GCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCA  >  NC_002947/5451137‑5451435
                                                                                                                                                     |                                                                                                                                                     
gCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCt                                                                                                                                                       >  1:2461436/1‑150 (MQ=255)
  tCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGc                                                                                                                                                     <  1:3132650/150‑1 (MQ=255)
    ggCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGcgc                                                                                                                                                   >  1:2006486/1‑148 (MQ=255)
    ggCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTc                                                                                                                                                   <  1:3946453/150‑1 (MQ=255)
    ggCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTc                                                                                                                                                   >  1:1162509/1‑150 (MQ=255)
      ccGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGa                                                                                                                                                 <  2:1999653/150‑1 (MQ=255)
       cGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAg                                                                                                                                                <  2:3978717/150‑1 (MQ=255)
         tCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGcc                                                                                                                                              >  1:5167080/1‑149 (MQ=255)
            ggCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGc                                                                                                                                           >  1:5552911/1‑150 (MQ=255)
             gCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCGGGCTCGACCAGGcc                                                                                                                                          >  1:3613724/1‑150 (MQ=255)
              cAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGccc                                                                                                                                         <  2:4549292/150‑1 (MQ=255)
              cAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCCGGccc                                                                                                                                         >  1:620848/1‑150 (MQ=255)
                atcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCcgc                                                                                                                                       >  2:3930564/1‑150 (MQ=255)
                 tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCACGCCcgcg                                                                                                                                      >  2:1459588/1‑150 (MQ=255)
                 tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGGCGGCCGGCTCGAGCAGGCCcgcg                                                                                                                                      >  2:428545/1‑150 (MQ=255)
                    tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGccc                                                                                                                                   >  2:311010/1‑150 (MQ=255)
                    tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTCGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGGCGGCGCGGGCCGGCCCCCGccc                                                                                                                                   >  1:3449454/1‑150 (MQ=255)
                      ggCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCcccccc                                                                                                                                 >  1:3092641/1‑148 (MQ=255)
                      ggCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCGCGccccc                                                                                                                                 >  1:2594978/1‑148 (MQ=255)
                       gCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCCACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAgg                                                                                                                                >  2:1274972/1‑150 (MQ=255)
                       gCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCacg                                                                                                                                >  1:3709751/1‑148 (MQ=255)
                          ccACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTgg                                                                                                                             >  2:3408120/1‑149 (MQ=255)
                          ccACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCacggga                                                                                                                             >  1:5377777/1‑145 (MQ=255)
                           cACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGaa                                                                                                                            <  2:1751576/150‑1 (MQ=255)
                           cACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCCGGTGaa                                                                                                                            >  2:2599403/1‑150 (MQ=255)
                                gCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGACACTGc                                                                                                                       >  2:3855558/1‑150 (MQ=255)
                                   aCGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGCAACTGCCGa                                                                                                                    >  2:4502775/1‑150 (MQ=255)
                                    cGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCCCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGAt                                                                                                                   <  2:4587229/150‑1 (MQ=255)
                                       cAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCCGGCCCGCGCCCCGGGGAGACTGCCGAgcgc                                                                                                                >  1:4181132/1‑150 (MQ=255)
                                          ccTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCCGGCCCGCGCCCCCCGGAAACAGCCGAGCGCCga                                                                                                             >  1:3042513/1‑149 (MQ=255)
                                           cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGTGAAACTGCCGATCGCCggg                                                                                                            >  1:3646918/1‑148 (MQ=255)
                                            tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGa                                                                                                           >  2:784085/1‑150 (MQ=255)
                                             ttcgcTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCACGCCCCCGCCCAGGTGAAACTGCCGATCGCCGTGAt                                                                                                          >  2:4155612/5‑150 (MQ=255)
                                              gAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGCAACTGCCGATCGCCGTGATc                                                                                                         >  1:4136974/1‑150 (MQ=255)
                                               aaCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCGATTCTGTCACCAAGCCGGGCGCGCCCGCTGCAAACGTCGGCCGGCGCGAGCAGGCCCCCGCCCAGGGGAAACTGCCGCTCGCCGTGATCg                                                                                                        >  2:1274200/1‑150 (MQ=255)
                                               aaCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCCCCCAGGTGAAACTGCCGATCGCCGTGATCg                                                                                                        >  1:372551/1‑150 (MQ=255)
                                                  gggCCGCCCCGGCTGCCTGCGATGGCGCCTGCTACTTGGCCTTTGGCCGCTTCTTCATGCCTGTCACCACGCGGGGCGCGCCCGCTGCCAAGGTCGGCGTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTg                                                                                                     <  1:1374843/149‑1 (MQ=255)
                                                        cccAGGCGGCCCGCGAAGGCGCCAGCTACGGGGCCTTGGGCGGCCTCTCCCCTCCGGTCCCCACGCGGGGGGCGCCCGCCGCCACCGCCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCa                                                                                               <  2:2511721/150‑1 (MQ=255)
                                                          cAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCAcc                                                                                             >  2:2410050/1‑150 (MQ=255)
                                                           aGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGGGAAACTGCCGAGCGCCGTGGTCGGTGGCATCAccc                                                                                            >  1:2531299/1‑150 (MQ=255)
                                                             gCTGCCAGCGACGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCCCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTc                                                                                          <  2:1675312/150‑1 (MQ=255)
                                                               tGCCAGCGACGGCGCCAGCTACGTGGCCTTTGGCCGCTTCGGCACTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGCTCGCCGTGAGCGGTGGCATCACCCTCGa                                                                                        >  2:4073626/1‑150 (MQ=255)
                                                               tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGa                                                                                        >  2:5582852/1‑150 (MQ=255)
                                                               ggCCCGCGATGGCGCCAGCTAGGTGGCCTTTGGCCGCTTCTTCAATGCCGCCCCCAAGCCGGGCGCGCCCGCTGCCACGGCCGGCCGGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGa                                                                                        <  2:473381/149‑1 (MQ=255)
                                                                 ccAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCCGCTGAAACAGCCGCTCGCCGTGATCGGTGGCCTCACCCTCGACa                                                                                      >  1:5474395/1‑150 (MQ=255)
                                                                   aGCGACGGCGCCAGCTACGTGGCCTTTGGCCGATTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAc                                                                                    <  1:5139353/150‑1 (MQ=255)
                                                                   aGCGAAGGCGCCAGCTAGGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAc                                                                                    <  2:1169686/150‑1 (MQ=255)
                                                                   aGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATAGGTGGCATCACCCTCGACAAc                                                                                    >  2:1660411/1‑150 (MQ=255)
                                                                   aGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGATTCTTCATTTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAc                                                                                    <  2:1169664/150‑1 (MQ=255)
                                                                     cGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCa                                                                                                                                  <  2:2250997/102‑1 (MQ=255)
                                                                              cAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGTGCAGGCCCGCGCCCACGTGAAACAGCCGATCGCCCGGGTCGGTGGCATCACCCTCGACAACCCCGCCCCAca                                                                         >  1:2099780/1‑149 (MQ=255)
                                                                              cAGCTACGTGGCCATGGGCCGCTTCTCCAATTCTGTCCCCAACCCGGGCGCGCCCGCTGCCACCGGCCGCCGGCCCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACt                                                                         <  1:2251480/150‑1 (MQ=255)
                                                                                  tACGGGGCCTTTGGCCGCTTCTTCAATTCGGTCCCCAAGCCGGGCGCGCCCGCTGCCACCGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTc                                                                     <  2:332100/150‑1 (MQ=255)
                                                                                    cGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACAGGTCGc                                                                   >  2:2794998/1‑150 (MQ=255)
                                                                                        gCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCACGCCCCCGCCCCGGTGAAACTGCCCATCGCCGGGATAGGTGGCATCACCCACGACAACGCCCCCCCACAGGTCGCCCAc                                                               >  1:2197917/1‑150 (MQ=255)
                                                                                           tttGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGgc                                                            >  2:1375558/1‑150 (MQ=255)
                                                                                           tttGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACAGGTCGCCCACGgc                                                            >  2:5414382/1‑150 (MQ=255)
                                                                                           tatcgCCGATTCTACAATTCTGTCACCAAGCCGGGCGCGCACGCAGCCAAAGTCGGCATGCACGAGCAGGCCCCCCCCCAGGAGAAACTGCCGAACGCCCTGAGCCCTGGCATCAACATCGACAACCCCGCCCCACACTTCGCCCACCgc                                                            >  2:3780804/5‑150 (MQ=255)
                                                                                                ccGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGa                                                       >  2:1636609/1‑150 (MQ=255)
                                                                                                ccGCTTCTTCAATTCTGACACCAAGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGa                                                       <  1:3742589/150‑1 (MQ=255)
                                                                                                ccGCTTATTCAATTCTGTCCCCCCGCCGGGCGCGCCCGCGGCCACCGCCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGa                                                       <  1:3812407/150‑1 (MQ=255)
                                                                                                 cGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGTGAAAATGCCGATCGCCGTGATCGGTGGCAACACCCACGACAACGCCGCCCCACACGTCGCCCACGGCGCCCAc                                                      >  1:2326513/1‑150 (MQ=255)
                                                                                                   cttcttCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCGTGCTCGTGCCGGCGCGCGCCCAGGTGAAACTGCCGTTGGCCGTGATCGGTGGCCTCACCCTCGACAACGCCGCCCCACTGGTGGCCCACGGCGCCGACct                                                    >  2:5283317/1‑150 (MQ=255)
                                                                                                       ttcAATTCTGTCACCCCGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                <  2:3307075/150‑1 (MQ=255)
                                                                                                       ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTTCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                >  2:4089216/1‑150 (MQ=255)
                                                                                                       ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGTGAAACAGCCGATCGCCGTGATCGGTGGCCTCACCCTCCACAACGCCCCCCCCCAGGACGCCCACCGCGCCGACctcctg                                                >  1:728449/1‑150 (MQ=255)
                                                                                                       ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGGATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                >  2:4314710/1‑150 (MQ=255)
                                                                                                       ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                >  2:4681175/1‑150 (MQ=255)
                                                                                                       ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCCCCCCACTGGTCGCCCACGGCGCCGACctgcag                                                >  2:4934543/1‑148 (MQ=255)
                                                                                                        tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTgg                                               >  2:2069391/1‑150 (MQ=255)
                                                                                                          aaTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCg                                             >  2:2144269/1‑150 (MQ=255)
                                                                                                          aaTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCGCGACAACGCCGCCCCACAGGTCGCCCACGGCGCCGACCAGCAGGCg                                             >  1:1782599/1‑150 (MQ=255)
                                                                                                            ttCTGTCACCCAGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGt                                           <  2:4985648/150‑1 (MQ=255)
                                                                                                               tGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGcc                                                                                                               >  1:879331/1‑79 (MQ=255)
                                                                                                               tGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGcc                                                                                                               <  2:879331/79‑1 (MQ=255)
                                                                                                                     cacgCCGGGCGCGCCCGCTGCCACCGCCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  <  2:2188630/147‑1 (MQ=255)
                                                                                                                     cAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  >  2:4920744/1‑150 (MQ=255)
                                                                                                                     cAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  >  2:3995236/1‑150 (MQ=255)
                                                                                                                     cAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCCGGTGAAACTGCCGATCGCCGTGATCGGTGGCCTCACCCTCGACAACCCCCCCCCACAGGTCGCCCACGGCGCCGCCCTGCGGGCGGTGATCCACgg                                  >  1:4845371/1‑150 (MQ=255)
                                                                                                                            ggCGCGCCCGCGGCCCCCGGCGGCCGGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTc                           <  1:2221124/150‑1 (MQ=255)
                                                                                                                            ggCGCGCCCGCGGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTc                           <  2:3659533/150‑1 (MQ=255)
                                                                                                                              cgcgcCCGCTGCCAACGGCGGCGTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCgg                         <  1:101344/150‑1 (MQ=255)
                                                                                                                               gcgcCGGCGGCCACGGTGGGCCGGGGCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGt                        <  1:4310374/150‑1 (MQ=255)
                                                                                                                               gcgcCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGt                        <  2:3182884/150‑1 (MQ=255)
                                                                                                                               gcgcCCGCGGCCCCCGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGt                        <  2:4808043/150‑1 (MQ=255)
                                                                                                                                  cccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGcc                     >  2:5448196/1‑150 (MQ=255)
                                                                                                                                  cccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGcc                     >  2:5390048/1‑150 (MQ=255)
                                                                                                                                   ccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCg                    >  1:2914774/1‑150 (MQ=255)
                                                                                                                                   ccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCg                    >  2:5249618/1‑150 (MQ=255)
                                                                                                                                     gCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGAc                  >  1:4960607/1‑150 (MQ=255)
                                                                                                                                          cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc             >  1:2297336/1‑150 (MQ=255)
                                                                                                                                           aaCGTCGGCCTGCTCGAGCAGGCACGCGCCCAGGTGAAACTGCTGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgcg            >  2:4367615/1‑150 (MQ=255)
                                                                                                                                             cGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCa          <  2:1568640/150‑1 (MQ=255)
                                                                                                                                              gcgggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAg         >  1:5097131/4‑150 (MQ=255)
                                                                                                                                                 ggCCTGCCCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGaa      <  1:2146487/150‑1 (MQ=255)
                                                                                                                                                  gcgtGCGCGAGCTGGCCCGCGCCCAGGTGATACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGTCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAg     <  2:3613724/147‑1 (MQ=255)
                                                                                                                                                    cTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGACGTc   >  2:5071630/1‑150 (MQ=255)
                                                                                                                                                    cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTc   >  2:1813029/1‑150 (MQ=255)
                                                                                                                                                    cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTc   >  1:4401017/1‑150 (MQ=255)
                                                                                                                                                     tGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCa  >  1:3588107/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCA  >  NC_002947/5451137‑5451435

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: