Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,304,595 T→C 8.1% E886G (GAG→GGG)  recC ← multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,304,5950TC 8.1% 173.7 / 5.6 62E886G (GAG→GGG) recCmultifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit
Reads supporting (aligned to +/- strand):  ref base T (24/33);  new base C (3/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGG  >  NC_002947/5304447‑5304737
                                                                                                                                                    |                                                                                                                                              
GAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGct                                                                                                                                                <  6:158512‑M1/149‑1 (MQ=255)
 AAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGctc                                                                                                                                               >  1:136132‑M1/1‑149 (MQ=255)
 AAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGctc                                                                                                                                               >  3:50764‑M1/1‑149 (MQ=255)
 AAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGctc                                                                                                                                               >  5:290146‑M1/1‑149 (MQ=255)
     cGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCg                                                                                                                                           <  8:204354/149‑1 (MQ=255)
       gCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTc                                                                                                                                            <  3:58018/146‑1 (MQ=255)
       gCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTc                                                                                                                                            >  4:58018/1‑146 (MQ=255)
            cgcgcccaggcccgacgcctgcaggcctgcggtctgctgccgctagccggtttcggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCggg                                                                                                                                     <  7:28620‑M2/90‑1 (MQ=255)
             ggCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGt                                                                                                                                   <  4:178998/149‑1 (MQ=255)
                agcagACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCgg                                                                                                                                >  7:184970/1‑149 (MQ=255)
                agcagACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCgg                                                                                                                                <  8:54319/149‑1 (MQ=255)
                agcagACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCgg                                                                                                                                <  4:234186/149‑1 (MQ=255)
                     aCCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCtt                                                                                                                           <  8:184970/149‑1 (MQ=255)
                      ccGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTc                                                                                                                          <  8:48006/149‑1 (MQ=255)
                       cGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCg                                                                                                                                       >  6:346438/1‑135 (MQ=255)
                       cGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCg                                                                                                                                       <  5:346438/135‑1 (MQ=255)
                        gCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTc                                                                                                                          >  5:210141/1‑147 (MQ=255)
                            gCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAgcgc                                                                                                                    >  3:62459/1‑149 (MQ=255)
                             ccTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGcc                                                                                                                   >  6:272425/1‑149 (MQ=255)
                              cTGCAGGCGTCAGGACTGGGCGCGCAAGACTTGTTCGACGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGcc                                                                                                                   >  8:54521/1‑148 (MQ=255)
                                  aGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCaaa                                                                                                              <  5:158512/149‑1 (MQ=255)
                                      gTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTac                                                                                                          <  8:276619/149‑1 (MQ=255)
                                          ggtctgctgccgctagccggtttcggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCt                                                                                                      <  6:23924‑M2/121‑1 (MQ=255)
                                                  cgcgcAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAgcgc                                                                                              <  1:12474/149‑1 (MQ=255)
                                                     gcAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctgg                                                                                           >  4:38806/1‑149 (MQ=255)
                                                       aaGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggct                                                                                         >  6:54799/1‑149 (MQ=255)
                                                        aGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCgctggct                                                                                         >  5:181749/1‑148 (MQ=255)
                                                            ttGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTc                                                                                                    <  7:295840/133‑1 (MQ=255)
                                                             tGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCgg                                                                                                                                      >  4:159203/1‑98 (MQ=255)
                                                             tGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCgg                                                                                                                                      <  3:159203/98‑1 (MQ=255)
                                                               ttcggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGt                                                                                 <  5:54799‑M2/142‑1 (MQ=255)
                                                                  ggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGGGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGc                                                                               <  6:181749‑M2/144‑1 (MQ=255)
                                                                  ggcgTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGc                                                                               <  1:125661‑M2/144‑1 (MQ=255)
                                                                         tCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCaa                                                                                                  >  1:71205/1‑122 (MQ=255)
                                                                         tCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCaa                                                                                                  <  2:71205/122‑1 (MQ=255)
                                                                           cGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGt                                                                     >  4:110577/1‑149 (MQ=255)
                                                                                    cAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaag                                                             >  5:70962/1‑148 (MQ=255)
                                                                                     aaGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCTCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaagaa                                                           <  7:54521/149‑1 (MQ=255)
                                                                                     aaGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCaagaa                                                           <  6:356723/149‑1 (MQ=255)
                                                                                          ggcgCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCgg                                                                       <  4:147993/132‑1 (MQ=255)
                                                                                          ggcgCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCgg                                                                       >  3:147993/1‑132 (MQ=255)
                                                                                             gcCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCa                                                   <  3:38806/149‑1 (MQ=255)
                                                                                                 aGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCAtt                                                     >  5:339629/1‑143 (MQ=255)
                                                                                                 aGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCAtt                                                     <  6:339629/143‑1 (MQ=255)
                                                                                                      gCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTg                                            >  3:248078/1‑147 (MQ=255)
                                                                                                        tttCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGa                                        >  2:341015/1‑149 (MQ=255)
                                                                                                           cACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACt                                      <  2:139488/148‑1 (MQ=255)
                                                                                                            aCTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCa                                    <  6:70962/149‑1 (MQ=255)
                                                                                                             cTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAg                                   <  3:128735/149‑1 (MQ=255)
                                                                                                                              cAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGAt                  <  5:272425/149‑1 (MQ=255)
                                                                                                                                gggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCa                            >  8:271316/1‑137 (MQ=255)
                                                                                                                                gggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCa                            <  7:271316/137‑1 (MQ=255)
                                                                                                                                 ggTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAg               >  3:167868/1‑149 (MQ=255)
                                                                                                                                      cccggACGAAGAGCCCTTCGTCCTGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCa                                                       >  7:287968/4‑104 (MQ=255)
                                                                                                                                      cccggACGAAGAGCCCTTCGTCCTGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCa                                                       <  8:287968/101‑1 (MQ=255)
                                                                                                                                       gAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCa                            <  3:232060/130‑1 (MQ=255)
                                                                                                                                       gAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGCGCTTCa                            >  4:232060/1‑130 (MQ=255)
                                                                                                                                       gAGGACGAAGGGCTCTTCGTCC‑GGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCACCTGGTCCAGACTCAGGCCTTCAGCGTCCGGATAAggcggc          >  6:90234/1‑147 (MQ=255)
                                                                                                                                        cggACGAAGAGCCCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATaa                <  6:69371/140‑1 (MQ=255)
                                                                                                                                        cggACGAAGAGCCCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATaa                >  5:69371/2‑141 (MQ=255)
                                                                                                                                        aGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCt                     <  3:294552/136‑1 (MQ=255)
                                                                                                                                        aGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCt                     >  4:294552/1‑136 (MQ=255)
                                                                                                                                          gACGAAGAGCCCTTCGTCCGGGGTCAGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGc      >  3:263219/1‑149 (MQ=255)
                                                                                                                                            cGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGccc    >  4:91974/1‑149 (MQ=255)
                                                                                                                                             gAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCg   <  6:27367/149‑1 (MQ=255)
                                                                                                                                              aaGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCgg  <  3:110577/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GAAACCGGCTAGCGGCAGCAGACCGCAGGCCTGCAGGCGTCGGGCCTGGGCGCGCAAGGCTTGTTCGGCGTTGTCCGGTTCGGCCAAGGCGGCGCCAAGCAGGCTTTCACTGGCGCCGTAGCGTTGCAGGGTATCGAGGACGAAGGGCTCTTCGTCCGGGGTCGGCGCTTCCAGCGCCTCAAAGTACACCTTCAAGCGCTGGCTGAAGAAGTGCCGTACCGGGTGGCGCAAGAAGTCATTCAGCTGGTGCAGACTCAGGGCTTCAGCGTCCTGATAAGGCGGCAAGCCCGG  >  NC_002947/5304447‑5304737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑