Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,069,170 G→T 6.0% L997M (CTG→ATG)  mdtB ← multidrug transporter membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,069,1700GT 6.0% 221.3 / 3.0 67L997M (CTG→ATG) mdtBmultidrug transporter membrane protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base T (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GTCGCCACCGGGCGACGAATGAAAGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGG  >  NC_002947/4069023‑4069316
                                                                                                                                                   |                                                                                                                                                  
gTCGCCACCGGGCGACGAATGAAAGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAg                                                                                                                                                    <  8:238629/148‑1 (MQ=255)
       gcGGGCGACGAATGAAAGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACccgc                                                                                                                                            <  4:286655/148‑1 (MQ=255)
        cGGGCGACGAATGAAAGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACccgcc                                                                                                                                           <  2:131095/149‑1 (MQ=255)
         gggCGACGAATGAAAGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACccgccg                                                                                                                                          <  4:131177/149‑1 (MQ=255)
            cgacgaATGAAAGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGAc                                                                                                                                        >  8:129910/1‑148 (MQ=255)
                     aaaGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAgg                                                                                                                              <  7:129910/149‑1 (MQ=255)
                     aaaGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAg                                                                                                                               <  1:369041/148‑1 (MQ=255)
                         gCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGcc                                                                                                                            >  3:132790/1‑147 (MQ=255)
                             ggACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGtgg                                                                                                                      <  2:246998/149‑1 (MQ=255)
                                   gTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGgcg                                                                                                                <  1:249005/149‑1 (MQ=255)
                                   gTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGgcg                                                                                                                >  2:238107/1‑149 (MQ=255)
                                   gTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTCACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGgcg                                                                                                                >  5:355283/1‑149 (MQ=255)
                                     tCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGgcg                                                                                                                >  3:70894/1‑147 (MQ=255)
                                      cATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAg                                                                                                             >  2:119081/1‑149 (MQ=255)
                                         gCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTc                                                                                                          >  5:387870/1‑149 (MQ=255)
                                           cTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCAt                                                                                                                                     >  1:360900/1‑120 (MQ=255)
                                           cTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCAt                                                                                                                                     <  2:360900/120‑1 (MQ=255)
                                              ggCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGgcgc                                                                                                     <  1:238107/149‑1 (MQ=255)
                                                 cTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGccc                                                                                                   >  3:400686/1‑148 (MQ=255)
                                                    cTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGaa                                                                                                <  4:132790/148‑1 (MQ=255)
                                                    cTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAAc                                                                                               <  2:262128/149‑1 (MQ=255)
                                                    cTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAAc                                                                                               >  6:262712/1‑149 (MQ=255)
                                                     ttCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAAcc                                                                                              >  7:38022/1‑149 (MQ=255)
                                                     ttCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACACGTAGATGACCCGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCACGGCCCCCGAAcc                                                                                              >  2:316706/1‑149 (MQ=255)
                                                       cACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCg                                                                                             <  1:119081/148‑1 (MQ=255)
                                                         cGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTg                                                                                          <  5:262712/149‑1 (MQ=255)
                                                          gTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTgg                                                                                         >  8:237282/1‑149 (MQ=255)
                                                          gTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTgg                                                                                         >  5:343213/1‑149 (MQ=255)
                                                           tCAGTGGCCGGGCGCAGCCGGCGGGCCAGGCGCCCGAAGTACCGGCCGCCGACCGGACTCGTGACCCGCCTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCCCCCGGCCCAGTGGCTGGCGCAGCGCGGCGCCCGAACCGGTg                                                                                          <  2:310247/147‑1 (MQ=255)
                                                             aGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTg                                                                                          >  1:310247/1‑145 (MQ=255)
                                                             aGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGc                                                                                        >  6:148812/1‑147 (MQ=255)
                                                              gTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGATCAGCGTCAGCTCCTGGCTGACCATCTACCCGCCGACGATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAg                                                                                     >  7:271728/1‑149 (MQ=255)
                                                              gTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCATCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAg                                                                                     >  7:277295/1‑149 (MQ=255)
                                                                  ccGGGCGCAGGCGGCGGGCCAGGTGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATc                                                                                 <  2:274326/149‑1 (MQ=255)
                                                                              ggcggGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGgcaccgca                                                                     <  4:70894/149‑1 (MQ=255)
                                                                              ggcggGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAAGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGgcaccgca                                                                     >  6:413966/1‑149 (MQ=255)
                                                                                  ggCCAGGCGATCGAAGTACAGGTACATGACCGGAGTGGTGAACAGCGTCAGCATCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCa                                                                  >  2:53791/1‑148 (MQ=255)
                                                                                          gATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGCCCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGcgcg                                                         >  3:16204/1‑149 (MQ=255)
                                                                                          gATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGcgcg                                                         >  5:310311/1‑149 (MQ=255)
                                                                                           aTCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCgg                                                        <  4:400686/149‑1 (MQ=255)
                                                                                                aGTACAGGTAGATGACCGGAGTGGTGTACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTGGGCGCCCGAACCGTTGGCAAGCATCAGCGGCACCGCACCAAACAGCGAGGCCAgg                                                   <  1:316706/149‑1 (MQ=255)
                                                                                                  tACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGggt                                                 <  6:310311/149‑1 (MQ=255)
                                                                                                    cAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGggtgg                                               >  2:278484/1‑149 (MQ=255)
                                                                                                     aGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGggtgg                                               >  1:349541/1‑148 (MQ=255)
                                                                                                       gTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGtca                                            >  1:276153/1‑149 (MQ=255)
                                                                                                          gaagACCTGAGTGGTGAACAGCGTCAGCACCTGTCTGACCATCAACCCGCCTACCATCAACACGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGtcatca                                         <  8:271728/146‑1 (MQ=255)
                                                                                                          gATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCATCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGtcatca                                         <  8:277295/149‑1 (MQ=255)
                                                                                                            tGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAgg                                       <  5:176105/149‑1 (MQ=255)
                                                                                                             gACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCCCGGCCAGGGTGGTCATCAGGa                                      >  6:101901/1‑149 (MQ=255)
                                                                                                               ccGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATc                                    <  4:120812/149‑1 (MQ=255)
                                                                                                                 ggAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCgg                                  <  6:343213/149‑1 (MQ=255)
                                                                                                                    gtggtgAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGtcatca                                         <  3:343941/139‑1 (MQ=255)
                                                                                                                    gtggtgAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGtcatca                                         >  4:343941/1‑139 (MQ=255)
                                                                                                                         gAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGggt                                                 >  7:133078/1‑126 (MQ=255)
                                                                                                                         gAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGggt                                                 <  8:133078/126‑1 (MQ=255)
                                                                                                                          aaCAGTGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAgcg                         >  6:243488/1‑149 (MQ=255)
                                                                                                                            cAGCGTCAGCATCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCgca                       <  1:53791/149‑1 (MQ=255)
                                                                                                                            cAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCCCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAgcg                         >  7:93143/1‑147 (MQ=255)
                                                                                                                             aGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCgcag                      <  5:148812/149‑1 (MQ=255)
                                                                                                                                 tCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAAc                  >  7:5946/1‑149 (MQ=255)
                                                                                                                                       ccTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACgccgcc            >  7:369405/1‑149 (MQ=255)
                                                                                                                                       ccTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACgccgcc            <  8:38022/149‑1 (MQ=255)
                                                                                                                                           gCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTg          <  5:243488/147‑1 (MQ=255)
                                                                                                                                           gCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGt        <  8:5946/149‑1 (MQ=255)
                                                                                                                                             tGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTgg         >  5:327418/1‑146 (MQ=255)
                                                                                                                                              gACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGtcatca                                         <  2:201348/113‑1 (MQ=255)
                                                                                                                                                 cAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATgg  >  2:220109/1‑149 (MQ=255)
                                                                                                                                                 cAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCACAGCAACGCCGCCTGGTAGATgg  >  6:280004/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCGCCACCGGGCGACGAATGAAAGGCCCGGACAGGTTCATGCCTGGGCCTGCTTCACGTCAGTGGCCGGGCGCAGGCGGCGGGCCAGGCGATCGAAGTACAGGTAGATGACCGGAGTGGTGAACAGCGTCAGCACCTGGCTGACCAGCAACCCGCCGACCATCACCAGGCCCAGTGGCTGGCGCAGCTCGGCGCCCGAACCGGTGGCAAGCATCAGCGGCACCGCACCAAACAGCGCGGCCAGGGTGGTCATCAGGATCGGCCGGAAGCGCAGCAACGCCGCCTGGTAGATGG  >  NC_002947/4069023‑4069316

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: