Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,304,635 T→C 7.6% S370G (AGC→GGC)  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,6350TC 7.6% 442.2 / 17.4 145S370G (AGC→GGC) recBChi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base T (84/50);  new base C (2/9);  total (86/59)
Fisher's exact test for biased strand distribution p-value = 7.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 147 total aligned reads displayed.

GAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTGCAGTTGCTGC  >  NC_002947/5304486‑5304783
                                                                                                                                                     |                                                                                                                                                    
gAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCt                                                                                                                                                      >  2:2588440/1‑150 (MQ=255)
 aaCAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTg                                                                                                                                                     >  2:5140871/1‑150 (MQ=255)
  aCAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCGTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTgg                                                                                                                                                    >  1:3065925/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)
         ggTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCAGCGAGGCCGGCGAgcgc                                                                                                                                             <  1:3915566/150‑1 (MQ=255)
            ttCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtg                                                                                                                                          <  1:2226524/150‑1 (MQ=255)
            ttCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtg                                                                                                                                          >  1:2611168/1‑150 (MQ=255)
              cTGCGCGATTTTCGCTGACCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtgct                                                                                                                                        <  1:4645535/150‑1 (MQ=255)
               tGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtgctg                                                                                                                                       >  1:4590797/1‑150 (MQ=255)
                gcgcGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTAACGGATCAGGCTGGCCAGGCGCTCGCCAGCGAGGCCGGCGAGCGCCtggcca                                                                                                                                      <  2:4366514/150‑5 (MQ=255)
                gcgcGATTTTCGCTGATCTGGTAGAGGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGtgctgc                                                                                                                                      <  1:3788248/150‑1 (MQ=255)
                 cgcgATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCa                                                                                                                                     >  1:5446933/1‑150 (MQ=255)
                 cgcgATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGACCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCa                                                                                                                                     >  2:149840/1‑150 (MQ=255)
                    gATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAgc                                                                                                                                  <  1:4376566/150‑1 (MQ=255)
                          cGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGa                                                                                                                            >  2:3914266/1‑150 (MQ=255)
                            cTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGc                                                                                                                          <  1:1708580/150‑1 (MQ=255)
                                  tGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCTCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcg                                                                                                                    >  2:593007/1‑150 (MQ=255)
                                  tGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcg                                                                                                                    >  1:2074033/1‑150 (MQ=255)
                                  tGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcg                                                                                                                    >  2:3370191/1‑150 (MQ=255)
                                  tGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcg                                                                                                                    >  2:1397686/1‑150 (MQ=255)
                                   ggTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATgtcgt                                                                                                                   >  2:1240436/1‑150 (MQ=255)
                                     tAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGtcgtcg                                                                                                                 <  2:4778499/150‑1 (MQ=255)
                                       gATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGaa                                                                                                               <  2:3911888/150‑1 (MQ=255)
                                         tGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGc                                                                                                             >  2:2797216/1‑150 (MQ=255)
                                            gcTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCa                                                                                                          >  2:3821274/1‑150 (MQ=255)
                                            gcTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCa                                                                                                          >  1:2845063/1‑150 (MQ=255)
                                            gcTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCa                                                                                                          >  2:2609389/1‑150 (MQ=255)
                                              tCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCAtt                                                                                                        >  1:2391098/1‑150 (MQ=255)
                                              tCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCAtt                                                                                                        >  1:4803849/1‑150 (MQ=255)
                                              tCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCAtt                                                                                                        >  2:2425285/1‑150 (MQ=255)
                                               cGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCAttt                                                                                                       <  1:888218/150‑1 (MQ=255)
                                                gAAGATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTc                                                                                                      >  1:2201587/1‑150 (MQ=255)
                                                        cGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGcc                                                                                              >  1:812812/1‑150 (MQ=255)
                                                             tGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGgcgc                                                                                         <  1:957438/150‑1 (MQ=255)
                                                             tGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCACGCGCTCGCCAGCGAGGCCGGCGAGCGCCtg                                                                                                                                          >  1:2467668/1‑101 (MQ=37)
                                                             tGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCACGCGCTCGCCAGCGAGGCCGGCGAGCGCCtg                                                                                                                                          <  2:2301677/101‑1 (MQ=37)
                                                               aaCCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCtt                                                                                       >  2:4163917/1‑150 (MQ=255)
                                                                  cGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTATCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTtctc                                                                                    >  1:2363003/1‑150 (MQ=255)
                                                                  cGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTtctc                                                                                    >  1:2363059/1‑150 (MQ=255)
                                                                     gTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCAc                                                                                 >  1:1325878/1‑150 (MQ=255)
                                                                         gtgtCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCg                                                                             >  1:53814/1‑150 (MQ=255)
                                                                            tCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGaaa                                                                          >  1:3896253/1‑150 (MQ=255)
                                                                              ntGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAAcg                                                                        >  1:4452990/2‑150 (MQ=255)
                                                                              cTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAAcg                                                                        >  2:2517996/1‑150 (MQ=255)
                                                                               tGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTAGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAAcgc                                                                       >  1:5378455/1‑150 (MQ=255)
                                                                                ggAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAAcgcg                                                                      <  1:1256170/150‑1 (MQ=255)
                                                                                 gAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACgcgc                                                                     >  2:4005212/1‑150 (MQ=255)
                                                                                     tcgtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGAt                                                                 <  2:1008213/150‑1 (MQ=255)
                                                                                       gtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCgccac                                                               >  1:3947066/1‑146 (MQ=255)
                                                                                       gtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCgccac                                                               >  1:2648980/1‑146 (MQ=255)
                                                                                       gtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATAc                                                               >  2:3342782/1‑150 (MQ=255)
                                                                                       gtcgATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATAc                                                               >  2:1976972/1‑150 (MQ=255)
                                                                                           aTCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCaa                                                           >  2:3762164/1‑150 (MQ=255)
                                                                                            tCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAAc                                                          <  1:4089453/150‑1 (MQ=255)
                                                                                            tCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAAc                                                          >  2:161660/1‑150 (MQ=255)
                                                                                                cgcTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggc                                                      >  1:3893677/1‑150 (MQ=255)
                                                                                                 gcTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcg                                                     >  1:3155491/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)
                                                                                                     ccGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGt                                                 <  1:3370191/150‑1 (MQ=255)
                                                                                                      cGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTg                                                >  1:808280/1‑150 (MQ=255)
                                                                                                           aCTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGa                                           <  1:911146/150‑1 (MQ=255)
                                                                                                           aCTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGa                                           <  1:1236861/150‑1 (MQ=255)
                                                                                                                tCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGagcagc                                      >  2:1987531/1‑150 (MQ=255)
                                                                                                                  aCGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCAACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCgg                                    >  2:2449568/1‑150 (MQ=255)
                                                                                                                   cGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCggg                                   <  2:4970396/150‑1 (MQ=255)
                                                                                                                    ggATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCgggg                                  <  1:3762164/150‑1 (MQ=255)
                                                                                                                       tCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCAt                               >  1:656592/1‑150 (MQ=255)
                                                                                                                         aGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATcc                             >  2:3197159/1‑150 (MQ=255)
                                                                                                                           gCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCgg                           <  1:2797216/150‑1 (MQ=255)
                                                                                                                           gCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCgg                           >  2:5220400/1‑150 (MQ=255)
                                                                                                                           gCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCgg                           <  1:2797273/150‑1 (MQ=255)
                                                                                                                           gCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCgg                           >  2:4858977/1‑150 (MQ=255)
                                                                                                                            cTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCggg                          <  2:2611168/150‑1 (MQ=255)
                                                                                                                             tGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGc                         >  2:3424247/1‑150 (MQ=255)
                                                                                                                             tGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGc                         >  1:4449109/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)
                                                                                                                                     cgcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAACCGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAgg                 <  2:4056312/150‑1 (MQ=255)
                                                                                                                                      gcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcg                                                     >  1:2804596/1‑113 (MQ=255)
                                                                                                                                      gcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggcg                                                     <  2:2804596/113‑1 (MQ=255)
                                                                                                                                      gcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAggg                >  1:4284430/1‑150 (MQ=255)
                                                                                                                                      gcTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAggg                >  2:3506983/1‑150 (MQ=255)
                                                                                                                                        tCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTggc                                                      >  2:4595219/1‑110 (MQ=255)
                                                                                                                                          gCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTg            >  2:4595302/1‑150 (MQ=255)
                                                                                                                                                 tCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTGCAGTtgc     <  1:4920880/150‑1 (MQ=255)
                                                                                                                                                 tCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACGTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTGCAGTtgc     >  1:648583/1‑150 (MQ=255)
                                                                                                                                                   gCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTGCAGTtgctg   <  2:3065925/150‑1 (MQ=255)
                                                                                                                                                    cTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTGCAGTtgctgc  <  1:3506983/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAACAGACCGGTTTCTGCGCGATTTTCGCTGATCTGGTAGATGCGCTCGAAAATGCCGTACTGAACCGGGTCGGTGTCCTGGAATTCGTCGATCAGCGCTACCGGGAACTGCTCACGGATCAGGCTGGCCAGGCGCTCGCCGGCCTCGCTGGCGAGCGCGTGCTGCAAGCGCACGAGCATGTCGTCGAAGCCCATTTCTGCACGCCGGCGCTTCTCCACTTCGAAACGCGCCGATACCCAACTGGCGGCGTGCTCGAGCAGCGGGGCATCCGGGCTGTCCAGGGCCTGCAGTTGCTGC  >  NC_002947/5304486‑5304783

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: