Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,787,629 C→G 6.1% intergenic (+76/‑49) ahpC → / → ahpF peroxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,787,6290CG 6.1% 208.6 / 2.5 64intergenic (+76/‑49)ahpC/ahpFperoxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base G (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

TGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGC  >  NC_002947/2787483‑2787769
                                                                                                                                                  |                                                                                                                                            
tgtgCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCt                                                                                                                                            >  7:337326/1‑149 (MQ=255)
        gCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTATGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACGTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCggg                                                                                                                                    <  8:179329/149‑1 (MQ=255)
        gCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTggg                                                                                                                                    <  2:190822/149‑1 (MQ=255)
        gCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTggg                                                                                                                                    >  2:86386/1‑149 (MQ=255)
        gCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTggg                                                                                                                                    <  3:316098/149‑1 (MQ=255)
                  aagaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGCGtt                                                                                                                                <  6:134858/143‑1 (MQ=255)
                  aagaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGCGtt                                                                                                                                >  5:134858/1‑143 (MQ=255)
                   agaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTAt                                                                                                                         >  7:312930/1‑149 (MQ=255)
                   agaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTAt                                                                                                                         >  5:53445/1‑149 (MQ=255)
                    gaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATc                                                                                                                        >  2:292991/1‑149 (MQ=255)
                    gaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATc                                                                                                                        <  5:373335/149‑1 (MQ=255)
                          gAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtt                                                                                                                  <  1:292991/149‑1 (MQ=255)
                            aGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtttt                                                                                                                <  5:228182/149‑1 (MQ=255)
                              ccACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTg                                                                                                               >  2:112628/1‑148 (MQ=255)
                                aCTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaa                                                                                                             >  7:125248/1‑148 (MQ=255)
                                    tGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                          <  8:312930/147‑1 (MQ=255)
                                    tGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                          <  8:125248/147‑1 (MQ=255)
                                    tGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAg                                                                                                         >  5:360668/1‑148 (MQ=255)
                                    tGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCAACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                          <  3:260543/147‑1 (MQ=255)
                                       cgcCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAAt                                                                                                     <  5:194869/149‑1 (MQ=255)
                                          cGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGc                                                                                                  >  6:56857/1‑149 (MQ=255)
                                              gCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCg                                                                                               >  7:91030/1‑148 (MQ=255)
                                                  gACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGt                                                                                          <  8:28355/149‑1 (MQ=255)
                                                  gACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGt                                                                                          >  3:339632/1‑149 (MQ=255)
                                                    ccTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTg                                                                                        <  5:56857/149‑1 (MQ=255)
                                                      tCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGa                                                                                      >  3:325772/1‑149 (MQ=255)
                                                             cAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACACCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACg                                                                                    >  8:8502/1‑144 (MQ=255)
                                                             cAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACACCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACg                                                                                    <  7:8502/144‑1 (MQ=255)
                                                              aaGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGAATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGc                                                                              >  3:316088/1‑149 (MQ=255)
                                                                     aaGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATc                                                                       >  8:36699/1‑149 (MQ=255)
                                                                      aGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATc                                                                       <  3:38048/148‑1 (MQ=255)
                                                                        gTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGc                                                                     >  2:47186/1‑148 (MQ=255)
                                                                            tGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTg                                                                <  4:325772/149‑1 (MQ=255)
                                                                                                    gccgccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCAcc                                        >  8:270825/1‑149 (MQ=255)
                                                                                                       gccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCa                                                                                 <  2:66089/105‑1 (MQ=255)
                                                                                                       gccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAg                                     >  4:125439/1‑149 (MQ=255)
                                                                                                          aaaGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGcc                                   >  3:63758/1‑148 (MQ=255)
                                                                                                               tGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                               >  1:68006/1‑147 (MQ=255)
                                                                                                                 aaGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGt                                            >  6:234310/1‑132 (MQ=255)
                                                                                                                 aaGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGt                                            <  5:234310/132‑1 (MQ=255)
                                                                                                                 aaGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCgaga                           <  6:360668/149‑1 (MQ=255)
                                                                                                                 aaGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCga                             <  4:63758/147‑1 (MQ=255)
                                                                                                                 aaGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCga                             <  4:339632/147‑1 (MQ=255)
                                                                                                                   gTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                               >  5:52764/1‑143 (MQ=255)
                                                                                                                      cctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGt                       <  7:36699/148‑1 (MQ=255)
                                                                                                                       ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                     <  6:53445/149‑1 (MQ=255)
                                                                                                                       ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                     <  1:112628/149‑1 (MQ=255)
                                                                                                                       ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGAATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                     <  4:316088/149‑1 (MQ=255)
                                                                                                                         accGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGc                    <  3:267893/148‑1 (MQ=255)
                                                                                                                         accGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGc                    >  4:267893/1‑148 (MQ=255)
                                                                                                                         accGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCt                   >  3:73611/1‑149 (MQ=255)
                                                                                                                          ccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCtt                  >  4:15757/1‑149 (MQ=255)
                                                                                                                               ccGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTc             <  1:47186/149‑1 (MQ=255)
                                                                                                                                 gTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcga           >  3:171195/1‑149 (MQ=255)
                                                                                                                                 gTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCAGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcga           >  8:321643/1‑149 (MQ=255)
                                                                                                                                  tCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcgac          >  4:280426/1‑149 (MQ=255)
                                                                                                                                  tCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcgac          <  4:73611/149‑1 (MQ=255)
                                                                                                                                     aaCGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCgacgac       <  7:321643/149‑1 (MQ=255)
                                                                                                                                      aCGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGacgacg      <  3:280426/149‑1 (MQ=255)
                                                                                                                                        gCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                     <  7:48957/132‑1 (MQ=255)
                                                                                                                                        gCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                     >  8:48957/1‑132 (MQ=255)
                                                                                                                                        gCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTccc               <  8:195599/138‑1 (MQ=255)
                                                                                                                                        gCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTccc               >  7:195599/1‑138 (MQ=255)
                                                                                                                                          ccaggCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAAATGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGgcgc  >  6:57393/4‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGC  >  NC_002947/2787483‑2787769

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: ‑