Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,304,627 G→A 7.7% A372A (GCC→GCT recB ← Chi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,304,6270GA 7.7% 428.6 / 17.5 143A372A (GCC→GCTrecBChi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base G (84/48);  new base A (2/9);  total (86/57)
Fisher's exact test for biased strand distribution p-value = 7.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CCGATCATGAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTG  >  NC_002947/5304478‑5304773
                                                                                                                                                     |                                                                                                                                                  
ccGATCATGAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCg                                                                                                                                                    <  1:4430608/150‑1 (MQ=255)
ccGATCATGAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCg                                                                                                                                                    <  2:2841941/150‑1 (MQ=255)
ccGATCATGAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCg                                                                                                                                                    <  2:2841954/150‑1 (MQ=255)
  gATCATGAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGc                                                                                                                                                  <  1:4612287/150‑1 (MQ=255)
        gAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCt                                                                                                                                            >  2:3914921/1‑150 (MQ=255)
         aaCAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTg                                                                                                                                           >  2:5140871/1‑150 (MQ=255)
          aCAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTgg                                                                                                                                          >  2:911146/1‑150 (MQ=255)
          aCAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTgg                                                                                                                                          >  1:4509333/1‑150 (MQ=255)
           cAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGc                                                                                                                                         <  2:493126/150‑1 (MQ=255)
              aCCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCAGCGAGGCCGGCGAg                                                                                                                                      <  2:5101173/150‑1 (MQ=255)
               ccGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCACGCGCTCGCCAGCGAGGCCGGCGAgc                                                                                                                                     <  2:1331815/150‑1 (MQ=255)
                    ttCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtg                                                                                                                                >  1:2611168/1‑150 (MQ=255)
                    ttCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtg                                                                                                                                <  1:2226524/150‑1 (MQ=255)
                      cTGCGCGATTTTCGCTGACCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtgct                                                                                                                              <  1:4645535/150‑1 (MQ=255)
                        gcgcGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTAACGGATCAGGCTGGCCAGGCGCTCGCCAGCGAGGCCGGCGAGCGCCtggcca                                                                                                                            <  2:4366514/150‑5 (MQ=255)
                        gcgcGATTTTCGCTGATCTGGTAGAGGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtgctgc                                                                                                                            <  1:3788248/150‑1 (MQ=255)
                         cgcgATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCa                                                                                                                           >  2:4091996/1‑150 (MQ=255)
                         cgcgATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGACCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCa                                                                                                                           >  2:149840/1‑150 (MQ=255)
                              ttttCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAgcgc                                                                                                                      <  1:446674/150‑1 (MQ=255)
                                  cGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGa                                                                                                                  >  2:3914266/1‑150 (MQ=255)
                                    cTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGc                                                                                                                <  1:1708580/150‑1 (MQ=255)
                                     tGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCa                                                                                                               >  1:5020659/1‑150 (MQ=255)
                                      gATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCAt                                                                                                              <  1:3914266/150‑1 (MQ=255)
                                          tGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCTCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcg                                                                                                          >  2:593007/1‑150 (MQ=255)
                                          tGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcg                                                                                                          >  2:3370191/1‑150 (MQ=255)
                                          tGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcg                                                                                                          >  1:2074033/1‑150 (MQ=255)
                                          tGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcg                                                                                                          >  2:1397686/1‑150 (MQ=255)
                                           ggTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcgt                                                                                                         >  2:1240436/1‑150 (MQ=255)
                                            gTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcgtc                                                                                                        >  1:824918/1‑150 (MQ=255)
                                             tAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGtcgtcg                                                                                                       <  2:4778499/150‑1 (MQ=255)
                                               gATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGaa                                                                                                     >  1:691372/1‑150 (MQ=255)
                                               gATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGaa                                                                                                     <  2:3911832/150‑1 (MQ=255)
                                                 tGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGc                                                                                                   >  2:2797216/1‑150 (MQ=255)
                                                 tGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGc                                                                                                   >  2:2797273/1‑150 (MQ=255)
                                                    gcTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCa                                                                                                >  1:2845063/1‑150 (MQ=255)
                                                    gcTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCa                                                                                                >  1:2845004/1‑150 (MQ=255)
                                                    gcTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCa                                                                                                >  2:2609389/1‑150 (MQ=255)
                                                    gcTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCa                                                                                                >  2:3821274/1‑150 (MQ=255)
                                                      tCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCAtt                                                                                              >  2:2425285/1‑150 (MQ=255)
                                                      tCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCAtt                                                                                              >  1:2391098/1‑150 (MQ=255)
                                                      tCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCAtt                                                                                              >  1:4803849/1‑150 (MQ=255)
                                                       cGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCAttt                                                                                             <  1:888218/150‑1 (MQ=255)
                                                        gAAGATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTc                                                                                            >  1:2201587/1‑150 (MQ=255)
                                                         aaaaTGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCt                                                                                           >  2:2045361/1‑150 (MQ=255)
                                                           aaTGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGc                                                                                         <  1:435394/150‑1 (MQ=255)
                                                                     tGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGgcgc                                                                               <  1:957438/150‑1 (MQ=255)
                                                                     tGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCACGCGCTCGCCAGCGAGGCCGGCGAGCGCCtg                                                                                                                                >  1:2467668/1‑101 (MQ=37)
                                                                     tGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCACGCGCTCGCCAGCGAGGCCGGCGAGCGCCtg                                                                                                                                <  2:2467668/101‑1 (MQ=37)
                                                                          cGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTATCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTtctc                                                                          >  1:2363003/1‑150 (MQ=255)
                                                                          cGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTtctc                                                                          >  1:2363059/1‑150 (MQ=255)
                                                                             gTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCAc                                                                       >  1:1325878/1‑150 (MQ=255)
                                                                                 gtgtCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCg                                                                   >  1:3993012/1‑150 (MQ=255)
                                                                                 gtgtCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCg                                                                   >  1:53814/1‑150 (MQ=255)
                                                                                    tCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCt                                                                                           <  1:3292530/123‑1 (MQ=255)
                                                                                    tCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGaaa                                                                >  1:3896253/1‑150 (MQ=255)
                                                                                      cTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAAcg                                                              >  2:2517996/1‑150 (MQ=255)
                                                                                             tcgtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGAt                                                       <  2:1008213/150‑1 (MQ=255)
                                                                                               gtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCgccac                                                     >  1:2648980/1‑146 (MQ=255)
                                                                                               gtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCgccac                                                     >  1:3947066/1‑146 (MQ=255)
                                                                                               gtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATAc                                                     >  2:1976972/1‑150 (MQ=255)
                                                                                               gtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATAc                                                     >  2:3342782/1‑150 (MQ=255)
                                                                                                  gATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCa                                                  >  1:3103444/1‑150 (MQ=255)
                                                                                                    tCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAAc                                                >  2:161660/1‑150 (MQ=255)
                                                                                                        cgcTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggc                                            >  1:3893677/1‑150 (MQ=255)
                                                                                                         gcTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcg                                           >  1:3155404/1‑150 (MQ=255)
                                                                                                         gcTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcg                                           <  2:5406066/150‑1 (MQ=255)
                                                                                                          cTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcgg                                          <  1:1111344/150‑1 (MQ=255)
                                                                                                          cTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcgg                                          <  1:1111127/150‑1 (MQ=255)
                                                                                                           tACCGGGAACTGCTCACGGATCAGGCTGGCCGGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcggc                                         >  2:3281756/1‑150 (MQ=255)
                                                                                                           tACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcggc                                         >  2:4372866/1‑150 (MQ=255)
                                                                                                           tACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcggc                                         <  1:5403301/150‑1 (MQ=255)
                                                                                                             ccGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGt                                       <  1:3370191/150‑1 (MQ=255)
                                                                                                              cGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTg                                      >  1:808280/1‑150 (MQ=255)
                                                                                                                   aCTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGa                                 <  1:1236861/150‑1 (MQ=255)
                                                                                                                     tGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGagc                               <  2:197689/150‑1 (MQ=255)
                                                                                                                          aCGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCgg                          >  2:356798/1‑150 (MQ=255)
                                                                                                                          aCGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCGGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCgg                          >  2:2589374/1‑150 (MQ=255)
                                                                                                                           cGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCggg                         <  2:4970396/150‑1 (MQ=255)
                                                                                                                            ggATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCgggg                        <  1:3762164/150‑1 (MQ=255)
                                                                                                                                 aGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATcc                   >  2:3197159/1‑150 (MQ=255)
                                                                                                                                   gCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCgg                 >  2:4858977/1‑150 (MQ=255)
                                                                                                                                   gCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCgg                 >  2:5220400/1‑150 (MQ=255)
                                                                                                                                   gCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCgg                 <  1:2797273/150‑1 (MQ=255)
                                                                                                                                   gCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCgg                 <  1:2797216/150‑1 (MQ=255)
                                                                                                                                    cTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCggg                >  2:2375400/1‑150 (MQ=255)
                                                                                                                                     tGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGc               >  1:4449109/1‑150 (MQ=255)
                                                                                                                                     tGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGc               >  2:3424247/1‑150 (MQ=255)
                                                                                                                                        ccAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGt            >  2:707070/1‑150 (MQ=255)
                                                                                                                                        ccAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGt            >  2:2283675/1‑150 (MQ=255)
                                                                                                                                         cAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTc           >  1:4679042/1‑150 (MQ=255)
                                                                                                                                          aGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTcc          >  1:580433/1‑150 (MQ=255)
                                                                                                                                            gcgcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCa                                                                        <  2:63020/86‑1 (MQ=255)
                                                                                                                                            gcgcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCa                                                                        >  1:63020/1‑86 (MQ=255)
                                                                                                                                             cgcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAACCGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAgg       <  2:4056312/150‑1 (MQ=255)
                                                                                                                                              gcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAggg      >  2:3506983/1‑150 (MQ=255)
                                                                                                                                              gcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAggg      >  1:4284430/1‑150 (MQ=255)
                                                                                                                                                tCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggc                                            >  2:4595219/1‑110 (MQ=255)
                                                                                                                                                tCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggc                                            <  1:4595219/110‑1 (MQ=255)
                                                                                                                                                 cGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCt   <  1:1397686/150‑1 (MQ=255)
                                                                                                                                                  gCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTg  >  2:4595302/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGATCATGAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTG  >  NC_002947/5304478‑5304773

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: