Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 541,447 G→C 15.6% L1274L (CTG→CTC 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,4470GC15.6% 577.2 / 87.1 211L1274L (CTG→CTCrpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (88/90);  new base C (21/12);  total (109/102)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 215 total aligned reads displayed.

TGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGT  >  NC_002947/541299‑541596
                                                                                                                                                    |                                                                                                                                                     
  gTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggt                                                                                                                                                    <  2:3766754/150‑1 (MQ=255)
      tGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTCACCCAGCGGCTGCTGggtaacc                                                                                                                                                <  1:4439297/150‑5 (MQ=255)
       gTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggtgacca                                                                                                                                               <  2:2277045/150‑5 (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)
            aTGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcgc                                                                                                                                          >  2:3461823/1‑150 (MQ=255)
                  gTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCg                                                                                                                                    <  2:2775619/150‑1 (MQ=255)
                     ccGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggtg                                                                                                                                 <  2:3098647/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)
                                 tCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCg                                                                                                                     <  2:845219/150‑1 (MQ=255)
                                  cGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                                                    >  2:3951022/1‑150 (MQ=255)
                                  cGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                                                    <  1:1714482/150‑1 (MQ=255)
                                  cGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgg                                                                                                                    >  1:4061547/1‑150 (MQ=255)
                                   gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                                                   >  2:2376870/1‑150 (MQ=255)
                                   gAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                                                   <  1:962444/150‑1 (MQ=255)
                                    aGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCgggg                                                                                                                  <  1:1925035/150‑1 (MQ=255)
                                          aggtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTGCGGCGAAATgg                                                                                                            >  2:246249/3‑150 (MQ=255)
                                           tgtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGa                                                                                                           <  1:3307265/150‑1 (MQ=255)
                                            gtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGag                                                                                                          <  1:2879363/150‑2 (MQ=255)
                                            gtgACCATGGGCTACATCTATATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGag                                                                                                          <  2:2839629/150‑2 (MQ=255)
                                            gtgACCATGGGCTACATCTAGATGTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGag                                                                                                          <  1:1256792/150‑2 (MQ=255)
                                              gACCTTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCTTTCCACTGGTTCCTACAGCCTGGTTACCCAGAAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGAGTTTCGGGGAGATGGAAgt                                                                                                        <  2:255767/150‑1 (MQ=255)
                                                ccGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAgtgt                                                                                                      <  2:4876375/150‑1 (MQ=255)
                                                 catgggctacatctatatgTTGAAGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGtctg                                                                                                     <  2:626990/143‑1 (MQ=255)
                                                 cGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGtgtg                                                                                                     >  2:3838992/1‑150 (MQ=255)
                                                  gTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTgg                                                                                                    <  2:4831357/150‑1 (MQ=255)
                                                       ttACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCt                                                                                               <  2:2125482/150‑1 (MQ=255)
                                                           atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                           <  2:2280683/150‑1 (MQ=255)
                                                           atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                           >  2:5318073/1‑150 (MQ=255)
                                                              tacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCCGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCa                                                                                        >  2:4886764/1‑150 (MQ=255)
                                                               acatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCAt                                                                                       >  1:3790667/1‑150 (MQ=255)
                                                                 atgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCggt                                                                                                                                  <  1:2854657/105‑1 (MQ=255)
                                                                 atgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTgg                                                                                             >  1:806014/1‑142 (MQ=255)
                                                                     tCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcg                                                                                 >  1:5434859/1‑150 (MQ=255)
                                                                     tCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcg                                                                                 >  1:1268887/1‑150 (MQ=255)
                                                                     tCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGCTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcg                                                                                 >  2:213065/1‑150 (MQ=255)
                                                                      cAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcgc                                                                                <  1:4989222/150‑1 (MQ=255)
                                                                      cAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcgc                                                                                >  2:1475905/1‑150 (MQ=255)
                                                                      cAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAATTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGgcgc                                                                                <  2:5406008/150‑1 (MQ=255)
                                                                       aaGCTCCACCACCTGGTCGACGACAATATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGtgcg                                                                               >  1:2956731/1‑150 (MQ=255)
                                                                       aaGCTCCACCACCTGGTCGACGACAATATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGtgcg                                                                               >  1:2956879/1‑150 (MQ=255)
                                                                         gctccaccaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTCCTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGc                                                                             >  1:937130/6‑150 (MQ=255)
                                                                         gCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGc                                                                             >  2:3911093/1‑150 (MQ=255)
                                                                         gCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGc                                                                             >  2:1163108/1‑150 (MQ=255)
                                                                          cTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCa                                                                            >  2:2426362/1‑150 (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:5128798/150‑1 (MQ=255)
                                                                                 acctggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTacacgc                                                                     >  1:2523015/4‑148 (MQ=255)
                                                                                   ttGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTg                                                                   <  2:1347004/150‑1 (MQ=255)
                                                                                   ctggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTg                                                                   >  2:3018278/2‑150 (MQ=255)
                                                                                     ggtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCa                                                                 >  1:520338/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)
                                                                                            gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAGGCATACGGCGCGGCATACACCCTGCAAGAAATg                                                          >  2:2022904/1‑150 (MQ=255)
                                                                                              cAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCt                                                        >  1:3827878/1‑150 (MQ=255)
                                                                                                aGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTca                                                      <  2:1871848/150‑1 (MQ=255)
                                                                                                 gATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTcac                                                     <  1:5316635/150‑1 (MQ=255)
                                                                                                  aTGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTcaca                                                    <  1:1014456/150‑1 (MQ=255)
                                                                                                     cACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTg                                                 <  2:15741/150‑1 (MQ=255)
                                                                                                      aCGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGa                                                <  2:520338/150‑1 (MQ=255)
                                                                                                        gcgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAg                                              <  2:4061547/150‑1 (MQ=255)
                                                                                                          gcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTc                                            >  1:692686/1‑150 (MQ=255)
                                                                                                          gcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTc                                            <  1:849163/150‑1 (MQ=255)
                                                                                                               ccACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacg                                       >  1:2018347/1‑150 (MQ=255)
                                                                                                                 aCTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCggg                                                                                                                   >  1:1758428/1‑72 (MQ=255)
                                                                                                                  cTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGacgacg                                    >  2:4910745/1‑150 (MQ=255)
                                                                                                                   tGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGt                                   >  1:3577088/1‑150 (MQ=255)
                                                                                                                   tGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGt                                   >  2:3279767/1‑150 (MQ=255)
                                                                                                                   tGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGt                                   >  2:4315175/1‑150 (MQ=255)
                                                                                                                   tGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGATGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGt                                   >  2:646807/1‑150 (MQ=255)
                                                                                                                     gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGa                                 >  1:616601/1‑150 (MQ=255)
                                                                                                                           aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGcc                           >  2:5032982/1‑150 (MQ=255)
                                                                                                                               ccTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTca                                                      <  2:2184288/119‑1 (MQ=255)
                                                                                                                                 tGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCa                     >  1:232776/1‑150 (MQ=255)
                                                                                                                                  taccaCCCAGCGGCTGCTGGGTAGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCaa                    >  1:4231095/5‑150 (MQ=255)
                                                                                                                                  cgtcaCCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCaa                    >  2:4149745/5‑150 (MQ=255)
                                                                                                                                   gTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTca                                                      <  1:301187/115‑1 (MQ=255)
                                                                                                                                     tACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGAt                 >  1:215750/1‑150 (MQ=255)
                                                                                                                                      aCCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAg                                              <  1:306171/120‑1 (MQ=255)
                                                                                                                                        ccAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTa              >  2:3548499/1‑150 (MQ=255)
                                                                                                                                        ccAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTa              >  1:4810025/1‑150 (MQ=255)
                                                                                                                                        ccAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTCGCAGGCCGTACCAAggtcta              >  2:426382/1‑145 (MQ=255)
                                                                                                                                         cagcagCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTAc             >  2:3612351/1‑150 (MQ=255)
                                                                                                                                           gcagcCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACaa           >  2:3162939/1‑150 (MQ=255)
                                                                                                                                            cagcCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAAt                                                           >  2:873836/1‑101 (MQ=255)
                                                                                                                                             agcCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGa         >  2:2997234/1‑150 (MQ=255)
                                                                                                                                             agcCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAGGTCTAcgaag         <  2:2932450/150‑5 (MQ=255)
                                                                                                                                              gcCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGaa        <  2:1509088/150‑1 (MQ=255)
                                                                                                                                              gcCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAGGTCTAcgaagc        <  1:4721466/150‑6 (MQ=255)
                                                                                                                                                   tGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCg   >  2:4388699/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGT  >  NC_002947/541299‑541596

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: