Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,215,672 C→G 8.0% L529L (CTC→CTG) ‡ mexD → multidrug RND transporter MexD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,215,6720CG 8.0% 190.4 / 3.8 61L529L (CTC→CTG) ‡mexDmultidrug RND transporter MexD
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base G (2/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGCCCTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTGGATGTGCAGTTGC  >  NC_002947/3215537‑3215817
                                                                                                                                       |                                                                                                                                                   
ccTGGCCCTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGTTGAGCGGCGaatag                                                                                                                                         <  2:273588/148‑4 (MQ=255)
 cTGGCCCTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACtggt                                                                                                                                       <  1:203285/149‑1 (MQ=255)
   ggCCCTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACtggtgg                                                                                                                                    >  6:222532/1‑149 (MQ=255)
      ccTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGTTGAGCAGCGaatagcgtt                                                                                                                                     <  7:231879/146‑8 (MQ=255)
      ccTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGTTGAGCAGCGaatagcgtt                                                                                                                                     >  8:231879/1‑139 (MQ=255)
      ccTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTc                                                                                                                                                     >  4:2500/1‑130 (MQ=255)
      ccTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTc                                                                                                                                                     <  3:2500/130‑1 (MQ=255)
                  cccGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACa                                                                                                                                                 <  8:61835/122‑1 (MQ=255)
                  cccGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACa                                                                                                                                                 >  7:61835/1‑122 (MQ=255)
                  cccGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGaata                                                                                                                                          >  7:175505/1‑127 (MQ=255)
                  cccGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGaata                                                                                                                                          <  8:175505/129‑3 (MQ=255)
                   ccGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTg                                                                                                                   >  5:305562/1‑149 (MQ=255)
                    cGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTgg                                                                                                                  <  4:305320/149‑1 (MQ=255)
                          tgtgCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTg                                                                                                            >  8:93960/1‑149 (MQ=255)
                              cgcCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTgccggcc                                                                                                                       <  7:25448/134‑1 (MQ=255)
                              cgcCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTgccggcc                                                                                                                       >  8:25448/1‑134 (MQ=255)
                              cgcCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCg                                                                                                         >  7:99037/1‑148 (MQ=255)
                              cgcCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCg                                                                                                         >  6:284185/1‑148 (MQ=255)
                                ccACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTAt                                                                                                      >  3:78234/1‑149 (MQ=255)
                                  aCGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTAtgt                                                                                                    <  2:59537/149‑1 (MQ=255)
                                  aCGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTAtgt                                                                                                    <  5:62739/149‑1 (MQ=255)
                                    gctgctCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTgg                                                                                                  <  7:93960/149‑1 (MQ=255)
                                     ctgctCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGCGCACGTGCCGGCCGCTGGATGCTGGCGTATGTgg                                                                                                  <  1:44417/149‑1 (MQ=255)
                                        ctCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGt                                                                                                        >  3:63396/1‑139 (MQ=255)
                                        ctCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGt                                                                                                        <  4:63396/139‑1 (MQ=255)
                                         tCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATc                                                                                             >  3:276910/1‑149 (MQ=255)
                                          cAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATcc                                                                                            >  4:110570/1‑149 (MQ=255)
                                           aaGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATc                                                                                             <  8:99037/147‑1 (MQ=255)
                                           aaGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCt                                                                                           >  5:202427/1‑149 (MQ=255)
                                               cGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCtggtg                                                                                       <  5:284185/149‑1 (MQ=255)
                                                     cccATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTg                                                                                 <  6:305562/149‑1 (MQ=255)
                                                      ccATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTgg                                                                                >  4:267044/1‑149 (MQ=255)
                                                       cATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTggg                                                                               >  4:28920/1‑149 (MQ=255)
                                                            ccaccaTGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTAtt                                                                          <  4:276910/149‑1 (MQ=255)
                                                             caccaTGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTc                                                                         >  8:115054/1‑149 (MQ=255)
                                                               ccaTGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTcc                                                                        >  6:183526/1‑148 (MQ=255)
                                                                caTGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTa                                                                      >  6:187001/1‑149 (MQ=255)
                                                                    aaaaGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCt                                                                   <  7:115054/148‑1 (MQ=255)
                                                                     aaaGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTg                                                                  <  4:120763/148‑1 (MQ=255)
                                                                       aGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTgcgc                                                               >  2:193098/1‑149 (MQ=255)
                                                                                          cTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGACATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTgc                                            >  1:122938/1‑149 (MQ=255)
                                                                                           ttCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCAGTGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTgc                                            >  1:258024/1‑149 (MQ=255)
                                                                                           ttCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTgcc                                           >  6:218390/1‑149 (MQ=255)
                                                                                            tCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTgccg                                          >  2:223135/1‑149 (MQ=255)
                                                                                            tCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGAGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTgccg                                          >  5:217279/1‑149 (MQ=255)
                                                                                                  gTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAg                                    <  5:183526/149‑1 (MQ=255)
                                                                                                     gCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGAc                                 <  5:187001/149‑1 (MQ=255)
                                                                                                      cTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCtt                                                             <  1:74022/120‑1 (MQ=255)
                                                                                                      cTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCtt                                                             >  2:74022/1‑120 (MQ=255)
                                                                                                                gtCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTAc                                                                     >  4:241585/1‑102 (MQ=255)
                                                                                                                gtCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTAc                                                                     <  3:241585/102‑1 (MQ=255)
                                                                                                                    ccGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGcgtc                  >  6:256745/1‑149 (MQ=255)
                                                                                                                      gAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGACATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGcgtcgt                <  2:122938/149‑1 (MQ=255)
                                                                                                                        aCGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTgg              <  5:222532/149‑1 (MQ=255)
                                                                                                                        aCGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTgg              <  2:258024/149‑1 (MQ=255)
                                                                                                                        aCGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTgg              <  3:81965/149‑1 (MQ=255)
                                                                                                                         cGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTGGa             >  3:246462/1‑149 (MQ=255)
                                                                                                                          gCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTGGAt            <  3:267044/149‑1 (MQ=255)
                                                                                                                            agttCGCTGTTGAGCAGCGAAATGGTGTCACGTGCCGGCCGCTGGATTCTGGCGTATGTGTGCATCCTGGTGGTGGTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTGTAGCGTCGTgg              <  6:16532/143‑1 (MQ=255)
                                                                                                                            agttCGCTGTTGACCAGCGAACTGGTTGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCATTCGTGCCGGCCGAGGACCTTGGCTGTAGAGTCGTGGaag           >  5:16532/3‑146 (MQ=255)
                                                                                                                                 gctgctCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTGGATGTGCAGt     >  5:251816/1‑149 (MQ=255)
                                                                                                                                     ctCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTGGATGTGCAGTTg   <  3:110570/147‑1 (MQ=255)
                                                                                                                                     ctCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTGGATGTGCAGTTGc  >  6:273523/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
CCTGGCCCTGACCTTCACCCCGGCGCTGTGCGCCACGCTGCTCAAGCCGGTGCCCCATGGCCACCATGAAAAGGCCGGTTTCTTCGGCGCCTTCAACCGTGGCTTCGCCCGTGTCACCGAACGCTATTCGCTGCTCAACAGCGAACTGGTGGCACGTGCCGGCCGCTGGATGCTGGCGTATGTGGGCATCCTGGTGGTGCTGGGCTATTCCTACCTGCGCTTGCCGGAGGCCTTCGTGCCGGCCGAGGACCTTGGCTATAGCGTCGTGGATGTGCAGTTGC  >  NC_002947/3215537‑3215817

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: