Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,787,623 G→A 13.0% intergenic (+70/‑55) 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,6230GA13.0% 168.9 / 14.4 69intergenic (+70/‑55)ahpC/ahpFperoxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base A (4/5);  total (30/39)
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.10e-01

GGCGAAGTGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCG  >  NC_002947/2787476‑2787768
                                                                                                                                                   |                                                                                                                                                 
ggCGAAGTGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCgg                                                                                                                                                  <  7:148923/149‑1 (MQ=255)
   gAAGTGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGc                                                                                                                                                >  3:284489/1‑148 (MQ=255)
     aGTGTGCCCGGCCAAGTGGAAAGAAGTCGAATCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGAc                                                                                                                                             <  4:212442/149‑1 (MQ=255)
     aGTGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGAc                                                                                                                                             <  2:92808/149‑1 (MQ=255)
             cGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTg                                                                                                                                     >  8:374013/1‑149 (MQ=255)
              ggCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                    <  4:192453/149‑1 (MQ=255)
                    gTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttt                                                                                                                              <  8:255029/149‑1 (MQ=255)
                    gTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGCGttt                                                                                                                              <  3:165325/149‑1 (MQ=255)
                       gAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGc                                                                                                                                  >  7:342645/1‑142 (MQ=255)
                       gAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGc                                                                                                                                  <  8:342645/142‑1 (MQ=255)
                        aaaGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGCGtttgacg                                                                                                                          <  7:47155/149‑5 (MQ=255)
                         aagaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGt                                                                                                                          <  4:298040/148‑1 (MQ=255)
                         aagaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTa                                                                                                                         >  5:49083/1‑149 (MQ=255)
                          agaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTAt                                                                                                                        >  1:1715/1‑149 (MQ=255)
                           gaagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATc                                                                                                                       <  4:345515/149‑1 (MQ=255)
                            aagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTAt                                                                                                                        <  4:362166/147‑1 (MQ=255)
                            aagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTAt                                                                                                                        >  3:356333/1‑147 (MQ=255)
                            aagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTAt                                                                                                                        >  3:362166/1‑147 (MQ=255)
                            aagGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTAt                                                                                                                        <  4:356333/147‑1 (MQ=255)
                             agGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                    >  8:33126/1‑134 (MQ=255)
                             agGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                    <  7:33126/134‑1 (MQ=255)
                                 gAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                            <  7:84136/138‑1 (MQ=255)
                                 gAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                            >  8:84136/1‑138 (MQ=255)
                                 gAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGtt                                                                                                                 <  5:132078/149‑1 (MQ=255)
                                     ccACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGa                                                                                                             <  1:91316/149‑1 (MQ=255)
                                        ctctGGCGCCGTCGCTGGACCTCGTAGGCAAGATCAAACGTCATGCCCATAACGCGGGCAGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaa                                                                                                          <  5:22378/149‑1 (MQ=255)
                                           tGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTATGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                         <  8:26717/147‑1 (MQ=255)
                                           tGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                         <  2:1715/147‑1 (MQ=255)
                                             gcgcCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGaa                                                                                                     <  8:74469/149‑1 (MQ=255)
                                                      cTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTa                                                                                            <  6:89496/149‑1 (MQ=255)
                                                      cTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTa                                                                                            >  6:49038/1‑149 (MQ=255)
                                                               gTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                    <  2:82317/100‑1 (MQ=255)
                                                               gTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTgg                                                                                                                                    >  1:82317/1‑100 (MQ=255)
                                                               gTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                         <  8:238899/127‑1 (MQ=255)
                                                               gTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                         >  7:238899/1‑127 (MQ=255)
                                                               gTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGAc                                                                                    <  4:284489/148‑1 (MQ=255)
                                                                     aaGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGCATGTTGGACGCCACGc                                                                             <  7:374013/149‑1 (MQ=255)
                                                                            aaGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATcc                                                                                                                      >  2:61758/1‑101 (MQ=255)
                                                                            aaGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATcc                                                                                                                      <  1:61758/101‑1 (MQ=255)
                                                                             aGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGTCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTACATc                                                                      <  5:285364/148‑1 (MQ=255)
                                                                               gTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGc                                                                    >  7:354612/1‑148 (MQ=255)
                                                                                tCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCaa                                                                  <  6:49083/149‑1 (MQ=255)
                                                                                      ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGaa                                                             >  3:57988/1‑148 (MQ=255)
                                                                                      ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACCCCCGTGCGACCTAGCCTGTGCGTTTTTGTATAAGAGTTTTGAAAAAGTAAACGCCAGTATGTTGGACGCCACCCTTAAATAGCAACTGaa                                                             >  6:58703/1‑148 (MQ=255)
                                                                                       cATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGaaaa                                                           <  1:237266/149‑1 (MQ=255)
                                                                                           aCGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTa                                                       >  1:305318/1‑149 (MQ=255)
                                                                                              cgggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACtt                                                    <  5:49038/149‑1 (MQ=255)
                                                                                                              gccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAg                                    >  1:20721/1‑149 (MQ=255)
                                                                                                               ccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGc                                   >  1:18199/1‑149 (MQ=255)
                                                                                                                cAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGccc                                                                                               >  5:344050/1‑88 (MQ=255)
                                                                                                                cAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGccc                                                                                               <  6:344050/88‑1 (MQ=255)
                                                                                                                 aaaGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGcc                                  >  5:311893/1‑148 (MQ=255)
                                                                                                                  aaGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGa                                <  2:21369/149‑1 (MQ=255)
                                                                                                                    gCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAATTAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                              <  5:58703/149‑1 (MQ=255)
                                                                                                                             cctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGtt                     <  2:18199/149‑1 (MQ=255)
                                                                                                                              ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                    <  2:305318/149‑1 (MQ=255)
                                                                                                                                accGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCt                  >  2:232716/1‑149 (MQ=255)
                                                                                                                                 ccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAAc                                                          <  3:232884/108‑1 (MQ=255)
                                                                                                                                 ccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAAc                                                          >  4:232884/1‑108 (MQ=255)
                                                                                                                                   gCCCCGTCAAACGCCCGGGCGACCTCGACTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGATAACCGACTTGGAGCGGGTCACCAAGCAGATCGAGATCGTTGCTTcc               >  2:41383/1‑149 (MQ=255)
                                                                                                                                       acaaaaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTc            >  2:369611/5‑148 (MQ=255)
                                                                                                                                         aaaaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCaa                                                                  >  5:136469/3‑92 (MQ=37)
                                                                                                                                         aaaaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCaa                                                                  <  6:136469/90‑1 (MQ=37)
                                                                                                                                           aaaCGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCt                                                        >  5:248971/1‑100 (MQ=255)
                                                                                                                                           aaaCGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCt                                                        <  6:248971/100‑1 (MQ=255)
                                                                                                                                             aCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTg                                                   <  4:2303/103‑1 (MQ=39)
                                                                                                                                             aCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTg                                                   >  3:2303/1‑103 (MQ=39)
                                                                                                                                               gCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTc                <  5:273516/136‑1 (MQ=255)
                                                                                                                                               gCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTc                >  6:273516/1‑136 (MQ=255)
                                                                                                                                                cccGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAACCCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGgcg  >  5:14646/1‑149 (MQ=255)
                                                                                                                                                   gggCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                    <  8:116349/128‑1 (MQ=255)
                                                                                                                                                   gggCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTg                    >  7:116349/1‑128 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCGAAGTGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCG  >  NC_002947/2787476‑2787768

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