Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,877,159 C→A 10.2% A83D (GCC→GAC)  PP_3425 → RND family transporter MFP subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,877,1590CA10.2% 200.4 / 10.1 68A83D (GCC→GAC) PP_3425RND family transporter MFP subunit
Reads supporting (aligned to +/- strand):  ref base C (28/33);  new base A (2/5);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGCGT  >  NC_002947/3877011‑3877304
                                                                                                                                                     |                                                                                                                                                  
gCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGgcgc                                                                                                                                                    <  8:258020/149‑1 (MQ=255)
gCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGgcgc                                                                                                                                                    >  8:293854/1‑149 (MQ=255)
     gcCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCTTCGGTGAAGGCGACCTgg                                                                                                                                               <  1:319000/149‑1 (MQ=255)
      ccGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCTTCGGTGAAGGCGACCTGGt                                                                                                                                              <  3:310324/149‑1 (MQ=255)
      ccGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGt                                                                                                                                              >  8:93948/1‑149 (MQ=255)
      ccGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACAGTCGAAGTGCGCCCGCGCGTCTCCGGCCAGAACGACACGGTCGACTTCAGCGAAGGCGCACAGGt                                                                                                                                              >  6:346645/1‑149 (MQ=255)
         aaGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCaa                                                                                                                                           <  7:73304/149‑1 (MQ=255)
            gTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACAGGCCGCCTCGAAGCCCCGGAGACAGTCGAAGTGCGCCCGCGGTTCTCCGGCCAGATCGAACAGGTCGCCTTCACCGAAGTCGCCCAGGTCAAGa                                                                                                                                         >  3:225677/1‑148 (MQ=255)
             tAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGaaa                                                                                                                                      <  7:179629/149‑1 (MQ=255)
               aGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAg                                                                                                                                     <  5:200636/148‑1 (MQ=255)
               aGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAg                                                                                                                                     <  6:114722/148‑1 (MQ=255)
                   tgGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGAc                                                                                                                                <  2:137704/149‑1 (MQ=255)
                   tgGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGAc                                                                                                                                <  7:93948/149‑1 (MQ=255)
                     gCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCTTCGGTGAAGGCGACCTGGTc                                                                                                                                             <  5:138087/135‑1 (MQ=38)
                     gCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCTGGGCGCCTTCGGTGAAGGCGACCTGGTc                                                                                                                                             >  6:138087/1‑135 (MQ=38)
                         aGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACctgctg                                                                                                                          <  2:21544/149‑1 (MQ=255)
                          ggTGATCGAACAACCGATCACCAAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACctgctg                                                                                                                          >  7:55323/1‑148 (MQ=255)
                              aTCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCa                                                                                                                     <  8:45958/149‑1 (MQ=255)
                                 ggttgttcGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGa                                                                                                                   <  3:257266/141‑1 (MQ=255)
                                 gAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGTCACCGACGAAGTGCGCCCCCGCGTCTCCGGCCCGATCGACCACGTCGCATTCGCCGAAGGCGCGCAGGTCATCGAAGGCTCCGTGCTGGTCCAgg                                                                                                                   >  8:335625/1‑147 (MQ=255)
                                 gAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCTTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGa                                                                                                                   <  7:142788/148‑1 (MQ=255)
                                   acaacCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTGTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATc                                                                                                                 <  1:79818/149‑1 (MQ=255)
                                         gATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCcgc                                                                                                          >  3:63615/1‑149 (MQ=255)
                                              ccGACTGGGACGAGTGCACCGGCCGCCTCGAATCCCCGGAAACCTTCGAAGTGCGCCCGCGTGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTc                                                                                                                                             <  6:167849/110‑1 (MQ=255)
                                              ccGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGACACCCACGAAGTGCGCCCGCGGGTCGGCGGCGCGATCGACCCGCTCCCAGTCACCGACGTCGCCGACGTCAAGAAACGCGACCGCCTGTTCCAGATAGACCCGCGACCg                                                                                                     >  6:213420/1‑149 (MQ=255)
                                              ccGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCg                                                                                                     >  3:347220/1‑149 (MQ=255)
                                              ccGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGACAGAGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTc                                                                                                                                             >  5:167849/1‑110 (MQ=255)
                                                 aaTGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTc                                                                                                  <  7:293854/149‑1 (MQ=255)
                                                     ggACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAgg                                                                                              <  8:55323/149‑1 (MQ=255)
                                                          aGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGaa                                                                                         <  2:145336/149‑1 (MQ=255)
                                                              cACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTcc                                                                                     >  4:315445/1‑149 (MQ=255)
                                                               aCCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTg                                                                                           >  3:183860/1‑142 (MQ=255)
                                                               aCCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTg                                                                                           <  4:183860/142‑1 (MQ=255)
                                                                ccGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCCTCACCGAAGGCGCCCAGGTCAAGAAACGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTccgc                                                                                   <  6:163061/149‑1 (MQ=255)
                                                                  ggCCGCCTCGAAGCCCCGGAACCCCTCGAAGTGCGCCCGCGGCTCGCCGGCCAGATCGACCAGGTCGGGGTCCCCGAAGGCGCCCAGGGCAAGCAACGCGCCCTGCAGTTCCAGCACGACCCGCGCCCTGTCCAGGCTGAAGAcctccg                                                                                 >  2:221526/1‑149 (MQ=255)
                                                                    ccgccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCgccgcc                                                                               <  3:315445/149‑1 (MQ=255)
                                                                     cgccTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGCCCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCt                                                                              >  5:222134/1‑149 (MQ=255)
                                                                           gAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGTTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCAGTTCCAGGATGAATTCCGCCGCCTCGaa                                                                          <  2:146346/147‑1 (MQ=255)
                                                                           gAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGaa                                                                          <  4:88718/147‑1 (MQ=255)
                                                                           gAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGc                                                                        >  2:115552/1‑149 (MQ=255)
                                                                                 ccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCg                                                                                                                                  >  5:341983/1‑85 (MQ=255)
                                                                                 ccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCg                                                                                                                                  <  6:341983/85‑1 (MQ=255)
                                                                                 ccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGaa                                                                          <  4:3252/141‑1 (MQ=255)
                                                                                 ccGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGaa                                                                          >  3:3252/1‑141 (MQ=255)
                                                                                       cgggcataaGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACa                                                            <  3:53087/142‑1 (MQ=255)
                                                                                              aaGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaag                                                     >  8:325937/1‑149 (MQ=255)
                                                                                                gTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaaggc                                                   >  8:251411/1‑149 (MQ=255)
                                                                                                    gcCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTTAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCAcc                                               <  1:44023/149‑1 (MQ=255)
                                                                                                         cgGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGCCCCGCGCCCGTTCCAGGCGGAAGTCCGCCGCCTCGAAGCCCAACTGCAACag                                                           >  1:91652/1‑132 (MQ=255)
                                                                                                          gggTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGAGTCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATc                                         <  4:63615/149‑1 (MQ=255)
                                                                                                             tctcCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCg                                       >  3:319169/1‑148 (MQ=255)
                                                                                                                ccGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAgc                                   >  2:77787/1‑149 (MQ=255)
                                                                                                                ccGGCCAGATCGACCAGGTCACCTTCACCGAAGGCGCCCAGGTCAAGAAGAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAg                                    >  1:85214/1‑149 (MQ=255)
                                                                                                                   gCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATc                                         >  4:189511/1‑140 (MQ=255)
                                                                                                                   gCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATc                                         <  3:189511/140‑1 (MQ=255)
                                                                                                                     cAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGTTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCCCGCCCGTTCCAGGCTTAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCaa                              <  5:346645/149‑1 (MQ=255)
                                                                                                                            aCCTGGGCGCCTTCGGTGAAGGCGACCTGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAgc                                   >  4:112416/1‑137 (MQ=255)
                                                                                                                            aCCTGGGCGCCTTCGGTGAAGGCGACCTGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAgc                                   <  3:112416/137‑1 (MQ=255)
                                                                                                                            aCCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGcc                       >  3:39553/1‑149 (MQ=255)
                                                                                                                                    gCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTgg               <  1:77787/149‑1 (MQ=255)
                                                                                                                                      cTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaa                                                      >  1:178180/1‑108 (MQ=255)
                                                                                                                                      cTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACaggccaa                                                      <  2:178180/108‑1 (MQ=255)
                                                                                                                                      cTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGc              >  2:211915/1‑148 (MQ=255)
                                                                                                                                       ttCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGc              >  8:58799/1‑147 (MQ=255)
                                                                                                                                          cACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAg           <  1:211915/148‑1 (MQ=255)
                                                                                                                                           aCCGAAGGCGCCCATGTCAAGAAAGGCGACCTGCTGTTCCAGATCGAGCGGCGCGATTTGCAGGCTGAATTCGGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAgc          <  2:87953/148‑1 (MQ=255)
                                                                                                                                           aCCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAgcg         <  7:251411/149‑1 (MQ=255)
                                                                                                                                               aaGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTg     >  4:114031/1‑149 (MQ=255)
                                                                                                                                                aGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGc    <  7:58799/149‑1 (MQ=255)
                                                                                                                                                aGGCGACCAGGTCAAGAAAGTCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGc    <  6:222134/149‑1 (MQ=255)
                                                                                                                                                  gcgcCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGCGt  >  2:192789/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCTGCGCCGAAGGTAAGTGTGGCCAAGGTGATCGAACAACCGATCACCGAATGGGACGAGTTCACCGGCCGCCTCGAAGCCCCGGAAACCGTCGAAGTGCGCCCGCGGGTCTCCGGCCAGATCGACCAGGTCGCCTTCACCGAAGGCGCCCAGGTCAAGAAAGGCGACCTGCTGTTCCAGATCGACCCGCGCCCGTTCCAGGCTGAAGTCCGCCGCCTCGAAGCCCAACTGCAACAGGCCAAGGCCACCGCCATCCGCAGCGCCAACGAAGCCCGCCGTGGCGAGCGCCTGCGT  >  NC_002947/3877011‑3877304

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: