Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,313,279 C→G 7.8% intergenic (‑598/+116) PP_16SG ← / ← msrQ 16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,313,2790CG 7.8% 190.6 / 4.9 64intergenic (‑598/+116)PP_16SG/msrQ16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base G (3/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 = 8.78e-01

GCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGC  >  NC_002947/5313133‑5313416
                                                                                                                                                   |                                                                                                                                         
gCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCcgcgc                                                                                                                                         >  3:177638/1‑149 (MQ=255)
 cTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGATGAAATATAAGGGGAAATACCGTATAGGGGCTGCCAAGTGTTAGGCACAACAGATTAAGCGCCTATGAGATCGAGCGCCGCCcgcgc                                                                                                                                         >  1:103089/1‑149 (MQ=14)
    gAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCgcggcg                                                                                                                                     >  4:245282/1‑149 (MQ=255)
         ggggCGCTTTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCGAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                 <  5:57697/148‑1 (MQ=255)
          gggCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAt                                                                                                                               >  2:226459/1‑149 (MQ=255)
                aTTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTcaca                                                                                                                         <  2:103089/149‑1 (MQ=255)
                       ggggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCa                                                                                                                  <  2:187253/149‑1 (MQ=255)
                        gggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCag                                                                                                                 <  7:22175/149‑1 (MQ=255)
                          gATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCa                                                                                                                  >  2:136335/1‑146 (MQ=255)
                               gagaGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                   >  5:50831/1‑124 (MQ=255)
                               gagaGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                   <  6:50831/124‑1 (MQ=255)
                               gagaGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAAc                                                                                                           >  2:158/1‑148 (MQ=255)
                                 gagGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACaaa                                                                                                        <  1:121237/149‑1 (MQ=255)
                                      cGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACATCAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                   <  3:153643/149‑1 (MQ=255)
                                      cGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                   <  6:6444/149‑1 (MQ=255)
                                       gTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCggg                                                                                                  <  1:226459/149‑1 (MQ=255)
                                            ccTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTACCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCctc                                                                                             >  5:325463/1‑149 (MQ=255)
                                             cTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCaga                                                                                                                >  8:241667/1‑129 (MQ=255)
                                             cTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCaga                                                                                                                <  7:241667/129‑1 (MQ=255)
                                             cTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCctct                                                                                            >  6:7976/1‑149 (MQ=255)
                                              tttAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCtctc                                                                                           >  2:129529/1‑149 (MQ=255)
                                              tttAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCtctc                                                                                           >  4:87434/1‑149 (MQ=255)
                                                tAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCtctc                                                                                           <  4:177638/147‑1 (MQ=255)
                                                  aTTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCcaca                                                                                       >  6:57602/1‑149 (MQ=255)
                                                  aTTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCcac                                                                                        <  5:7976/148‑1 (MQ=255)
                                                  aTTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCcac                                                                                        <  3:87434/148‑1 (MQ=255)
                                                      cAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCa                                                                                   >  2:215017/1‑149 (MQ=255)
                                                           aaG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCt                                                                              >  7:67785/1‑149 (MQ=255)
                                                                gAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCTCCTATGAGATCCAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtct                                                                          >  2:139731/1‑149 (MQ=255)
                                                                 aaaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtct                                                                          <  1:136335/148‑1 (MQ=255)
                                                                  aaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtctc                                                                         >  2:326679/1‑148 (MQ=255)
                                                                  aaTATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCa                                                                        <  5:57602/149‑1 (MQ=255)
                                                                   atatAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCaa                                                                       <  8:67785/149‑1 (MQ=255)
                                                                         ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCca                                                                                         >  6:103197/1‑125 (MQ=255)
                                                                         ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCca                                                                                         <  5:103197/125‑1 (MQ=255)
                                                                          gggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGc                                                                <  3:245282/149‑1 (MQ=255)
                                                                          gggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTACCTCAAGACAGGc                                                                >  5:120887/1‑149 (MQ=255)
                                                                              aaTACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                   <  4:344532/78‑1 (MQ=255)
                                                                              aaTACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                   >  3:344532/1‑78 (MQ=255)
                                                                                tACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTgg                                                          >  3:116943/1‑149 (MQ=255)
                                                                                  ccGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTgg                                                          >  6:93518/1‑147 (MQ=255)
                                                                                                       ttCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGc                                                      >  4:238202/1‑130 (MQ=255)
                                                                                                       ttCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGc                                                      <  3:238202/130‑1 (MQ=255)
                                                                                                         cGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCg                                  <  1:215017/148‑1 (MQ=255)
                                                                                                           ccACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTACCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGccc                               <  6:120887/149‑1 (MQ=255)
                                                                                                             acaacaGTTTCAGCTCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCATGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTa                             <  1:70002/149‑1 (MQ=255)
                                                                                                             acaacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTa                             >  8:334799/1‑149 (MQ=255)
                                                                                                               aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTg                                                           >  4:164473/1‑117 (MQ=255)
                                                                                                               aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTg                                                           <  3:164473/117‑1 (MQ=255)
                                                                                                                      tCAGCGCCTATGAGATCGCGCGCCGCCCGCGCGGCGCTCGATCTGCCAGACGCGAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGGTATCTCa                                                                        >  7:249508/1‑97 (MQ=37)
                                                                                                                      tCAGCGCCTAAGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCa                                                                        <  8:249508/97‑1 (MQ=255)
                                                                                                                       cAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTAc                    <  1:139731/148‑1 (MQ=255)
                                                                                                                          cgcCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTgg                >  8:278089/1‑149 (MQ=255)
                                                                                                                          cgcCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTgg                <  1:326679/149‑1 (MQ=255)
                                                                                                                             cTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCaa             >  5:146020/1‑149 (MQ=255)
                                                                                                                               aTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCCCAGGCGCCAGAGAAACACAGGCGGGCCCTCTCCACGGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCGCCCCCCCCCTCTCAAACCGCCCTACCCTGCTACCTGGCAAcc           >  3:205085/1‑149 (MQ=38)
                                                                                                                                tGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAGCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAAc            <  1:158/147‑1 (MQ=255)
                                                                                                                                  agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCTGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTg        <  1:129529/149‑1 (MQ=255)
                                                                                                                                   gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTgg                >  6:214086/1‑140 (MQ=255)
                                                                                                                                   gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTgg                <  5:214086/140‑1 (MQ=255)
                                                                                                                                     tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTggg      >  2:13261/1‑148 (MQ=255)
                                                                                                                                      cGAGCGCCACGCGGGCGGCGCTCGATCTCATAGGCGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACCGGCCCTTGGCAGCCCCCCCCCACTCTCAAGCCGCCCAACCCTGGTACCTGGCAACCTTggg      >  7:157091/1‑147 (MQ=255)
                                                                                                                                       gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTAcc                   >  3:41727/1‑133 (MQ=255)
                                                                                                                                       gAGCGCCGCGCGGGCGGCGCTCGATCTCACACACGCCAGAGAAACACACGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTAcc                   <  4:41727/133‑1 (MQ=39)
                                                                                                                                       gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTAcc                   <  4:36799/133‑1 (MQ=255)
                                                                                                                                        aGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGTCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgc  >  8:180834/1‑149 (MQ=255)
                                                                                                                                               cccgcGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTg        <  6:224956/136‑1 (MQ=255)
                                                                                                                                               cccgcGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTg        >  5:224956/1‑136 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAG‑TGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGC  >  NC_002947/5313133‑5313416

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