Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,350,525 T→C 6.0% L911P (CTG→CCG)  PP_2065 → multidrug RND transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,350,5250TC 6.0% 196.2 / 2.5 67L911P (CTG→CCG) PP_2065multidrug RND transporter
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

CTGGCGGCGCAGTACGAAAGCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGACTTTTGGCGTG  >  NC_002947/2350377‑2350673
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                               >  6:207811/1‑149 (MQ=255)
                                       aTTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCaa                                                                                                               >  5:62674/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                                                                                                                                            >  8:287975/1‑116 (MQ=255)
                                           tGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCagcaag                                                                                                                                            <  7:287975/116‑1 (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                                                                                                          <  7:270277/141‑1 (MQ=255)
                                                    tGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACgaga                                                                                                          >  8:270277/1‑141 (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)
                                                                                                                                                    tGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTCCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGACTTTTGGCGtg  >  5:184128/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGGCGGCGCAGTACGAAAGCCTGGTCGACCCCCTGATCATTCTGGTGACCGTGCCGCTGTCGATCTGCGGTGCGCTGCTGCCGCTGGCGCTGGGCTGGGCCACGCTGAACATCTACACACAGATTGGCCTGGTGACGCTGATCGGCCTGATCAGCAAGCACGGTATCCTGATGGTGGCGTTCGCCAACGAGATTCAGGTGCGCGACAACCTCGATCGTGCTGCCGCCATTGTTCGTGCGGCGCAAATACGCTTGCGGCCAGTGCTGATGACCACGGCGGCGATGACTTTTGGCGTG  >  NC_002947/2350377‑2350673

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: