Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,350,523 C→G 6.0% G910G (GGC→GGG)  PP_2065 → multidrug RND transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,350,5230CG 6.0% 217.0 / 3.4 67G910G (GGC→GGG) PP_2065multidrug RND transporter
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

GCTGGCGGCGCAGTACGAAAGCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGACTTTTG  >  NC_002947/2350376‑2350668
                                                                                                                                                   |                                                                                                                                                 
gCTGGCGGCGCAGTACGAAAGCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGcc                                                                                                                                                  <  4:372443/149‑1 (MQ=255)
 cTGGCGGCGCAGTACGAAAGCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCt                                                                                                                                                 <  7:186980/149‑1 (MQ=255)
 cTGGCGGCGCAGCACGAAAGCCTGCTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCt                                                                                                                                                 <  2:305171/149‑1 (MQ=255)
      ggcgCAGTACGAAAGCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGAt                                                                                                                                              >  2:61587/1‑147 (MQ=255)
        cgcAGTACGAAAGCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCAGGCCGATCagc                                                                                                                                          <  5:262670/149‑1 (MQ=255)
                    gCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATc                                                                                                                              <  7:7978/149‑1 (MQ=255)
                    gCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATc                                                                                                                              <  7:2111/149‑1 (MQ=255)
                         gTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGAt                                                                                                                         >  3:206453/1‑149 (MQ=255)
                            gACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGAtgg                                                                                                                       <  8:116832/148‑1 (MQ=255)
                            gACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGAtgg                                                                                                                       <  4:206453/148‑1 (MQ=255)
                               acccTGATCATTCTGGTGAACGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCATGATGGTGGc                                                                                                                   >  8:291695/2‑149 (MQ=255)
                                 ccTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGAGCAGCAAGCACGGTATCCTGATGGTGGCGt                                                                                                                 >  8:103291/1‑149 (MQ=255)
                                        aTTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCaa                                                                                                          >  5:62674/1‑149 (MQ=255)
                                        aTTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCaa                                                                                                          >  6:207811/1‑149 (MQ=255)
                                         ttCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAAc                                                                                                         <  5:376864/149‑1 (MQ=255)
                                         ttCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAAc                                                                                                         >  5:202580/1‑149 (MQ=255)
                                         ttCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAAc                                                                                                         >  5:16433/1‑149 (MQ=255)
                                            tGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCagcaag                                                                                                                                       <  7:287975/116‑1 (MQ=255)
                                            tGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCagcaag                                                                                                                                       >  8:287975/1‑116 (MQ=255)
                                              gTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACgaga                                                                                                     <  6:16433/148‑1 (MQ=255)
                                              gTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGAt                                                                                                    >  1:11873/1‑149 (MQ=255)
                                                 aCCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCa                                                                                                 <  1:336242/149‑1 (MQ=255)
                                                 aCCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACACTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGAtt                                                                                                   >  1:154690/1‑147 (MQ=255)
                                                     tGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACgaga                                                                                                     >  8:270277/1‑141 (MQ=255)
                                                     tGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACgaga                                                                                                     <  7:270277/141‑1 (MQ=255)
                                                        cGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTgcgc                                                                                          <  5:292115/149‑1 (MQ=255)
                                                         gCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGcgcg                                                                                         <  1:211159/149‑1 (MQ=255)
                                                          cTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGa                                                                                        >  6:259323/1‑149 (MQ=255)
                                                             tCGATCTGCGGTGCGCTGCTGCTTTTGTTTATGGTTTGGGCCACGATGAAGATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCg                                                                                                              <  2:15485/124‑1 (MQ=25)
                                                             tCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCg                                                                                                              >  1:15485/1‑124 (MQ=255)
                                                               gATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAAcc                                                                                   >  4:132387/1‑149 (MQ=255)
                                                                   tGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCt                                                                                  >  5:6999/1‑146 (MQ=255)
                                                                   tGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCt                                                                                  <  6:6999/146‑1 (MQ=255)
                                                                     cGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATc                                                                             <  5:207811/149‑1 (MQ=255)
                                                                       gTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGt                                                                           >  4:364118/1‑149 (MQ=255)
                                                                              gctgcCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTgccg                                                                    >  1:260915/1‑149 (MQ=255)
                                                                               ctgcCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTgccgc                                                                   >  4:94209/1‑149 (MQ=255)
                                                                                  ccGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCATGTGCGCGACAACCTCGATCGTGCTGCCGCCAt                                                                <  6:62674/149‑1 (MQ=255)
                                                                                   cgctggcgctggGCTGGGCCACGCTGACCATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCAtt                                                               <  4:55285/149‑1 (MQ=255)
                                                                                       ggcgctggGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTc                                                           <  5:259323/149‑1 (MQ=255)
                                                                                         cgctggGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGt                                                         <  3:364118/149‑1 (MQ=255)
                                                                                          gctgggctggGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTg                                                        >  4:44853/1‑149 (MQ=255)
                                                                                           cttggctgggcGACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTgc                                                       <  8:236445/146‑1 (MQ=255)
                                                                                              ggctgggcCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTgcggc                                                    <  2:260915/149‑1 (MQ=255)
                                                                                                  gggcCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCa                                                 >  2:253282/1‑148 (MQ=255)
                                                                                                    gcCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCaaa                                               >  8:20009/1‑148 (MQ=255)
                                                                                                    gcCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAAt                                              <  1:61587/149‑1 (MQ=255)
                                                                                                      ccgacTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATAc                                            <  2:319777/145‑1 (MQ=255)
                                                                                                             aaCATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTg                                       <  3:132387/147‑1 (MQ=255)
                                                                                                              aCATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGt                                                             <  2:65452/124‑1 (MQ=255)
                                                                                                              aCATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGt                                                             >  1:65452/1‑124 (MQ=255)
                                                                                                                    acacacAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCa                                 <  1:403991/146‑1 (MQ=255)
                                                                                                                    acacacAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCa                                 >  2:403991/1‑146 (MQ=255)
                                                                                                                     cacacaGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCa                                 <  6:202580/145‑1 (MQ=255)
                                                                                                                      acacaGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCt                            >  5:235004/1‑149 (MQ=255)
                                                                                                                        acaGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCtga                          >  3:59832/1‑149 (MQ=255)
                                                                                                                        acaGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGGGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCtga                          <  8:287769/149‑1 (MQ=255)
                                                                                                                            aTTGGCCAGGTCACGCTGAACCGCCTGATCACCAAGCACGCTAACCTGATCGTGCCCATCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATAGTTCGTGCGGCGCAAATCACCAGGCGGCCAGTGCTGAAGAc                      >  4:272205/1‑149 (MQ=255)
                                                                                                                                 ccTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACGAcgg                 >  4:87437/1‑149 (MQ=255)
                                                                                                                                   tGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGTTGATGACCAcggcg               >  8:385926/1‑149 (MQ=255)
                                                                                                                                      gcttgCTGATCAGGCCGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACggcggc             >  8:70645/5‑148 (MQ=255)
                                                                                                                                        aCGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGAt          <  3:94209/149‑1 (MQ=255)
                                                                                                                                          gCTGATCAGGCCGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCt                                                                                  <  2:123311/75‑1 (MQ=37)
                                                                                                                                          gCTGATCAGGCCGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCt                                                                                  >  1:123311/1‑75 (MQ=37)
                                                                                                                                           cTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGAc       >  6:42184/1‑149 (MQ=255)
                                                                                                                                               tCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGACtttt   >  2:208130/1‑149 (MQ=255)
                                                                                                                                                 ggCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGACTTTTg  <  3:44853/148‑1 (MQ=255)
                                                                                                                                                 ggCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGACTTTTg  <  3:256194/148‑1 (MQ=255)
                                                                                                                                                   ccTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTg                              >  6:321249/1‑118 (MQ=255)
                                                                                                                                                   ccTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTg                              <  5:321249/118‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCTGGCGGCGCAGTACGAAAGCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGACTTTTG  >  NC_002947/2350376‑2350668

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