Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 541,426 C→G 12.9% S1267R (AGC→AGG) ‡ rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947541,4260CG12.9% 588.2 / 66.8 201S1267R (AGC→AGG) ‡rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (77/98);  new base G (17/9);  total (94/107)
Fisher's exact test for biased strand distribution p-value = 5.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 205 total aligned reads displayed.

CTGCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCC  >  NC_002947/541280‑541571
                                                                                                                                                  |                                                                                                                                                 
cTGCCAGAGAGCGGCCAGATGGTGCTGTTCGCTGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCCCTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTg                                                                                                                                                <  2:3316086/150‑1 (MQ=255)
 tGCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTgg                                                                                                                                               >  1:40633/1‑150 (MQ=255)
  gCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGt                                                                                                                                              >  1:3139054/1‑150 (MQ=255)
   ccAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGtt                                                                                                                                             >  1:15741/1‑150 (MQ=255)
          gCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCcagc                                                                                                                                      <  2:1793780/150‑1 (MQ=255)
          gCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCcagc                                                                                                                                      <  1:1710358/150‑1 (MQ=255)
           cGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCcagca                                                                                                                                     <  1:4369674/150‑1 (MQ=255)
           cGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCcagca                                                                                                                                     <  2:2566336/150‑1 (MQ=255)
                     gTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGGTCCTACAGCCTGGTTACCCAGCAGCCGCTGggt                                                                                                                           <  2:2800911/150‑1 (MQ=255)
                         tGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAATCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTa                                                                                                                       <  2:3891192/150‑1 (MQ=255)
                         tGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTCACCCAGCGGCTGCTGggtaacc                                                                                                                       <  1:4439297/150‑5 (MQ=255)
                          gTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCg                                                                                                                                 <  2:3710823/139‑1 (MQ=255)
                          gTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggtgacca                                                                                                                      <  2:2277047/150‑5 (MQ=255)
                           ttCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAg                                                                                                                     >  1:563311/1‑150 (MQ=255)
                           ttCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAg                                                                                                                     <  1:944281/150‑1 (MQ=255)
                           ttCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACGTGGTGGCCGACCAGATGCACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCCGCAGCCGCGGGGTGGTagg                                                                                                                     >  2:3223418/1‑148 (MQ=255)
                             cGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTggg                                                                                                                            >  1:2457740/1‑141 (MQ=255)
                             cGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgc                                                                                                                   <  2:4452652/150‑1 (MQ=255)
                               aTGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcgc                                                                                                                 >  1:1509088/1‑150 (MQ=255)
                                tGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCa                                                                                                                <  1:3380726/150‑1 (MQ=255)
                                        ccGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtg                                                                                                        <  2:3098647/150‑1 (MQ=255)
                                        ccGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtg                                                                                                        <  1:3223418/150‑1 (MQ=255)
                                         cGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtgg                                                                                                       >  2:1742820/1‑150 (MQ=255)
                                          ggCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtggt                                                                                                      >  1:4515045/1‑150 (MQ=255)
                                          ggCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACACGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtggt                                                                                                      <  2:924524/150‑1 (MQ=255)
                                           gCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGc                                                                                                                                                   <  1:5372359/104‑1 (MQ=255)
                                             aacaaGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAg                                                                                                   >  2:3611924/1‑150 (MQ=255)
                                                aaGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGt                                                                                                >  2:4946355/1‑150 (MQ=255)
                                                      gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                          <  2:4384530/150‑1 (MQ=255)
                                                      gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                          <  1:962444/150‑1 (MQ=255)
                                                      gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                          <  1:4797084/150‑1 (MQ=255)
                                                      gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                          >  2:2376870/1‑150 (MQ=255)
                                                      gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                          <  2:2166608/150‑1 (MQ=255)
                                                       aGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgggg                                                                                         <  1:2856950/150‑1 (MQ=255)
                                                       aGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgggg                                                                                         <  1:1925035/150‑1 (MQ=255)
                                                             aggtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATgg                                                                                   >  1:4181326/3‑150 (MQ=255)
                                                             aggtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTGCGGCGAAATgg                                                                                   >  2:246249/3‑150 (MQ=255)
                                                               gtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGaa                                                                                 <  2:5366346/150‑1 (MQ=255)
                                                               gtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGag                                                                                 <  1:2879363/150‑2 (MQ=255)
                                                               gtgACCATGGGCTACATCTAGATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGag                                                                                 <  1:1256792/150‑2 (MQ=255)
                                                                   ccGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAgtgt                                                                             <  2:4876794/150‑1 (MQ=255)
                                                                    catgggctacatctatatgTTGAAGCTACACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCTCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGtctg                                                                            <  2:626986/143‑1 (MQ=255)
                                                                    cGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGtgtg                                                                            <  1:816680/150‑1 (MQ=255)
                                                                     gTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTgg                                                                           <  2:4831357/150‑1 (MQ=255)
                                                                          ttACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCt                                                                      >  1:5043219/1‑150 (MQ=255)
                                                                          ttACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCt                                                                      <  2:2125482/150‑1 (MQ=255)
                                                                          ttACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCt                                                                      <  1:784211/150‑1 (MQ=255)
                                                                          ctacatctatatgTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGctct                                                                      >  2:51238/2‑150 (MQ=255)
                                                                             catgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTGACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGAGGGAAGTGTGGGCGCTGGa                                                                   <  1:1099095/150‑1 (MQ=255)
                                                                              atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                  >  2:5318073/1‑150 (MQ=255)
                                                                              atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGGTCCGCTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                  <  2:2243847/150‑1 (MQ=255)
                                                                                gtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggt                                                                                                                           >  2:377200/1‑91 (MQ=255)
                                                                                gtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggt                                                                                                                           <  1:377200/91‑1 (MQ=255)
                                                                                 tacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCa                                                               >  2:4886773/1‑150 (MQ=255)
                                                                                 tacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCa                                                               <  1:2771475/150‑1 (MQ=255)
                                                                                    atgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggt                                                                                                         >  2:2854657/1‑105 (MQ=255)
                                                                                    atgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTgg                                                                    <  2:806014/142‑1 (MQ=255)
                                                                                    atgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATAc                                                            <  2:1949914/150‑1 (MQ=255)
                                                                                      gCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACgg                                                          <  1:1781359/150‑1 (MQ=255)
                                                                                        tCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCcagca                                                                                                                                     <  2:4167400/73‑1 (MQ=255)
                                                                                        tCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcg                                                        >  1:1268887/1‑150 (MQ=255)
                                                                                        tCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGCTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcg                                                        >  2:213065/1‑150 (MQ=255)
                                                                                         cAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcgc                                                       >  2:1475905/1‑150 (MQ=255)
                                                                                         cAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcgc                                                       >  2:4306030/1‑150 (MQ=255)
                                                                                          aaGCTCCACCACCTGGTCGACGACAATATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGtgcg                                                      >  1:2956731/1‑150 (MQ=255)
                                                                                          aaGCTCCACCACCTGGTCGACGACAATATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGtgcg                                                      >  1:2956879/1‑150 (MQ=255)
                                                                                            gctccaccaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTCCTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGc                                                    >  1:937130/6‑150 (MQ=255)
                                                                                            gctccaccaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGc                                                    <  1:3826552/145‑1 (MQ=255)
                                                                                             cTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCa                                                   >  2:2426362/1‑150 (MQ=255)
                                                                                              tccaccaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGctt                                                  >  1:643617/4‑148 (MQ=255)
                                                                                                aaCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATac                                                >  1:4099268/1‑150 (MQ=255)
                                                                                                 accaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTaca                                               >  1:3699655/1‑150 (MQ=255)
                                                                                                 aCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATaca                                               <  1:1475905/150‑1 (MQ=255)
                                                                                                  ccACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATacac                                              >  2:2652752/1‑150 (MQ=255)
                                                                                                  ccACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATacac                                              >  2:4661266/1‑150 (MQ=255)
                                                                                                   cACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACAcc                                             <  1:4913834/150‑1 (MQ=255)
                                                                                                    acctggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTacacgc                                            >  1:2523015/4‑148 (MQ=255)
                                                                                                      ttGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTg                                          <  2:3528817/150‑1 (MQ=255)
                                                                                                      ttGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTg                                          >  1:5179444/1‑150 (MQ=255)
                                                                                                      ctggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTg                                          >  2:3018278/2‑150 (MQ=255)
                                                                                                       tggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGc                                         <  2:643617/150‑1 (MQ=255)
                                                                                                       tggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGc                                         >  1:5223808/1‑150 (MQ=255)
                                                                                                        ggtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCa                                        >  2:849163/1‑150 (MQ=255)
                                                                                                         gtggACGACAAGATGCACGCGCGGTCCACTGGGTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                       <  2:3358999/150‑1 (MQ=255)
                                                                                                           ggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGa                                     <  1:958044/150‑1 (MQ=255)
                                                                                                            gacgacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGaa                                    <  2:4509018/150‑1 (MQ=255)
                                                                                                                     aTGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTcaca                           <  1:1014456/150‑1 (MQ=255)
                                                                                                                       gCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGt                         <  1:1742820/150‑1 (MQ=255)
                                                                                                                        cACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTg                        <  2:15741/150‑1 (MQ=255)
                                                                                                                         aCGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGa                       <  2:520338/150‑1 (MQ=255)
                                                                                                                         aCGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGa                       >  1:4021865/1‑150 (MQ=255)
                                                                                                                             gcgGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTc                   <  2:3927640/150‑1 (MQ=255)
                                                                                                                                  ccACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacg              >  1:2018347/1‑150 (MQ=255)
                                                                                                                                      tGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGt          >  1:3577088/1‑150 (MQ=255)
                                                                                                                                      tGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGt          >  2:3279767/1‑150 (MQ=255)
                                                                                                                                         ttCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGaa       >  2:4053517/1‑150 (MQ=255)
                                                                                                                                           ccTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACg     <  2:1165273/150‑1 (MQ=255)
                                                                                                                                             tACAGCCTGGCTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGc   >  2:5399575/1‑150 (MQ=255)
                                                                                                                                              aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGcc  <  1:3461823/150‑1 (MQ=255)
                                                                                                                                                  ccTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTca                             >  1:2184288/1‑119 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTGCCAGAGAGCGGCCAGATGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCC  >  NC_002947/541280‑541571

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: