Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,640,830 T→G 5.9% S325R (AGC→CGC)  gabP‑IV ← gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,640,8300TG 5.9% 215.9 / 3.8 68S325R (AGC→CGC) gabP‑IVgamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base G (2/2);  total (35/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 = 8.89e-01

ATGGCGCCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTG  >  NC_002947/4640688‑4640977
                                                                                                                                              |                                                                                                                                                   
aTGGCGCCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAgcgc                                                                                                                                               <  7:36857/149‑1 (MQ=255)
   gcgcCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCt                                                                                                                                             <  5:121044/148‑1 (MQ=255)
     gcCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGAc                                                                                                                                          >  4:63778/1‑149 (MQ=255)
          tggtCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATTCGCTTGGTCACGCTCAGCGCCTGACCAcgc                                                                                                                                     <  1:243150/149‑1 (MQ=255)
                 cagcagGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGcc                                                                                                                                >  4:63562/1‑147 (MQ=255)
                    cagGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGc                                                                                                                            >  6:89237/1‑148 (MQ=255)
                     agGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCg                                                                                                                           >  7:138863/1‑148 (MQ=255)
                            cAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGc                                                                                                                   <  3:63778/149‑1 (MQ=255)
                             aaaCACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATaga                                                                                                                     >  6:271350/1‑146 (MQ=255)
                             aaaCACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATaga                                                                                                                     <  5:271350/146‑1 (MQ=255)
                             aaaCACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGc                                                                                                                   <  7:278709/148‑1 (MQ=255)
                             aaaCACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGc                                                                                                                   >  8:278709/1‑148 (MQ=255)
                             aaaCACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCGTGGTCAGGCGCTGAGCGTGACCAagcg                                                                                                                                    <  5:70698/131‑1 (MQ=38)
                             aaaCACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCGTGGTCAGGCGCTGAGCGTGACCAagcg                                                                                                                                    >  6:70698/1‑131 (MQ=38)
                                   tttGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGCGCATGcgcg                                                                                                            >  8:164249/1‑149 (MQ=255)
                                        ccGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGcc                                                                                                       >  3:55304/1‑149 (MQ=255)
                                         cGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCg                                                                                                      >  6:325494/1‑149 (MQ=255)
                                                 caTAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGc                                                                                              >  8:283461/1‑149 (MQ=255)
                                                  aTAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCt                                                                                                                                             <  2:122071/101‑1 (MQ=255)
                                                  aTAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCt                                                                                                                                             >  1:122071/1‑101 (MQ=255)
                                                  aTAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGcc                                                                                             >  4:2010/1‑149 (MQ=255)
                                                     gTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATgc                                                                                                               >  2:261448/1‑128 (MQ=255)
                                                     gTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATgc                                                                                                               <  1:261448/128‑1 (MQ=255)
                                                     gTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCg                                                                                          <  5:89237/149‑1 (MQ=255)
                                                     gTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAAGAGAGCATGCGCGACGCCGTGTAGAGCCCCg                                                                                          >  8:243705/1‑149 (MQ=255)
                                                      ttGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGa                                                                                         <  1:160420/149‑1 (MQ=255)
                                                         aCGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGtt                                                                                      <  3:2010/149‑1 (MQ=255)
                                                          cGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTc                                                                                     <  4:55304/149‑1 (MQ=255)
                                                           gAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCa                                                                                    <  8:138863/149‑1 (MQ=255)
                                                           gAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCa                                                                                    >  5:158712/1‑149 (MQ=255)
                                                            aaGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCAt                                                                                   <  3:63562/149‑1 (MQ=255)
                                                            aaGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGCGTAGAGCCCCGAGTTCAt                                                                                   <  2:31302/149‑1 (MQ=255)
                                                                 gccagcGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGcc                                                                                             <  2:354377/134‑1 (MQ=255)
                                                                 gccagcGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGcc                                                                                             >  1:354377/1‑134 (MQ=255)
                                                                  ccagcGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGcccc                                                                                           <  6:271260/135‑1 (MQ=255)
                                                                  ccagcGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGcccc                                                                                           >  5:271260/1‑135 (MQ=255)
                                                                   cagcGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACtt                                                                            <  7:283461/149‑1 (MQ=255)
                                                                            ccGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCgt                                                                                                     >  2:255701/1‑115 (MQ=255)
                                                                            ccGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGc                                                                                              >  4:262448/1‑122 (MQ=255)
                                                                            ccGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGc                                                                                              <  3:262448/122‑1 (MQ=255)
                                                                            ccGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCACCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCgt                                                                                                     <  1:255701/115‑1 (MQ=255)
                                                                             cGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGc                                                                                                                            <  7:6128/91‑1 (MQ=255)
                                                                             cGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGc                                                                                                                            >  8:6128/1‑91 (MQ=255)
                                                                                   ggTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCACCGAGTTCATGCAACTTGTCAACGCCACCATca                                                            >  4:188458/1‑149 (MQ=255)
                                                                                   ggTGGAGAAGATCACTGCCACGATCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGca                                                            >  7:224393/1‑149 (MQ=255)
                                                                                      ggAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGcacca                                                         >  4:166604/1‑149 (MQ=255)
                                                                                          aagaTCACTGCCACGGTCGGCACGCCGGAACCGGAGAAGCGCTTCGTCACGGTCAGCGCCTGACCACGCGCGACCAGCGAATAGAGCATGCGAGACCCCGTGTAGAGCACCCAGGACAAGCAACAAGTCCGctccaccagctccaccag                                                     >  4:279576/1‑149 (MQ=255)
                                                                                                aCTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAAc                                               <  6:158712/149‑1 (MQ=255)
                                                                                                  tGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAAc                                               >  5:141713/1‑147 (MQ=255)
                                                                                                    ccACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAg                                            >  5:334876/1‑148 (MQ=255)
                                                                                                          taggCACGCCTGAACCGGAGATGCTCTTGGTCACGCTAATCGCCTTACCACGCGCGCCAAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCTAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTg                                     <  3:224848/147‑1 (MQ=255)
                                                                                                          tCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTg                                     <  5:325494/149‑1 (MQ=255)
                                                                                                                 gCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGCGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCg                              >  2:206037/1‑149 (MQ=255)
                                                                                                                   cGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGc                            >  5:321044/1‑149 (MQ=255)
                                                                                                                       aCCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGt                        >  7:294609/1‑149 (MQ=255)
                                                                                                                          ggAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCa                     <  3:166604/149‑1 (MQ=255)
                                                                                                                           gagaTGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAg                    <  8:294609/149‑1 (MQ=255)
                                                                                                                              aTGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATg                                                                                  >  7:30752/1‑84 (MQ=255)
                                                                                                                              aTGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATg                                                                                  >  7:25268/1‑84 (MQ=255)
                                                                                                                              aTGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATg                                                                                  <  8:25268/84‑1 (MQ=255)
                                                                                                                              aTGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATg                                                                                  <  8:30752/84‑1 (MQ=255)
                                                                                                                                 cgcTTGGTCACGCTCATCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTc              <  6:141713/149‑1 (MQ=255)
                                                                                                                                 cgcTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTc              <  6:334876/149‑1 (MQ=255)
                                                                                                                                  gcgtGGTCAGGCGCTGAGCGTGCCCAAGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCAc                          >  5:104583/4‑136 (MQ=255)
                                                                                                                                  gcgtGGTCAGGCGCTGAGCGTGACCAAGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCAc                          <  6:104583/133‑1 (MQ=255)
                                                                                                                                       gTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCCACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTa         >  1:129755/1‑148 (MQ=255)
                                                                                                                                          aCGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACg     >  5:327069/1‑149 (MQ=255)
                                                                                                                                             cTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCTAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTg  <  8:224393/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
ATGGCGCCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTG  >  NC_002947/4640688‑4640977

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