Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,640,837 C→G 5.6% Q322H (CAG→CAC)  gabP‑IV ← gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,640,8370CG 5.6% 236.0 / 3.8 72Q322H (CAG→CAC) gabP‑IVgamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base C (36/32);  new base G (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

GCGCCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTGGAACGAC  >  NC_002947/4640691‑4640984
                                                                                                                                                  |                                                                                                                                                   
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                                                                                                                                                    >  1:122071/1‑101 (MQ=255)
                                               aTAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCt                                                                                                                                                    <  2:122071/101‑1 (MQ=255)
                                               aTAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGcc                                                                                                    >  4:2010/1‑149 (MQ=255)
                                                  gTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATgc                                                                                                                      <  1:261448/128‑1 (MQ=255)
                                                  gTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATgc                                                                                                                      >  2:261448/1‑128 (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                                                                                                     <  3:262448/122‑1 (MQ=255)
                                                                         ccGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGc                                                                                                     >  4:262448/1‑122 (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                                                                                         <  8:30752/84‑1 (MQ=255)
                                                                                                                           aTGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATg                                                                                         >  7:25268/1‑84 (MQ=255)
                                                                                                                           aTGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATg                                                                                         >  7:30752/1‑84 (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)
                                                                                                                                          cTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCTAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTg         <  8:224393/149‑1 (MQ=255)
                                                                                                                                             aGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTg         <  2:353752/146‑1 (MQ=255)
                                                                                                                                             aGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTg         >  1:353752/1‑146 (MQ=255)
                                                                                                                                              gcgcCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAg                  <  8:195266/136‑1 (MQ=255)
                                                                                                                                              gcgcCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAg                  >  7:195266/1‑136 (MQ=255)
                                                                                                                                                gcCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTGGAACGa   >  5:246576/1‑149 (MQ=255)
                                                                                                                                                 ccTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTGGAACGAc  >  2:313216/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGCCGGTGGTCGACAGCAGGAAGCCAAACACTTTGCCCGGGAACACATAGTTGACGAAGCAGCCAGCGAAGCCCGCCAGGGTGGAGAAGATCACTGCCACGGTCGGCACGCCGGAACCGGAGATGCGCTTGGTCACGCTCAGCGCCTGACCACGCGCGCCCAGCGAATAGAGCATGCGCGACGCCGTGTAGAGCCCCGAGTTCATGCAACTTGTCACCGCCACCAGCACCACCAGGTCAACCAGCAACTTGGCGCCCGGCACGTTCAGCACTTCCAGTACACGCTGGAACGAC  >  NC_002947/4640691‑4640984

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