Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,044,479 G→T 5.6% intergenic (+20/+29) PP_3565 → / ← PP_3566 LysE family transporter/MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,044,4790GT 5.6% 231.2 / 2.7 71intergenic (+20/+29)PP_3565/PP_3566LysE family transporter/MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base T (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

ATACCTGCTGGCGCGCCTGGCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCAC  >  NC_002947/4044336‑4044618
                                                                                                                                               |                                                                                                                                            
aTACCTGCTGGCGCGCCTGGCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGcc                                                                                                                                         >  3:238325/1‑149 (MQ=255)
   ccTGCTGGCGCGCCTGGCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCTTCGCCCGt                                                                                                                                      <  2:208058/149‑1 (MQ=255)
   ccTGCTGGCGCGCCTGGCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGt                                                                                                                                      >  7:315745/1‑149 (MQ=255)
       ctgGCGCGCCTGGCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGccc                                                                                                                                  >  6:45091/1‑149 (MQ=255)
         ggCGCGCCTGGCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTg                                                                                                                                >  2:322892/1‑149 (MQ=255)
                  ggCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGt                                                                                                                        >  4:3550/1‑148 (MQ=255)
                   gCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTaa                                                                                                                     <  8:340160/149‑1 (MQ=255)
                    cGTTTCGCATTCGGCCATGGCTGTCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGATACTGTTTGCCTTGATTGTTTTGGCGCTGCTGCTTGGACGTTAGCGGGCCCGCGAACACGGTGGAAGCCCGTGCCCTGCACCGCGTat                                                                                                                      <  4:387159/148‑2 (MQ=255)
                     gTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACg                                                                                                                   >  7:179632/1‑149 (MQ=255)
                     gTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACg                                                                                                                   <  7:178062/149‑1 (MQ=255)
                     gTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACg                                                                                                                   <  1:84266/149‑1 (MQ=255)
                     gTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACg                                                                                                                   <  1:265584/149‑1 (MQ=255)
                          gCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTc                                                                                                               <  4:238325/148‑1 (MQ=255)
                          gCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTc                                                                                                               >  6:153642/1‑148 (MQ=255)
                            aTTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAAc                                                                                                            <  2:187031/149‑1 (MQ=255)
                              tCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAAccc                                                                                                          >  2:268565/1‑149 (MQ=255)
                                ggCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACcgc                                                                                                                          <  4:194672/132‑1 (MQ=255)
                                ggCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACcgc                                                                                                                          >  3:194672/1‑132 (MQ=255)
                                    aTGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGCTTTGGCGCGGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGt                                                                                                    <  2:131743/149‑1 (MQ=255)
                                           gcCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCTCCGCGTAACGGTTCAACCCCGGCGTGCttt                                                                                               <  1:34647/147‑1 (MQ=255)
                                           gcCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCttt                                                                                               <  5:153642/147‑1 (MQ=255)
                                           gcCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGa                                                                                             <  4:310803/149‑1 (MQ=255)
                                            ccAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCttt                                                                                               >  6:169211/1‑146 (MQ=255)
                                            ccAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGa                                                                                             <  5:169211/148‑1 (MQ=255)
                                            ccAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAg                                                                                            <  3:3550/149‑1 (MQ=255)
                                                gggTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTAGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGa                                                                                        <  1:322892/149‑1 (MQ=255)
                                                   tGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAac                                                                                     <  1:268565/149‑1 (MQ=255)
                                                   tGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAac                                                                                     <  8:315745/149‑1 (MQ=255)
                                                     gCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAacac                                                                                   <  5:45091/149‑1 (MQ=255)
                                                      cAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAg                                                                                         <  7:126310/142‑1 (MQ=255)
                                                      cAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAg                                                                                         >  8:126310/1‑142 (MQ=255)
                                                       aaaGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCttt                                                                                               <  4:88307/135‑1 (MQ=255)
                                                       aaaGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCttt                                                                                               >  3:88307/1‑135 (MQ=255)
                                                       aaaGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACaa                                                                                 >  8:17946/1‑149 (MQ=255)
                                                       aaaGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACaa                                                                                 >  7:4284/1‑149 (MQ=255)
                                                       aaaGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCCGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACaa                                                                                 >  2:146093/1‑149 (MQ=255)
                                                        aaGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAAc                                                                                <  4:243698/149‑1 (MQ=255)
                                                            tttCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAgc                                                                            >  2:86315/1‑149 (MQ=255)
                                                                aaTCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCca                                                                        <  4:183106/149‑1 (MQ=255)
                                                                 aTCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCca                                                                        >  5:213783/1‑148 (MQ=255)
                                                                  tCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCcaca                                                                      <  5:229703/149‑1 (MQ=255)
                                                                       gctgcGGGTGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCgg                                                                 >  5:17702/1‑149 (MQ=255)
                                                                                 ccTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGcc                                                       >  6:37163/1‑149 (MQ=255)
                                                                                  cTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGccc                                                      >  1:139862/1‑149 (MQ=255)
                                                                                       tGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCccac                                                  <  2:139862/148‑1 (MQ=255)
                                                                                        gCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCCACCCCGccccccccc                                                >  5:383956/1‑148 (MQ=255)
                                                                                        gCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCccccca                                                >  3:373198/1‑149 (MQ=255)
                                                                                        gCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCccac                                                  >  2:241224/1‑147 (MQ=255)
                                                                                         ccTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCccacca                                                >  4:145348/1‑148 (MQ=255)
                                                                                             gATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAACCACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGc                                            >  3:215502/1‑148 (MQ=255)
                                                                                                  gTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACACGCCAACCCCCCCCCCACCAGGTCCACGc                                       >  4:254061/1‑148 (MQ=255)
                                                                                                   tttGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTATAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCt                                     >  3:239863/1‑149 (MQ=255)
                                                                                                   tttGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCCCCAGTGCCACGCCt                                     >  4:291966/1‑149 (MQ=255)
                                                                                                     tGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTg                                   <  3:67113/149‑1 (MQ=255)
                                                                                                      ggCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGGCCCCCCCCGGGCCCCGccgttc                                  >  8:247251/1‑145 (MQ=255)
                                                                                                      ggCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTg                                   <  3:291966/148‑1 (MQ=255)
                                                                                                       gcgcTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTgcc                                 <  3:145348/149‑1 (MQ=255)
                                                                                                       gcgcTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTgc                                  <  1:241224/148‑1 (MQ=255)
                                                                                                         gcTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTgccag                               >  7:161087/1‑149 (MQ=255)
                                                                                                                        aTAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAACCACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAAc                <  4:215502/149‑1 (MQ=255)
                                                                                                                           gCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAgc             <  8:161087/149‑1 (MQ=255)
                                                                                                                              ggCCCGCGAACACGGGCGAAGCGCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCa          <  8:179632/149‑1 (MQ=255)
                                                                                                                                cccGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCCCACCCGGCAACCCCGCCCCCACCCGTGGCACGCCTTGCCCGCCCAACGGGCTGAACAGCGCCCt        >  7:219956/1‑149 (MQ=255)
                                                                                                                                 ccGCGAACACGGGCGAAGCCCGTGTCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACACCGCAAACCCGCCCCCACCACTGCCACCCCTTGCCCTCCAAACAGGCTGCACAGCGCACTg       >  4:219700/1‑149 (MQ=255)
                                                                                                                                 ccGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTg       <  6:383293/149‑1 (MQ=255)
                                                                                                                                 ccGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTg       <  5:37163/149‑1 (MQ=255)
                                                                                                                                 ccGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTg       <  8:4284/149‑1 (MQ=255)
                                                                                                                                   gcgAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGAGCTTTGAGTAGAAACACAACCAGCTCCACACCCGGCAACCCCCCCCCCCCCAGTGCCACGCCTTTCCAGCCCAACTGGCGGAACAGCGCCCTGGc     >  7:251064/1‑149 (MQ=255)
                                                                                                                                      ggcacGGGCTTCGCCCGTGCCCTGCACCGCGTATACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTg       >  1:252456/3‑145 (MQ=255)
                                                                                                                                      ggcacGGGCTTCGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCCCGCCTTGCCCGCCAAACGGGCTGAACAGCGCACGTGCCAc  >  8:144681/3‑149 (MQ=255)
                                                                                                                                      ggcacGGGCTTCGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTg       <  2:252456/142‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
ATACCTGCTGGCGCGCCTGGCGTTTCGCATTCGGCCATGGCTGGCCAAGGGTGGCAAAGGTTTCAATCGCTGCTGCGGGAGCCTGTTTGCCTTGATTGGTTTGGCGCTGCCGCTTGGCCGATAGCGGGCCCGCGAACACGGGCGAAGCCCGTGCCCTGCACCGCGTAACGGTTCAACCCCGGCGTGCTTTGAGTAGAAACACAACCAGCGCCACACCCGGCAACCCCGCCCCCACCAGTGCCACGCCTTGCCAGCCAAACTGGCTGAACAGCGCACTGGCCAC  >  NC_002947/4044336‑4044618

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