Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,640,827 C→G 6.0% V326L (GTG→CTG)  gabP‑IV ← gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,640,8270CG 6.0% 219.5 / 3.8 67V326L (GTG→CTG) gabP‑IVgamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base G (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 = 8.86e-01

ATGGCGCCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACG  >  NC_002947/4640688‑4640974
                                                                                                                                           |                                                                                                                                                   
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                                                                                                                  <  5:271350/146‑1 (MQ=255)
                             aaaCACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATaga                                                                                                                  >  6:271350/1‑146 (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                                                                                                                                          >  1:122071/1‑101 (MQ=255)
                                                  aTAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCt                                                                                                                                          <  2:122071/101‑1 (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                                                                                 >  5:158712/1‑149 (MQ=255)
                                                           gAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCa                                                                                 <  8:138863/149‑1 (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                                                                                        >  5:271260/1‑135 (MQ=255)
                                                                  ccagcGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGcccc                                                                                        <  6:271260/135‑1 (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:30752/84‑1 (MQ=255)
                                                                                                                              aTGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATg                                                                               <  8:25268/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)
                                                                                                                                           |                                                                                                                                                   
ATGGCGCCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACG  >  NC_002947/4640688‑4640974

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