Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,451,295 A→C 15.5% Q148P (CAG→CCG)  thiE → thiamine‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,451,2950AC15.5% 614.0 / 23.3 317Q148P (CAG→CCG) thiEthiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (156/100);  new base C (27/31);  total (184/133)
Fisher's exact test for biased strand distribution p-value = 5.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 319 total aligned reads displayed.

TCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGCCGC  >  NC_002947/5451146‑5451443
                                                                                                                                                     |                                                                                                                                                    
tCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCa                                                                                                                                                      >  2:2552711/1‑150 (MQ=255)
  aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAgg                                                                                                                                                    >  2:368908/1‑150 (MQ=255)
  aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAgg                                                                                                                                                    <  1:8292411/150‑1 (MQ=255)
  aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAgg                                                                                                                                                    >  1:8812849/1‑150 (MQ=255)
   ggCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTAGAGCAGGc                                                                                                                                                   >  2:1617996/1‑150 (MQ=255)
    gCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCCGGCGCGCCCGCTGCCAACCTCGGCCAGAACGAGCACGcc                                                                                                                                                  >  1:3161071/1‑150 (MQ=255)
    gCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTGTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCGCCAGCAGGcc                                                                                                                                                  >  2:3985668/1‑150 (MQ=255)
     cAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGccc                                                                                                                                                 <  2:164325/150‑1 (MQ=255)
      aaTCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCCACGTCGGCCAGCGCGAGCCGGCCcg                                                                                                                                                >  2:5479551/1‑150 (MQ=255)
      aaTCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGcccc                                                                                                                                                >  2:7984235/1‑149 (MQ=255)
        tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCcgcg                                                                                                                                              <  2:644243/150‑1 (MQ=255)
        tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCcgcg                                                                                                                                              <  1:6240066/150‑1 (MQ=255)
        tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCACGAGCACGCCcgcg                                                                                                                                              >  2:5833292/1‑150 (MQ=255)
          atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGcc                                                                                                                                            <  2:4336475/150‑1 (MQ=255)
          atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGcc                                                                                                                                            >  1:1449470/1‑150 (MQ=255)
           tcGGCTCCCCCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCGGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGccc                                                                                                                                           <  2:4413354/150‑1 (MQ=255)
             ggCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCGCGAGCAGGCCCCCGccccg                                                                                                                                         >  1:7014552/1‑148 (MQ=255)
              gCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGGGCAGGCCCGCccccccg                                                                                                                                        >  1:5985202/1‑147 (MQ=255)
               cTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCCGGCCCCCGcccccgg                                                                                                                                       >  1:5489278/1‑146 (MQ=255)
               cTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGATACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGt                                                                                                                                       <  2:1444194/150‑1 (MQ=255)
                tCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAgggg                                                                                                                                      >  2:8055194/1‑148 (MQ=255)
                tCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCCCGCCCCGGTg                                                                                                                                      >  1:4942609/1‑150 (MQ=255)
                tCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGCGCCGGCCCGCGCCCCgggg                                                                                                                                      >  1:6334341/1‑148 (MQ=255)
                   aCCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCCGCGACGGCGCCAGCTACGGGGCCTTTGGCGGCTTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGCCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGaaa                                                                                                                                   <  2:1749547/150‑1 (MQ=255)
                   aCCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGACGGCGCCAGCTACGTGGCCTTTGGCCGATTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGaaa                                                                                                                                   <  2:2531545/150‑1 (MQ=255)
                   aCCTGCCACGCCAGCCTCGAACGGGCCGCCCTGGCTGCCAGCGACGGCGCCAGATAGGTGGCCTTTGGCCGCTTCTTCAATTCTGTCCCCACGCCGGGCGCGCCCGCTGCCAACGTGGGCCTGCTCGGGCAGGCCCGCGCCCAGGTGaaa                                                                                                                                   <  1:2465124/150‑1 (MQ=255)
                       gCCACGCCAGCCTCGACCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGATTATTCAATTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGc                                                                                                                               <  1:5859456/150‑1 (MQ=255)
                         cACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGCGCCGGCCCGCGCCCCGGTGAAACTGCCg                                                                                                                             >  1:4877369/1‑150 (MQ=255)
                          aCGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCAGGTGAAACTGCCgc                                                                                                                            >  2:3112036/1‑149 (MQ=255)
                          aCGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCCCCCACGTGAAAATGCCGa                                                                                                                            >  1:2697476/1‑150 (MQ=255)
                              cAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTACTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCTCGAGCAGGCCCGCGCCCACGTGAAACTGCCGATCGc                                                                                                                        >  2:6575960/1‑150 (MQ=255)
                                 ccTCGAACTGGCCGCCCAGGCTGCCAGCGACGGCGCCTGCGACGTGGCCTTTGGCCGCTTCTTCACTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGt                                                                                                                     <  2:6368976/150‑1 (MQ=255)
                                   tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGa                                                                                                                   >  2:5934105/1‑150 (MQ=255)
                                   tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCCCCCCGGGGAAACTGCCGGACGCCGTGa                                                                                                                   >  1:2342309/1‑150 (MQ=255)
                                   tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCGCCCCCACGTGAAACTGCCGATCGCCGGGa                                                                                                                   >  1:1687977/1‑150 (MQ=255)
                                    cGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGGCGGCTGGAGCAGGCCCCCGCCCCGGGGAAACTGCCGATCGCCGGGAt                                                                                                                  >  1:8231935/1‑150 (MQ=255)
                                      aaCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCCGGCCCCCGCCCAGGGGCAACTGCCCATCGcccggatag                                                                                                                >  1:770633/1‑143 (MQ=255)
                                      aaCTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCCGGCAGGCCCGCGCCCACCTGAAACTGCCCAGCGCCGGGAGCg                                                                                                                >  1:1806705/1‑150 (MQ=255)
                                        cgggCCGCCCCGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGt                                                                                                              <  2:6897932/148‑1 (MQ=255)
                                        cgggCCGCCCAGGCTGCCCGCTACGGCGCCCGCTACGTGGCCTTTGCCCGCTCCTTCACTTCTGCCCCCAAGCCGGGCGCGCCCGCCCCCACCGCCGGCCCGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGt                                                                                                              <  1:3437201/148‑1 (MQ=255)
                                        cTGGCCGCCCAGGCTGTCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACCTCGGCCTGCGCGAGCACGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGt                                                                                                              >  2:4529741/1‑150 (MQ=255)
                                           gccgccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGACCACGCCCGCGCCCCCGTGAAACTGCCGAGAGCCCTGATCGGTggg                                                                                                           >  1:3756666/1‑149 (MQ=255)
                                              gccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGGTCGGTGGCATc                                                                                                        >  2:5023249/1‑150 (MQ=255)
                                                   ggCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGCGCAGGCCCGCGCCCCGGTGTCACTGCCGCTCCCCGGGGTAGGGGGCCTCAccca                                                                                                   >  1:7572398/1‑149 (MQ=255)
                                                    gCTGCCCGCGACGGCGCCCGCTCCGTGGCCTTTGCCCGCTTCATCACTCCTGTCCCCACGCCGGGCGCGCCCGCTGCCAACGGCGGCCGGCTCGAGCAGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTc                                                                                                  <  1:2662518/150‑1 (MQ=255)
                                                      tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGa                                                                                                >  2:8178104/1‑150 (MQ=255)
                                                      tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATAGCCGTGGTCGGTGGCATCACCCGCGa                                                                                                >  1:674855/1‑150 (MQ=255)
                                                       gCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCACGGGGAACTGCCCATCGCCGTGATCGGTGGCATCACCCTCCAc                                                                                               >  1:6601589/1‑150 (MQ=255)
                                                        ccAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGAACGCCGTGATCGGTGGCATCACCCACGACa                                                                                              >  1:640802/1‑150 (MQ=255)
                                                              acggCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCACTTCTGTCACCACGCCGGGCGCGCCCGCGGCCACCGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcgccg                                                                                        <  2:6094565/148‑1 (MQ=255)
                                                                 gcgcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCACGTGAAACGGCCGCTCGCCGTGATAGGTGGCATCACCCTCGACAACGCCGccc                                                                                     >  1:1323574/1‑150 (MQ=255)
                                                                 gcgcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGccc                                                                                     >  1:8002049/1‑150 (MQ=255)
                                                                 gcgcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCCCGCCCACGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGcccccc                                                                                     >  1:2014495/1‑150 (MQ=255)
                                                                  cgcCAGCTATGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGACGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGcccc                                                                                    <  2:2786246/150‑1 (MQ=255)
                                                                         tACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTc                                                                             <  2:1310937/150‑1 (MQ=255)
                                                                           cGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCACGTGAAACTGCCGATAGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACAGGgcgc                                                                           >  1:7314276/1‑150 (MQ=255)
                                                                              ggCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCAGGTGAAAAAGCCCAACGCCGGGATCGGTGGCAAAACCCACCAAAAACCCCCCCCACACGTCGcccc                                                                        >  1:2560467/1‑149 (MQ=255)
                                                                                        ctcttcttCAATCCGGCCCCCCCGCCGGGCGCGCCCGCGGCCAACGCGGGCCGGCTCGCGCAGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAc                                                              <  1:572984/148‑1 (MQ=255)
                                                                                        cGCTTCTTCAATTCTGTCCCCACGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAc                                                              <  2:3458284/150‑1 (MQ=255)
                                                                                        cGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAc                                                              >  2:1079464/1‑150 (MQ=255)
                                                                                        cGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCAGGTGAAACTGCCGATCGCCGAGATCGGTTGCATCACCCACGACAACCCCGCCCCACAGGTCGCCCACGGCGCCCAc                                                              >  1:4900664/1‑150 (MQ=255)
                                                                                        cGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCGGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAc                                                              <  2:8334647/150‑1 (MQ=255)
                                                                                         gCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGCCCCCCGCCCAGCTGAAAAATCCGAAAGCCGTGATAGGTGGCAAAACCCACCACAACCCCCCCCCACACCTCGCCCACCGCGCCGAcc                                                             >  1:5624359/1‑150 (MQ=255)
                                                                                           ttcttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAAAAGCCGATCGCCGTGATAGGTGGCAACACCCAACACAACCCCCCCCCACAGGTCGCCCACCGCGCCGAccag                                                           >  1:5230474/1‑148 (MQ=255)
                                                                                              ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCAACCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGCCctgctg                                                        >  1:7345738/1‑150 (MQ=255)
                                                                                               tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTgg                                                       >  2:5184121/1‑150 (MQ=255)
                                                                                               tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCCATCGCCGTGATCGGTGGCATCACCCTCCACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTgg                                                       >  2:7524819/1‑150 (MQ=255)
                                                                                               tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGAAAGCCCTGATCGGTGGCAACACCCACGACAAAGCCGCCCCACAGGTCGCCCACCGCGCCGACCAGCTgg                                                       >  1:6978057/1‑150 (MQ=255)
                                                                                               tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAACAGCCGATCGCCCTGATCGGTGGCATCACCCACGACAACGCCGCCCCACAGGTCGCCCACCCCGCCCCCCAGCTgg                                                       >  1:2988421/1‑150 (MQ=255)
                                                                                               tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCACGAGCAGGCCCGCGCCCAGGTGAAACTGCCCATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTgg                                                       >  1:1173820/1‑150 (MQ=255)
                                                                                                 aaTTCTGTCCCCCCGCCGGGCGCGCCCGCTGCCAACGCCGGCCTGCGCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCg                                                     <  2:7865112/150‑1 (MQ=255)
                                                                                                 aaTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCCGGTGAAAATGCCGATCGCCCAGATCGGTGGCATCACCCACGCCACCGCCCCCCCACAGGTCGCCCCCGGCGCCGACCAGCTGGCg                                                     >  1:2244010/1‑150 (MQ=255)
                                                                                                 aaTTCGGCCCCCACGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCg                                                     <  2:3048701/150‑1 (MQ=255)
                                                                                                  aTTCGGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCgg                                                    <  2:4726910/150‑1 (MQ=255)
                                                                                                          aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                            >  2:4918602/1‑150 (MQ=255)
                                                                                                          aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                            >  2:4756790/1‑150 (MQ=255)
                                                                                                           ccAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACg                                           >  2:1821842/1‑150 (MQ=255)
                                                                                                            caccccGGCCCCGCCCGCCCCCCCCGCCGGCCGGCCCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                          <  1:3852379/146‑1 (MQ=255)
                                                                                                             aaGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCGCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGc                                         >  1:3587226/1‑150 (MQ=255)
                                                                                                                 cGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGt                                     >  2:4251227/1‑150 (MQ=255)
                                                                                                                  gcgcccccccGCTCCCCACGACCGCCACCTCCACCCGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGtt                                    <  1:1495166/143‑1 (MQ=255)
                                                                                                                   ggCGCGCCCGCGGCCCCCGGCGGCCGGCGCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTc                                   <  2:8109417/150‑1 (MQ=255)
                                                                                                                      gcgcCCGCGGCCACCGGCGGCGGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGt                                <  2:762805/150‑1 (MQ=255)
                                                                                                                        gcCCGCTGCCCACGGCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGc                              <  2:2249584/150‑1 (MQ=255)
                                                                                                                          ccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCg                            >  2:6301555/1‑150 (MQ=255)
                                                                                                                          cagCTGCCAACGTCGGCAAGCTCGAGCAGGCCCGCGCCCACGTGAAACTGCCCATCGCCGTTATCGGTGTCATCAACATAGAAAACGCCGCCCCACAGGTCGCCCACGGCGCCCAACTGCTTGCGGTGAAACAAGGCTTGTTATGTGCCg                            >  2:3308562/3‑150 (MQ=255)
                                                                                                                           cGCTGCCCACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGa                           <  2:8792107/150‑1 (MQ=255)
                                                                                                                               gCCACCGGCGGCCGGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgc                       <  2:3329200/150‑1 (MQ=255)
                                                                                                                                ccAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcg                      >  1:7182738/1‑150 (MQ=255)
                                                                                                                                 cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc                     >  2:5126731/1‑150 (MQ=255)
                                                                                                                                     gTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAg                 <  2:7169268/150‑1 (MQ=255)
                                                                                                                                     gTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAg                 >  2:2310430/1‑150 (MQ=255)
                                                                                                                                     gTCGGCCGGCTCGTGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAg                 <  1:4360887/150‑1 (MQ=255)
                                                                                                                                      gcgaccggCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAgg                <  2:2398890/143‑1 (MQ=255)
                                                                                                                                      gcGGCCTGCACGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAgg                <  1:6137957/149‑1 (MQ=255)
                                                                                                                                        ggCCTGCCCGCCCCGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGaa              <  1:4454198/150‑1 (MQ=255)
                                                                                                                                            tGCCCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCa          <  1:2482780/150‑1 (MQ=255)
                                                                                                                                               tCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAccc       >  2:1859502/1‑150 (MQ=255)
                                                                                                                                                  aGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACccgcc    <  2:125048/150‑1 (MQ=255)
                                                                                                                                                   gCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACccgccg   >  1:8638401/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGCCGC  >  NC_002947/5451146‑5451443

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: