Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,600 G→A 15.4% I446I (ATC→ATT rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,6000GA15.4% 427.7 / 61.3 175I446I (ATC→ATTrpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base G (64/84);  new base A (15/12);  total (79/96)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 177 total aligned reads displayed.

TCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTA  >  NC_002947/471454‑471749
                                                                                                                                                  |                                                                                                                                                     
tCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATc                                                                                                                                                    <  1:265350/150‑1 (MQ=255)
   tCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGaa                                                                                                                                                 >  2:1998942/1‑150 (MQ=255)
       gCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACggg                                                                                                                                             <  1:2719179/150‑1 (MQ=255)
            tCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATc                                                                                                                                        >  2:2307150/1‑150 (MQ=255)
               aTGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGcc                                                                                                                                     >  1:4839835/1‑150 (MQ=255)
                  ttttcagccAGTTCTTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGa                                                                                                                                  >  2:4252037/8‑150 (MQ=255)
                   gcTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGAc                                                                                                                                 <  2:4550410/150‑1 (MQ=255)
                   gcTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGAc                                                                                                                                 <  2:4550267/150‑1 (MQ=255)
                     tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACgg                                                                                                                               >  1:1693023/1‑150 (MQ=255)
                       agccAGTTCTTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAAt                                                                                                                             >  1:2236515/3‑150 (MQ=255)
                         ccAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCAGCGATCGAACGGGTGATCGCCTGACGGATcc                                                                                                                           >  2:178222/1‑150 (MQ=255)
                          cAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATcca                                                                                                                          >  1:2055879/1‑150 (MQ=255)
                           agttctTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACgaa                                                                                                                              >  1:1387548/4‑143 (MQ=255)
                           agttctTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATccac                                                                                                                         <  2:1959634/147‑1 (MQ=255)
                           aGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATccac                                                                                                                         <  2:962954/150‑1 (MQ=255)
                            gCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATccacc                                                                                                                        >  1:5103271/1‑150 (MQ=255)
                              tctTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAg                                                                                                                      <  2:2917047/150‑1 (MQ=255)
                              tctTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAg                                                                                                                      <  2:2917065/150‑1 (MQ=255)
                              tctTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAg                                                                                                                      >  1:787751/1‑150 (MQ=255)
                              tctTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAg                                                                                                                      >  1:4628271/1‑150 (MQ=255)
                                ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                                                                                                    >  2:2101830/1‑150 (MQ=255)
                                gtCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                                                                                                    <  2:5036328/149‑1 (MQ=255)
                                 tCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTg                                                                                                                   <  1:3519045/150‑1 (MQ=255)
                                 tCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGAACATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTg                                                                                                                   <  2:683124/150‑1 (MQ=255)
                                  ccGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTgg                                                                                                                  >  1:496560/1‑150 (MQ=255)
                                    ggAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCa                                                                                                                <  2:1439995/150‑1 (MQ=255)
                                       gttggTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAg                                                                                                             <  2:1171954/147‑1 (MQ=255)
                                        tCGGTTCGCGACCCAGTTCCTGCAGCCTCTGGCGGGAAATACGGTTGAGCTTGTGGCTCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAgg                                                                                                            <  2:3883250/150‑1 (MQ=255)
                                              cgcgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCgaga                                                                                                      <  1:3246117/150‑1 (MQ=255)
                                               gcgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGaa                                                                                                     <  2:1475334/150‑1 (MQ=255)
                                               gcgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGaa                                                                                                     <  2:1000879/150‑1 (MQ=255)
                                                  aCCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACtt                                                                                                  <  1:1649401/150‑1 (MQ=255)
                                                  aCCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACtt                                                                                                  <  1:3167685/150‑1 (MQ=255)
                                                  aCCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACtt                                                                                                  <  1:1649144/150‑1 (MQ=255)
                                                    ccATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGt                                                                                                >  1:2570065/1‑150 (MQ=255)
                                                    ccATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGGTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGt                                                                                                >  2:2707349/1‑150 (MQ=255)
                                                     gaTTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTa                                                                                               <  2:4001006/149‑1 (MQ=255)
                                                     cATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTa                                                                                               <  1:3504721/150‑1 (MQ=255)
                                                      aTTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAg                                                                                              <  2:4839835/150‑1 (MQ=255)
                                                      aTTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAg                                                                                              <  1:4129691/150‑1 (MQ=255)
                                                       tttCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGc                                                                                             >  1:66492/1‑150 (MQ=255)
                                                          ccTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCcgc                                                                                          <  1:1901053/150‑1 (MQ=255)
                                                          ccTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCcgc                                                                                          >  1:5079537/1‑150 (MQ=255)
                                                            tGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCg                                                                                                         <  1:3855467/133‑1 (MQ=255)
                                                            tGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGa                                                                                        >  2:36940/1‑150 (MQ=255)
                                                               agcaTCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCacgacga                                                                                     >  1:422142/1‑149 (MQ=255)
                                                                  aTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAt                                                                                  <  1:3604300/150‑1 (MQ=255)
                                                                  aTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAt                                                                                  <  1:1032565/150‑1 (MQ=255)
                                                                    cTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTc                                                                                <  1:4557610/150‑1 (MQ=255)
                                                                       gCGGAAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGaa                                                                             <  2:4335855/150‑1 (MQ=255)
                                                                        cGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAAc                                                                            <  1:458041/150‑1 (MQ=255)
                                                                        cGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAAc                                                                            <  2:1322153/150‑1 (MQ=255)
                                                                        cGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAAc                                                                            <  1:2762388/150‑1 (MQ=255)
                                                                        cGGGAAATACGGTTGAGCTTGTAGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAAc                                                                            <  1:2762599/150‑1 (MQ=255)
                                                                         gggAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACt                                                                           >  1:381974/1‑150 (MQ=255)
                                                                          ggAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACtt                                                                          >  1:3833862/1‑150 (MQ=255)
                                                                             aaTACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTc                                                                       <  2:810395/150‑1 (MQ=255)
                                                                             aaTACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTc                                                                       >  2:3368180/1‑150 (MQ=255)
                                                                                  gaacgatcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCg                                                                  >  2:4415920/8‑150 (MQ=255)
                                                                                  gaacgatcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCg                                                                  >  2:3140906/8‑150 (MQ=255)
                                                                                      gatcTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCctt                                                              >  1:4754418/4‑150 (MQ=255)
                                                                                       aGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtc                                                             <  1:3733230/150‑1 (MQ=255)
                                                                                       aGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtc                                                             <  2:3077301/150‑1 (MQ=255)
                                                                                        gCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtca                                                            >  2:2457623/1‑150 (MQ=255)
                                                                                          ttgttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcatc                                                          >  2:4027760/1‑150 (MQ=255)
                                                                                           tgttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTtcatca                                                         <  1:5032491/150‑1 (MQ=255)
                                                                                           acaagATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTtcatca                                                         >  1:1101340/5‑150 (MQ=255)
                                                                                            gttgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAg                                                        >  2:2467625/1‑150 (MQ=255)
                                                                                              tgATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAc                                                      <  2:4317414/150‑1 (MQ=255)
                                                                                               gATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGAcc                                                     >  2:4999756/1‑150 (MQ=255)
                                                                                                  cGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGAt                                                  <  1:2995531/150‑1 (MQ=255)
                                                                                                  cGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGAt                                                  <  1:2500638/150‑1 (MQ=255)
                                                                                                    tctcGATCATGTGCACTGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                >  1:4458993/1‑150 (MQ=255)
                                                                                                    tctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                >  2:3514725/1‑150 (MQ=255)
                                                                                                    tctcGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGt                                                >  2:1649865/1‑150 (MQ=255)
                                                                                                     ttcGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATAtt                                               <  1:3028926/149‑1 (MQ=255)
                                                                                                      tcGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTg                                              >  1:2332366/1‑150 (MQ=255)
                                                                                                        gatgatGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGcc                                            >  1:2905263/1‑150 (MQ=255)
                                                                                                         aTCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCt                                           <  1:1514430/150‑1 (MQ=255)
                                                                                                          tCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCtt                                          >  1:2160513/1‑150 (MQ=255)
                                                                                                          tCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCAACCGCTTTCATCAGACCGATGTTGCCtt                                          >  1:1233412/1‑150 (MQ=255)
                                                                                                             tgtgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCt                                       <  1:489692/150‑1 (MQ=255)
                                                                                                             tgtgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCt                                       >  2:2595459/1‑150 (MQ=255)
                                                                                                             tgtgCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCt                                       <  2:3847685/150‑1 (MQ=255)
                                                                                                               tgCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTgg                                     >  1:2150717/1‑150 (MQ=255)
                                                                                                                gCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGa                                    >  1:143161/1‑150 (MQ=255)
                                                                                                                  aCCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGACTTCCTGGATc                                  >  1:1343175/1‑150 (MQ=255)
                                                                                                                    cGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTgg                                     <  1:3694168/145‑1 (MQ=255)
                                                                                                                    cGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAg                                <  1:4013535/150‑1 (MQ=255)
                                                                                                                      gAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGt                              <  2:1059115/150‑1 (MQ=255)
                                                                                                                         tGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGa                           <  2:2177654/150‑1 (MQ=255)
                                                                                                                           cGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAgg                         >  1:3349699/1‑150 (MQ=255)
                                                                                                                            ggATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGa                        >  2:4129842/1‑150 (MQ=255)
                                                                                                                                          tGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACgg          <  2:3714058/150‑1 (MQ=255)
                                                                                                                                           ggTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGt         <  1:4336995/150‑1 (MQ=255)
                                                                                                                                             tCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTg       <  2:4501746/150‑1 (MQ=255)
                                                                                                                                             tCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTg       <  2:4950682/150‑1 (MQ=255)
                                                                                                                                               ggCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGACGATATTCGAACTTATCGACCGCCTTCATCAGACCGATATTGCCTTCCTGAATAAGATCGAGGAACTGAAGGCCACGGTTGgt     >  1:4566644/1‑150 (MQ=255)
                                                                                                                                                gCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgtg    <  1:2422943/150‑1 (MQ=255)
                                                                                                                                                gCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGgtg    >  1:3170611/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCTTTCATCAGACCGATGTTGCCTTCCTGGATCAGGTCGAGGAACTGCAGGCCACGGTTGGTGTA  >  NC_002947/471454‑471749

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: