Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,448,265 A→C 16.2% E147A (GAG→GCG)  thiE → thiamine‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,448,2650AC16.2% 539.5 / 20.1 263E147A (GAG→GCG) thiEthiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (143/67);  new base C (15/33);  total (162/101)
Fisher's exact test for biased strand distribution p-value = 5.23e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

CCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGC  >  NC_002947/5448116‑5448413
                                                                                                                                                      |                                                                                                                                                    
    cAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAg                                                                                                                                                  <  1:4111063/150‑1 (MQ=255)
    cAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAg                                                                                                                                                  <  2:8715249/150‑1 (MQ=255)
     aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCCACGTCGGCCGGCGCGCGCAgg                                                                                                                                                 >  1:6540613/1‑150 (MQ=255)
     aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAgg                                                                                                                                                 >  2:3973405/1‑150 (MQ=255)
      ggCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGc                                                                                                                                                <  2:5685099/150‑1 (MQ=255)
       gCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGcc                                                                                                                                               >  2:3237045/1‑150 (MQ=255)
        cAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGATGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTATATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAGCGTCGGCCTGCTCGAGCAGGccg                                                                                                                                              <  1:7284529/150‑2 (MQ=255)
        cAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAGCGTCGGCCTGCTCGAGCAGGccg                                                                                                                                              <  1:1004485/150‑2 (MQ=255)
         aaTCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCACGcccc                                                                                                                                             >  1:7318884/1‑149 (MQ=255)
         aaTCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCGGCCAACGTCGCCCTGCGCGCGCCGGcccc                                                                                                                                             >  2:6340188/1‑149 (MQ=255)
         aaTCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTCCAATTCGGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGccga                                                                                                                                             <  2:6504222/150‑3 (MQ=255)
           tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCATGCTGCCAGCGAAGGCGCCAGCTACGTGGCATTTGGCCGCTTCTTAAATTCTGTCACCAAGCCGGGGGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCTGGccgacg                                                                                                                                           <  2:2361789/150‑5 (MQ=255)
           tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCAGGCCcgcg                                                                                                                                           >  2:3055552/1‑150 (MQ=255)
           tcatcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGATCAGGCCcgcg                                                                                                                                           >  2:5985512/1‑150 (MQ=255)
             atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGccgacgtt                                                                                                                                         <  2:760663/150‑7 (MQ=255)
             atcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGcc                                                                                                                                         >  1:7864643/1‑150 (MQ=255)
              tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCCCCAACCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGccc                                                                                                                                        <  2:3493088/150‑1 (MQ=255)
              tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCGGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGccc                                                                                                                                        <  2:7377939/150‑1 (MQ=255)
              tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGACTGCTCGAGCAGGCCCGCGccc                                                                                                                                        >  2:784868/1‑150 (MQ=255)
              tcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGACGGCGGGCGCGAGACCGCCCGCcccc                                                                                                                                        >  1:5655658/1‑150 (MQ=255)
                 gCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCgg                                                                                                                                     >  2:5953349/1‑150 (MQ=255)
                 gCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCACGAGCAGGCCCGCGCCCAgg                                                                                                                                     >  2:4740906/1‑150 (MQ=255)
                   tCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGCGCCGGCCCGCGcccccggg                                                                                                                                   >  1:6873079/1‑145 (MQ=255)
                   tCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGCGCCGGCCCCCCCCCacggg                                                                                                                                   >  1:6252077/1‑146 (MQ=255)
                      aCCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGaaa                                                                                                                                <  1:7489501/150‑1 (MQ=255)
                      aCCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTGGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGaaa                                                                                                                                <  2:2040246/150‑1 (MQ=255)
                          gCCCCGCCTGCCTCGAACGGGCCGCCCCGGCTGCCTGCTATGGCGCCAGATACGGGGCGTTTGGCCGCTTCTTAAATTCTGTCACCAAGCCGGGCGCCCCCGCTGCAACCGTCGGCTGGCGCGCGCAGGCCCGCGCCCAGGTGAAACTGc                                                                                                                            <  2:3812922/150‑1 (MQ=255)
                           ccACGCCAGCCTCGAACTGGCCGCCCTGGCTGCCAGCGACGGGGCCAGCTACGTGGCCTTGGGCCGCTTCTTCAATTCTGTCACCAACCCGGGCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGcc                                                                                                                           <  2:7705322/150‑1 (MQ=255)
                               gCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATc                                                                                                                       >  2:5538495/1‑150 (MQ=255)
                                 cAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGc                                                                                                                     >  2:7314837/1‑150 (MQ=255)
                                  aGCCTCGAACGGGCCGCCCAGGCTGCCAGCGACGGCGCCAGCTACGGGGCCTTTGGCCGCTACTCCATTTCTGCCACCAAGCGGGGCGCGCCCGCGGCCACCGCCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGcc                                                                                                                    <  2:1771726/150‑1 (MQ=255)
                                     cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCggg                                                                                                                 >  2:134004/1‑148 (MQ=255)
                                      tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGTGACACTGCCGATCGCCGTgg                                                                                                                >  2:7161935/1‑149 (MQ=255)
                                       ttcgcTGGCAGCCTGGGCGGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCACTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGAt                                                                                                               >  1:2917098/5‑150 (MQ=255)
                                             ggCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTgg                                                                                                         >  2:5594789/1‑150 (MQ=255)
                                                cgccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGCGCAGGCCCGCGCCCACGTGCAACAGCCGATCGCCGTGGTCGGTGGCAt                                                                                                      >  1:7030950/1‑150 (MQ=255)
                                                       gCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGGGTGCCGTGGCATCAcccgc                                                                                               >  2:1936292/1‑148 (MQ=255)
                                                         tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGTGAAACTGCCGCTCGCCGTGATCGGTGGCATCACCCACGa                                                                                             >  2:2613940/1‑150 (MQ=255)
                                                              gCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGTGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGCTGAAACAGCCGCTCGCCGGGATCGGTGGCCTCACCCACGACAAcc                                                                                        >  1:3330653/1‑149 (MQ=255)
                                                              gCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCCCCCAGGTGAAAATGACGATCGCCGAGATCGGTGGCATCACCCAAGACAAcc                                                                                        >  1:1134291/1‑149 (MQ=255)
                                                                gAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGGGTAACAGCCGATCGCCGGGATCGGTGGCATCACCCACGACAAagcc                                                                                      >  1:8936315/1‑150 (MQ=255)
                                                                       ccAGCTTGGTGGCGTTGGGCGGCTTCGTCATTTCGGCCCCCACGCCGGGCGCGCCCGCTCCCAACGGGGGCCGGCTCGCCCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCAc                                                                               <  1:6946962/150‑1 (MQ=255)
                                                                         cggtctttGGGCTTGGGCCTTTTCTTCTTTTCGGTACCCAGGCGGGGGGCGCCCGCGGCCACGGGCGGCGTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTg                                                                             <  1:2802806/143‑1 (MQ=255)
                                                                          gCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGGGATCGGGGGCATCACCCACGACAAAGCCCCCCCACAgg                                                                            >  1:1879288/1‑150 (MQ=255)
                                                                               gTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCCCCCCACAGGTCGcc                                                                       >  1:8332249/1‑150 (MQ=255)
                                                                               gTGGCCTTGGGCCGCTTCTTCATTTCCGCCACCCCCCCGGGCGCGCCCGCGCCCAACGCCGGCCGGCCCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGcc                                                                       <  1:133555/150‑1 (MQ=255)
                                                                                     tttGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCCCGTGAAACAGCCCATCGCCCGGATAGGTGGCATCACCCACCACAACGCCGCCCCACAGGTCGCCCAcccc                                                                 >  1:7190191/1‑147 (MQ=255)
                                                                                        ggCCGCTTCTTCAAGTCTGTCCCCACGCCGGGCGCGCCCGCTGCCACCGTCGGCGTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGcc                                                              <  2:8614166/150‑1 (MQ=255)
                                                                                          ccGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGa                                                            >  2:4630803/1‑150 (MQ=255)
                                                                                           cGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAc                                                           >  2:6548069/1‑150 (MQ=255)
                                                                                            gCTTCTTCAATTCTGTCACCCAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGGGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCCCGGGGCGCCCACGGCGCCGAcc                                                          >  2:7214363/1‑150 (MQ=255)
                                                                                              ttcttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGAACGCCGAGATCGGTGGCATAACCCACGACAAAGCCCCCCCACAGGTCGCCCACCGCGCCCAccgg                                                        >  1:1397616/1‑148 (MQ=255)
                                                                                                 ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACACCGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                     >  2:6234839/1‑150 (MQ=255)
                                                                                                 ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                     >  2:4874865/1‑150 (MQ=255)
                                                                                                 ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTAGAGCAGGCCCGCGCCCAGGTGAAACGGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACAGGTCGCCCACGGCGCCGACcagctg                                                     >  2:1000462/1‑150 (MQ=255)
                                                                                                 ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCACGAGCAGGCCCGCGCCCAGCGGTCAATGCCGATCGCCGAGGGCCGTGGCAGCACCCACGAAAACCCCCCCCCACGGCTCGCCCACCGCGCCGAccgcctg                                                     >  1:935829/1‑145 (MQ=255)
                                                                                                 ttcAACTCGGCCACCACGCCGGGCCCGCCCGCTGCCACCGTCGGCCTGCCCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                     <  1:8131140/150‑1 (MQ=255)
                                                                                                  tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGAACGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGCCGCCGACCAGCTgg                                                    >  2:7402165/1‑150 (MQ=255)
                                                                                                  tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCCTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACAGGTCGCCCACCGCGCCGACCAGCTgg                                                    >  1:2879932/1‑150 (MQ=255)
                                                                                                  tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCCTGATCGGTGGCATCACCCACGACAAAGCCGCCCCACAGGTCGCCCCCGGCGCCGCCCAGCTgg                                                    >  1:8625629/1‑150 (MQ=255)
                                                                                                  tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCCCCCACGAGAAACTGCCCATCGCCGGGATCGGTGGCAACACCCACGACAACCCCCCCCCACACGACGCCCACGCCGCCCACCAGCAgg                                                    >  1:8467114/1‑150 (MQ=255)
                                                                                                  tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCTCGAGCAGGCCCGCGCCCACGGGTAACAGCCCATCGCCGTGAGCGGAGGCGTAAACCTCGACAACACCCCCCCACAGGGCGCCCCCGCCGCCCCCCGGCTgg                                                    >  1:2323/1‑150 (MQ=255)
                                                                                                    aaTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCCTCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCg                                                  >  2:8478760/1‑150 (MQ=255)
                                                                                                     aTTCTGTCCCCACGCCGGGCGCGCCCGCCGCCAACGTCGGCCGGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCgg                                                 <  2:8255982/150‑1 (MQ=255)
                                                                                                     aTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCCCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCCCCCCCCAGGTCGCCCACGGCCCCGACCGGCTGGCgg                                                 >  1:2262056/1‑150 (MQ=255)
                                                                                                     aTTCTGCCCCCACGCCGGGCGCGCCCGCTGCCAACGCCGGCCGGCGCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCgg                                                 <  1:4540075/150‑1 (MQ=255)
                                                                                                        ctcaccccACGCGGGGCGCCCCCGCTGCCACGTCCGGCCGGATCCACCACGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATAGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGa                                              <  1:2455899/144‑1 (MQ=255)
                                                                                                        cTGTCCCCACGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGa                                              <  1:2102642/150‑1 (MQ=255)
                                                                                                          gTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCCGGCCCGCGCCCCGGGGCAACAGCCGCTCGCCGGGATCGGTGGCCTCACCCGCGACAACGCCGCCCCCCGGGTCGCCCACCGCGCCCACCGGCTGGCGGGGATc                                            >  1:5047713/1‑150 (MQ=255)
                                                                                                             cccACGCCGGGCCCCCCCGCGGCCCCCGCCGGCCGGCCCGCGCCGGCCCGCGCCCAGGTGAAACTGCTGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                         <  1:6315251/149‑1 (MQ=255)
                                                                                                             aCCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                         <  1:5975704/150‑1 (MQ=255)
                                                                                                             aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACACCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                         >  2:4974883/1‑150 (MQ=255)
                                                                                                               cAAGCCGGTCGCGCCCGCTGCCAACGTCGGCCTTCTCGATCAGGCCCGCGCCCATTTGAAACTGCCGATAGCCGTTATATGTGGCATCAACCTAGAAAACGCCGCCCCACTTGTCGCCCACTGCGCCGACCTGATGGCGGTGATACACgg                                       >  2:5453922/1‑150 (MQ=255)
                                                                                                                aaGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGCGCAGGCCCGCGCCCACGGGAAACTGCCGATAGCCCTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACGGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGc                                      >  1:3984203/1‑150 (MQ=255)
                                                                                                                aaGCCGGGCGCGCCCGCCGCCACCGGCGGCCTGCTCGAGCAGACCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGc                                      <  2:4362392/150‑1 (MQ=255)
                                                                                                                  gCCGGTCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCtt                                    <  2:3033035/150‑1 (MQ=255)
                                                                                                                  gCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCtt                                    <  2:2739750/150‑1 (MQ=255)
                                                                                                                   ccGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCCCAGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTg                                   >  1:8787267/1‑150 (MQ=255)
                                                                                                                     gggCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGtt                                 >  2:1528428/1‑150 (MQ=255)
                                                                                                                         gcgcCCGCTGCCCACGGCGGCCTGCTCGAGCTGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGt                             <  2:870253/150‑1 (MQ=255)
                                                                                                                         gcgcCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCCCTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGt                             >  1:3090925/1‑150 (MQ=255)
                                                                                                                           gcCCGCCGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGc                           <  1:8217364/150‑1 (MQ=255)
                                                                                                                           gcCCCCCGCCACCGGCGGCCGGCCCCCCCAGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGc                           <  1:5589386/150‑1 (MQ=255)
                                                                                                                                 tGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACGGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAg                     >  1:4408646/1‑150 (MQ=255)
                                                                                                                                   ccAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTg                                   >  2:4818113/1‑134 (MQ=255)
                                                                                                                                    cAACGTCGGCCTGCCCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc                  <  2:5742134/150‑1 (MQ=255)
                                                                                                                                    cAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc                  <  2:1153266/150‑1 (MQ=255)
                                                                                                                                     aaCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgcg                 >  1:3232924/1‑150 (MQ=255)
                                                                                                                                     aaCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCCACAACGCCGCCCCACTGGTCGCCCACGGCGCCCACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgcg                 >  1:8589630/1‑150 (MQ=255)
                                                                                                                                         gcgggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAg              >  1:3486981/4‑150 (MQ=255)
                                                                                                                                         gcgggCCGGCGCGCGCCGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAg              <  1:4788031/147‑1 (MQ=255)
                                                                                                                                            ggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGaa           <  2:3762967/150‑1 (MQ=255)
                                                                                                                                             gCCTGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAg          <  1:2998452/150‑1 (MQ=255)
                                                                                                                                              ccTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGt         >  1:5238540/1‑150 (MQ=255)
                                                                                                                                              ccTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGt         >  1:1954683/1‑150 (MQ=255)
                                                                                                                                               cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTc        >  1:4511193/1‑150 (MQ=255)
                                                                                                                                                  cTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAcc     <  2:8948728/150‑1 (MQ=255)
                                                                                                                                                  cTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAcc     <  2:6825378/150‑1 (MQ=255)
                                                                                                                                                   tCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAccc    >  2:1998491/1‑150 (MQ=255)
                                                                                                                                                     gAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACccgc  >  1:1701202/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGC  >  NC_002947/5448116‑5448413

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: