Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,448,259 T→G 11.9% 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,448,2590TG11.9% 279.4 / 2.2 144L145R (CTG→CGG) thiEthiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (72/45);  new base G (14/7);  total (89/55)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCA  >  NC_002947/5448111‑5448408
                                                                                                                                                    |                                                                                                                                                      
 tCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGc                                                                                                                                                      <  2:2663606/150‑1 (MQ=255)
  cGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCGACGTGGCCTTTGGCCGCTTCTTCAATGCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTcggccggcg                                                                                                                                                     >  1:394675/1‑149 (MQ=255)
   ggCCGTCTGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTc                                                                                                                                                    <  1:3086109/150‑1 (MQ=255)
      cGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAg                                                                                                                                                 <  1:1286195/150‑1 (MQ=255)
      cGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAg                                                                                                                                                 <  2:4982116/150‑1 (MQ=255)
      cGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAg                                                                                                                                                 >  2:2931739/1‑150 (MQ=255)
       gTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGc                                                                                                                                                <  1:4126520/150‑1 (MQ=255)
       gTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGc                                                                                                                                                <  2:1829525/150‑1 (MQ=255)
       gTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGc                                                                                                                                                <  2:1829540/150‑1 (MQ=255)
       gTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCCGCTCGAGc                                                                                                                                                <  1:3884880/150‑1 (MQ=255)
         cAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAg                                                                                                                                              <  1:3647414/150‑1 (MQ=255)
         cAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACGGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAg                                                                                                                                              <  1:4963849/150‑1 (MQ=255)
           ggCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGc                                                                                                                                            <  2:4023281/150‑1 (MQ=255)
           ggCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCGCGAGCAGGc                                                                                                                                            >  2:1298350/1‑150 (MQ=255)
            gCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAg                                                                                                                                                 >  1:4679901/1‑144 (MQ=255)
            gCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAg                                                                                                                                                 <  2:4679901/144‑1 (MQ=255)
              aaTCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCTCGGGCAGGcccc                                                                                                                                        >  1:3533502/1‑149 (MQ=255)
                tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCcgcg                                                                                                                                      >  2:2040168/1‑150 (MQ=255)
                  atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGcc                                                                                                                                    >  2:2973136/1‑150 (MQ=255)
                   tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTCCTCGAGCCGGCCCGCGccc                                                                                                                                   >  2:4315027/1‑150 (MQ=255)
                        tCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTg                                                                                                                              >  2:1558920/1‑150 (MQ=255)
                          cACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTTGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGGCGGGACGAGCAGCGCCGCGCCCAGCTGaa                                                                                                                            >  1:1870192/1‑150 (MQ=255)
                          cACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCTCGAGCCGGCCCGCGCCCAGGTGaa                                                                                                                            >  2:1134856/1‑150 (MQ=255)
                           aCCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCGCGCCCAGGGGaaa                                                                                                                           >  1:2982539/1‑150 (MQ=255)
                               gCCACGCCAGCCTCGAACTGGCCGCCCCGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGc                                                                                                                       <  2:4989832/150‑1 (MQ=255)
                               gCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGTCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGACTGCTCGAGCCGGCCCCCCCCCAGGTGAAACttg                                                                                                                       >  2:2919977/1‑148 (MQ=255)
                                   cGCCAGCCTCGACCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGAt                                                                                                                   <  2:297022/150‑1 (MQ=255)
                                    gCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCCCCCCGGTGAAACTGCCGaac                                                                                                                  >  1:151640/1‑148 (MQ=255)
                                         ccTCGAACTGGCCGCCCAGGCTGCCAGCGACGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCCCGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGt                                                                                                             <  2:377410/150‑1 (MQ=255)
                                          cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCCCGCCCACGTGAAACTGCCGATCGCCGTg                                                                                                            >  1:2762132/1‑150 (MQ=255)
                                           tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGa                                                                                                           >  2:2466969/1‑150 (MQ=255)
                                           tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCCGGCCCGCGCCCCCGGGAAACAGCCGAACGCCCTGa                                                                                                           >  1:1809286/1‑150 (MQ=255)
                                              aaCTGGCCGCCCTGGCTGCCAGCGATGGCGCCAGCTAGGTGGCCTTGGGCCGCTTCTTCATTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCg                                                                                                        <  1:1544595/150‑1 (MQ=255)
                                              aaCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTGGCCCGCTTCTTCAATTCTGTCCCCACGCCGGGCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCg                                                                                                        <  2:920773/150‑1 (MQ=255)
                                                  ggCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTgg                                                                                                    >  2:3069898/1‑150 (MQ=255)
                                                     cgccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCGCGCCCCCGGGAAACTGCCCAAAGCCGGGAGCGGTGGCAt                                                                                                 >  1:2194905/1‑150 (MQ=255)
                                                      gccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCACGCCCGCGCCCACCGGAAACAGCCCATCGCCGTGATCGGGGGCATc                                                                                                >  1:4375521/1‑150 (MQ=255)
                                                                   gCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTGTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGGGGGGTCGAGCAGGCCCGCGCCCAGGTGAAAATGCCGGTCGCCGTGATAGGTGGGTTACGCCTCGACaaag                                                                                   >  1:4281931/1‑148 (MQ=255)
                                                                     gAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcgcc                                                                                 >  2:2116810/1‑150 (MQ=255)
                                                                     gAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGACGATCGCCGAGATCGGTGGCATCACCCTCGACAAcgcc                                                                                 >  2:2712586/1‑150 (MQ=255)
                                                                      aaGGCGCCAGCTACGTGGCCTTTGGCCGCTTCCTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCCCGCCCCCGGGCAACTGCCGCTCGCCGGGGTCGGTGGCCGCACCCGCCACAAccccc                                                                                >  1:2800264/1‑149 (MQ=255)
                                                                       aGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcgccgc                                                                               >  2:4356951/1‑150 (MQ=255)
                                                                         gcgcCAGCTCCGGGGCCTTTGGCGGCTTCGTCCTTTCTGTCACCACGCCGGGCGCGCCCGCTGCCCACGTCGGCCTGCGCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGccc                                                                             <  2:3193978/150‑1 (MQ=255)
                                                                         gcgcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGccc                                                                             >  2:4552200/1‑150 (MQ=255)
                                                                         gcgcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCCCGCAGGCCCGCCCCCCCGGGAAACAGCCCCTCGCCCGGAGCGGGGGCCTCACCCACCACACcccccccc                                                                             >  1:2282320/1‑150 (MQ=255)
                                                                          cgcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGcccc                                                                            >  2:1492498/1‑150 (MQ=255)
                                                                           gcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATAAACCTCGACAACGCCGccccc                                                                           >  1:1665042/1‑149 (MQ=255)
                                                                            ccAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGAGCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGcccccc                                                                          >  1:2058532/1‑148 (MQ=255)
                                                                              agctcggTGGCATTTGGCCGTTTATTAACTTCTGACACCCCCCGGGGCGCCCCCGGTGCCAACGTGGGCCGGCTCGCGCTGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTg                                                                        <  1:1090707/144‑1 (MQ=255)
                                                                               gCTAGGTGGCTTTTGGCGGCTTTTTCATTCCTGACCCCACGCGGGGCGCGCCCGCTGCCACCGTGGGCCGGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTgg                                                                       <  1:3104748/150‑1 (MQ=255)
                                                                               gCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACACCGCCGCCCCACTgg                                                                       >  2:1036006/1‑150 (MQ=255)
                                                                                      ggCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTCCTCGAGCCGGCCCGCGCCCACGTGAAACTGCCGATAGCCGTGAGCGGTGGCATCACCCTCGACAACCCCGCCCCACTGGTCGcccc                                                                >  2:3938758/1‑149 (MQ=255)
                                                                                         cTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCCACACCCTCGACAACGCCGCCCCACTGGTCGCCCAccg                                                             >  2:210708/1‑148 (MQ=255)
                                                                                         cTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCACGCCCCCGCCCACGAGAAACAGCCCATCGCCCTGATAGGTGGCATAACCCACCACAAACCCCCCCCACACGTCGCCCACgc                                                             >  1:878435/1‑149 (MQ=255)
                                                                                         cTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTCCCAACGACGGCCTGCACGACACACCACGCGCCCACCTTAAAATGCCCCTAGCCGTGATCCGTGGAATAACCAGCGACAACCCCGCCCCACGGGTCGCCCACgg                                                             >  1:809252/1‑150 (MQ=255)
                                                                                               ccGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGa                                                       <  1:2931739/150‑1 (MQ=255)
                                                                                                cGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAc                                                      <  2:3304522/150‑1 (MQ=255)
                                                                                                cGCTACTTCAATGCTCTCCCCCCGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAc                                                      <  2:2217464/150‑1 (MQ=255)
                                                                                                   ttcttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCAACCTCGACAACCCCCCCCCCCAGGTCGCCCACGGCGCCGAccgg                                                   >  1:73627/1‑148 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCCACAACCCCCCCCCACAGGTCGCCCACGGCGCCGACcagcag                                                >  1:2820047/1‑148 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCCCCCACGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCGCGACAACGCCCCCCCACAGGTCGCCCACCGCGCCGACcagctg                                                >  1:3458978/1‑150 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGGGAGCGGTGGCATCACCCACGACAACGCCGCCCCACAGGGCGCCCACCGCGCCGACccgctg                                                >  1:4905561/1‑150 (MQ=255)
                                                                                                       tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACAGGTCGCCCACGGCGCCGACCAGCTgg                                               >  1:2548289/1‑150 (MQ=255)
                                                                                                          aTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGTGAAACTGCCGAGCGCCGGGATCGGGGGCATCACCCACGACAACGCCGCCCCCCTGGTCGCCCACCGCGCCGCCCAgcgggcgg                                            >  1:4876593/1‑150 (MQ=255)
                                                                                                           ttCTCTCACCACGCCGGGCGCGCCCGCTGCAACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGt                                           <  2:472960/150‑1 (MQ=255)
                                                                                                              tGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGAt                                        >  1:4635057/1‑150 (MQ=255)
                                                                                                              tGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGAt                                        >  2:3300088/1‑150 (MQ=255)
                                                                                                                  cccccgCCGGGCGCGCCCGCGGCCACCGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                    <  2:2428108/145‑1 (MQ=255)
                                                                                                                  aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGGGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                    >  1:2496023/1‑150 (MQ=255)
                                                                                                                    cccgCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  <  1:3377241/147‑1 (MQ=255)
                                                                                                                    cacgCCGGGCGCCCCCGGCGCCACCGTCGGCCGGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  <  1:804715/147‑1 (MQ=255)
                                                                                                                    cAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  >  2:3569448/1‑150 (MQ=255)
                                                                                                                     ccggcGGGCGCCCCCGCCCCCTACGTGGGCCGGCTCGCGCTGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGc                                 <  1:4807990/146‑1 (MQ=255)
                                                                                                                     ccgcggggCGCCCCCGCGCCCCCCGCCGGCCTGCTCGCCCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGc                                 <  1:4808020/145‑1 (MQ=255)
                                                                                                                     aaGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCACGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACAGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGc                                 >  2:757288/1‑150 (MQ=255)
                                                                                                                       gCCGGGCGCGCCCGCTGCCACCGACGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCtt                               <  2:2982539/150‑1 (MQ=255)
                                                                                                                       cccGGGCCCGCCCGCGCCCCCCGGCGGCCGGCTCGCCCCGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCtt                               <  1:1799999/149‑1 (MQ=255)
                                                                                                                          gggCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGtt                            <  2:2035065/150‑1 (MQ=255)
                                                                                                                           ggCGCGCCCGCTGCCACCGTCGGCCTGCCCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTc                           <  2:1700800/150‑1 (MQ=255)
                                                                                                                           ggCGCGCCCGCTGCCAACGTCGGCCTGCTCGAACAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTc                           <  2:3652964/150‑1 (MQ=255)
                                                                                                                                gcCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCTCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGc                      >  1:4855572/1‑150 (MQ=255)
                                                                                                                                 cccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGcc                     >  1:5069011/1‑150 (MQ=255)
                                                                                                                                 cccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGcc                     >  1:5069018/1‑150 (MQ=255)
                                                                                                                                     ctcccAACGCCGGCCACCTCCCCCAGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACACCGCCGCCCCACTGGTCGCCAACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACa                 <  1:1399660/147‑1 (MQ=255)
                                                                                                                                     cgcccACCGTCGGCCTGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACa                 <  2:2188788/147‑1 (MQ=255)
                                                                                                                                     cTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACa                 >  1:2690890/1‑150 (MQ=255)
                                                                                                                                     cTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATAGGTGGCATCACCCTCGACAACGCCGCCCCACAGGTCGCCCACCGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACa                 >  1:2690893/1‑150 (MQ=255)
                                                                                                                                        ccAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcg              >  2:134029/1‑150 (MQ=255)
                                                                                                                                           ccggcGGCCTGCGCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgcgc           <  1:1604358/146‑1 (MQ=255)
                                                                                                                                           aCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgcgc           >  1:406659/1‑150 (MQ=255)
                                                                                                                                               gggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGa       >  2:4513773/2‑150 (MQ=255)
                                                                                                                                                ggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGaa      <  1:2973136/150‑1 (MQ=255)
                                                                                                                                                 gCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAg     >  2:5071817/1‑150 (MQ=255)
                                                                                                                                                 gCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAg     >  2:1410587/1‑150 (MQ=255)
                                                                                                                                                   cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATTACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTc   >  2:1066602/1‑150 (MQ=255)
                                                                                                                                                   cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTc   >  2:1066717/1‑150 (MQ=255)
                                                                                                                                                   cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTc   >  1:3353975/1‑150 (MQ=255)
                                                                                                                                                   cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTc   >  2:5112960/1‑150 (MQ=255)
                                                                                                                                                    tGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCa  <  2:1749733/150‑1 (MQ=255)
                                                                                                                                                    tGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCa  <  2:3825944/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCA  >  NC_002947/5448111‑5448408

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: