Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,448,259 T→G 12.8% 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,2590TG12.8% 362.0 / 6.7 192L145R (CTG→CGG) thiEthiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (89/70);  new base G (20/10);  total (111/81)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

CTCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCA  >  NC_002947/5448111‑5448408
                                                                                                                                                    |                                                                                                                                                     
 tCGGCCGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGc                                                                                                                                                     <  1:201744/150‑1 (MQ=255)
     ccGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGa                                                                                                                                                 <  2:1184586/150‑1 (MQ=255)
      cGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAg                                                                                                                                                <  2:3701759/150‑1 (MQ=255)
      cGTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAg                                                                                                                                                >  1:2321581/1‑150 (MQ=255)
       gTCAGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGc                                                                                                                                               >  1:396811/1‑150 (MQ=255)
          aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAgg                                                                                                                                            <  2:2670889/150‑1 (MQ=255)
          aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAgg                                                                                                                                            <  1:2502079/150‑1 (MQ=255)
          aGGCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAgg                                                                                                                                            >  2:51131/1‑150 (MQ=255)
            gCAATCATCGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGcc                                                                                                                                          >  1:4054872/1‑150 (MQ=255)
                 catcGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCgcgc                                                                                                                                     >  2:3493028/1‑150 (MQ=255)
                       cTCCACCTGCAACGCCTGCCTCGAACTGGCCGCCATGGCTGCCTGCTATGGCGCCTGCTAGGTGGCCTTTGGCCGCTTGTTCATTTGTGTCACAACGCCGGGCGCGCCCGCTGCCAAGGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGt                                                                                                                               <  2:560879/150‑1 (MQ=255)
                        tCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCAGCTAGAGCAGGCCCGCGCCCAGGTg                                                                                                                              >  2:690633/1‑150 (MQ=255)
                          cACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGaa                                                                                                                            >  2:2940340/1‑150 (MQ=255)
                           aCCTGCCACGCCAGCCTCGAACTGGCCGCCCTGGCTGCCAGCGACGGCGCCCGCTCCGTGGCCTTTGGCCGCTTCTTCAATTCGGGCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGaaa                                                                                                                           <  2:3634326/150‑1 (MQ=255)
                             cTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACt                                                                                                                         >  2:453353/1‑150 (MQ=255)
                               gCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCAGGTGAAACTGc                                                                                                                       >  1:1862069/1‑150 (MQ=255)
                                  aCGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCCGGCCCGCGCCCCCGTGAAACGGCCGa                                                                                                                    >  1:3621499/1‑150 (MQ=255)
                                  aCGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCCGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGa                                                                                                                    <  1:1049680/150‑1 (MQ=255)
                                       aGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCCCCCAGGTGCAACTGCCGATCGcc                                                                                                               >  1:1621820/1‑150 (MQ=255)
                                         ccTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGt                                                                                                             >  2:4406948/1‑150 (MQ=255)
                                         ccTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGt                                                                                                             >  2:2138973/1‑150 (MQ=255)
                                          cTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCCGGTGAAACGGCCGAGCGCCggg                                                                                                            >  2:2410779/1‑148 (MQ=255)
                                           tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCCCGCCCAGGTGAAACAGCCGATCGCCGTgc                                                                                                           >  1:4500908/1‑149 (MQ=255)
                                           tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCCGGCCCGCGCCCCCGTGCAACTGCCGATCGCCgggc                                                                                                           >  1:1566703/1‑147 (MQ=255)
                                             gAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATc                                                                                                         <  2:2321581/150‑1 (MQ=255)
                                               aCTGGCCGCCCTGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGCCACCAAGCCGGGCGCCCCCGCGGCCAACGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCgg                                                                                                       <  1:343305/150‑1 (MQ=255)
                                                cggcccgccCCGGCGGCCCGCGACGGCGCCCGCTACGGGGCCTTGGGCCGCTTCTTCCAGTCGGGCCCCACGCCGGGCGCGCCCGCGGCCCCCGGCGGCCGGCGCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGt                                                                                                      <  1:1764492/146‑1 (MQ=255)
                                                cTGGCAGCCTGGGCGGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGt                                                                                                      >  2:3298469/1‑150 (MQ=255)
                                                     cgccCTGGCTGCCCGCGACGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATGCTGCCACCAAGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCAt                                                                                                 <  2:2038789/150‑1 (MQ=255)
                                                      gccCCGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTGGGCCGCTTCTTCAATTCTGTCCCCACGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATc                                                                                                <  2:890227/150‑1 (MQ=255)
                                                       cccAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGGCTGCTCGAGCAGGCCCCCGCCCACGGGAAAATGCCGATCGCCCGGATCGGTGGCCTCa                                                                                               >  1:1723394/1‑150 (MQ=255)
                                                         cAGGCTGCCAGCGACGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCGGCCAACGCCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCAcc                                                                                             <  2:431041/150‑1 (MQ=255)
                                                          aGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCAccc                                                                                            >  2:1374419/1‑150 (MQ=255)
                                                            gCTGCCAGCTATGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTc                                                                                          <  2:1396175/150‑1 (MQ=255)
                                                              tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCCCCCACGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGa                                                                                        >  1:4396332/1‑150 (MQ=255)
                                                                ccAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCCGTGAAACTGCCGCTCGCCGGGATCGGTGGCAGCACCCGCGACa                                                                                      >  1:949919/1‑150 (MQ=255)
                                                                ccAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAAAAGCCGATAGCCGTGATCGGTGGCCTCACCCACGACa                                                                                      >  2:160687/1‑150 (MQ=255)
                                                                ccAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAAAAGCCGATCGCCGTGATCGGTGGCATCACCCACGACa                                                                                      >  1:4321751/1‑150 (MQ=255)
                                                                ccAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCCCGGGGAACTGCCGAGCGCCGTGTTCGGTGGCCTCACCCGCGACa                                                                                      >  1:402141/1‑150 (MQ=255)
                                                                 cAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCCGGCCCCCGCCCAGGGGAAACAGCCGCGCGCCGTGATCGGTGGCATCAACCACCACaa                                                                                     >  1:2104952/1‑150 (MQ=255)
                                                                  aGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAc                                                                                    >  2:4085420/1‑150 (MQ=255)
                                                                  aGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGGGAAACTGCCGAACGCCGGGATCGGTGGCCTCACCCTCGACAAc                                                                                    >  1:4121657/1‑150 (MQ=255)
                                                                  aGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCGTCTTCAATTCTGTCACCAATCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGACCCGGCCCCCGCCCACGGGAAACGTCCCATCCCCCTGAGCGGGGGCATCACCCAAGACAAc                                                                                    >  1:2801896/1‑150 (MQ=255)
                                                                   gCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCCGTGAAACGGCCGCGCGCCGTGATCGGTGGCATCACCCTCGACCAcg                                                                                   >  2:1293711/1‑150 (MQ=255)
                                                                   gCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAAcc                                                                                   >  2:3596213/1‑149 (MQ=255)
                                                                     gcaGGCGCCTGCTCCGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTGGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcgcc                                                                                 <  2:3201892/148‑1 (MQ=255)
                                                                      aaGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcgccg                                                                                >  2:3369695/1‑150 (MQ=255)
                                                                      aaGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCACGCCCCCCCCCCGGGGACACTGCCGGACGCCCTGGGCGGGGGCCTCACCCACCACAAcgccc                                                                                >  1:1124948/1‑149 (MQ=255)
                                                                        ggCGCCAGCTACGGGGCCTTTGGCCGATTCTTCACTTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACgccgcc                                                                              <  2:233210/150‑1 (MQ=255)
                                                                            ccAGCTACGTGGCCTTTGGCCGCTGCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGGGGCATCACCCTCGACACCGCCGCCCCAc                                                                          >  2:1651179/1‑150 (MQ=255)
                                                                             cctcTACGGGGCCTGTGGCCGCTTCTTCAATTCGGTCACCACGCGGGGCGCGCCCGCGGCCACCGGCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACt                                                                         <  1:3498806/147‑1 (MQ=255)
                                                                                cTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACTGGt                                                                      >  2:1854983/1‑150 (MQ=255)
                                                                                  aCGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCg                                                                    <  2:3744850/150‑1 (MQ=255)
                                                                                   ggTGGCCTTTGGCCGCTTCTACAATTCGCTCACCCCGCCGGGCGCGCCCGCCGCCACCGCCGGCCGGCGCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGc                                                                   <  1:916611/149‑1 (MQ=255)
                                                                                   cgcggCCTTGGGCCTCTTCTTCAATTCTGTCACCCCCCCGGCCGCGCCCGCGCCCAACGTCGGCCTGATCCCCCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGc                                                                   <  1:4351216/147‑1 (MQ=255)
                                                                                        cgtttGGCCGCTTCTTCATTGGTGTCACCCCGCCGGGCGCGCCCGCGGCCACCGTCGGCCGGCGCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACg                                                              <  1:3220939/148‑1 (MQ=255)
                                                                                        ccTTTGGCCGCTTCTTAAATTCTGTCACCAAGCGGGGCGCGCCCGCTGCCCAGGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACg                                                              <  2:4183551/150‑1 (MQ=255)
                                                                                         cTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCCCCCCACTGGTCGCCCACgg                                                             >  2:1926199/1‑150 (MQ=255)
                                                                                              gCCGCTTCATCAATCCGGTCACCCCGCCGGGCGCGCCCGCGGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCg                                                        <  2:3165726/150‑1 (MQ=255)
                                                                                                 gCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAcc                                                     >  2:257778/1‑150 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACAGGTAGCCCACGGCGCCGACctgctg                                                >  2:22491/1‑150 (MQ=255)
                                                                                                      ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAAGCCGCCCCACTGGTCGCCCACGGCGCCGACcagctg                                                >  1:3006546/1‑150 (MQ=255)
                                                                                                       tcAATTCTGTCCCCACGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTgg                                               <  2:2863752/150‑1 (MQ=255)
                                                                                                        cAATTCTGCCACCAAGCCGGGCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGc                                              <  1:1831101/150‑1 (MQ=255)
                                                                                                         caTTCTGTCACCCAGCCGGGCGCGCCCGCTGCCACCGGCGGCCTGCCCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGGCCTGCTGGCg                                             <  2:2822533/149‑1 (MQ=255)
                                                                                                         atctCTGCCCCCAAGCGGGGCGCGCCCGCGGCCCCCGTCGGCCCGCTCGCGCTGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCg                                             <  1:274722/147‑1 (MQ=255)
                                                                                                         aaTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCCCCCCCCAGGTCGCCCACGGCGCCGACCAGCTGGCg                                             >  1:1313583/1‑150 (MQ=255)
                                                                                                         aaTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCACGTGAAAAAGCCCATCGCCGTGATCGGTGGCATCAACCTAGACAACGCCCCCCCACAGGTCGCCCACGGCGCCCACCAGCTGGCg                                             >  1:106841/1‑150 (MQ=255)
                                                                                                          aTTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAAATGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCgg                                            >  1:1017139/1‑150 (MQ=255)
                                                                                                          aTTCGGGCACCCCGCCGGGCGCGCCCGCTGCCAACGGCGGCCTGCGCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCgg                                            <  1:1461571/150‑1 (MQ=255)
                                                                                                             cTGTCCCCACGCCGGGCGCGCCCGCGGCCAACGGCGGCCGGCCCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGa                                         <  2:3834323/150‑1 (MQ=255)
                                                                                                               gTCACCCCGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATc                                       <  2:2598177/150‑1 (MQ=255)
                                                                                                                tCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCCTCACCCACGACAACCCCCCCCCACAGGTCGCCCACGGCGCCGCCCTGCTGGCGGTGATca                                      >  1:4220712/1‑149 (MQ=255)
                                                                                                                  cccccgCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                    <  2:1219532/145‑1 (MQ=255)
                                                                                                                  aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                    >  2:1072946/1‑150 (MQ=255)
                                                                                                                   ccACGCCGGGCGCGCCCGCAGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACg                                   <  1:2583967/150‑1 (MQ=255)
                                                                                                                   ccAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACg                                   >  2:3540961/1‑150 (MQ=255)
                                                                                                                    cccgCCGGGCGCGCCCGCTCCCAACGCCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  <  1:695959/147‑1 (MQ=255)
                                                                                                                    cacgcggggCGCGCCCGCTGCCCCCGGCGGCCTGCCCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  <  2:205107/144‑1 (MQ=255)
                                                                                                                    cAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCGCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                  >  2:810748/1‑150 (MQ=255)
                                                                                                                     ccgCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGTGAAACGGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGc                                 <  1:3923844/148‑1 (MQ=255)
                                                                                                                     aaGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCCCCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGc                                 >  2:2149680/1‑150 (MQ=255)
                                                                                                                       gCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCtt                               <  2:2575588/150‑1 (MQ=255)
                                                                                                                          gggCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGtt                            <  2:3346629/150‑1 (MQ=255)
                                                                                                                            gcgcggCCGCTCCCAAGGGGGGCGTTCCCCCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCg                          <  1:1222535/150‑1 (MQ=255)
                                                                                                                                gcCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGc                      <  2:990818/150‑1 (MQ=255)
                                                                                                                                    gCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGAc                  >  2:314336/1‑150 (MQ=255)
                                                                                                                                      tGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAg                >  1:3746272/1‑150 (MQ=255)
                                                                                                                                        ccAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcg              >  1:2279013/1‑150 (MQ=255)
                                                                                                                                        ccAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcg              >  2:3015678/1‑150 (MQ=255)
                                                                                                                                            cgccGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGGGGCATCACCCTCGACAACGCCGCCCCACGGGTCGCCCCCGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCa          <  2:250122/147‑1 (MQ=255)
                                                                                                                                              gcGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAgg        <  1:3298469/149‑1 (MQ=255)
                                                                                                                                              gcGGCCTGCACGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAgg        <  1:2738493/149‑1 (MQ=255)
                                                                                                                                               gggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGa       >  2:1629305/2‑150 (MQ=255)
                                                                                                                                               gggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGa       >  2:2132839/2‑150 (MQ=255)
                                                                                                                                               gggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGa       >  1:2538473/2‑150 (MQ=255)
                                                                                                                                                ggCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGaa      <  2:296580/150‑1 (MQ=255)
                                                                                                                                                ggCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGaa      <  2:1707268/150‑1 (MQ=255)
                                                                                                                                                   cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTc   >  1:1804971/1‑150 (MQ=255)
                                                                                                                                                    tGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCa  <  2:3266056/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: ‑