Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,613,039 G→A 6.6% G215G (GGC→GGT)  mcoA ← Mn(II) copper oxidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,613,0390GA 6.6% 178.0 / 2.1 60G215G (GGC→GGT) mcoAMn(II) copper oxidase A
Reads supporting (aligned to +/- strand):  ref base G (24/32);  new base A (2/2);  total (26/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.86e-01

CACCGACGGCATAGTTGTGCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGCTTCCAAAGCGCTGCCACT  >  NC_002947/3612900‑3613185
                                                                                                                                           |                                                                                                                                                  
cACCGACGGCATAGTTGTGCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCg                                                                                                                                           >  3:38387/1‑149 (MQ=255)
 aCCGACGGCATAGTTGTGCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCgg                                                                                                                                          <  1:141757/149‑1 (MQ=255)
 aCCGACGGCATAGTTGTGCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCAAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCgg                                                                                                                                          >  7:166430/1‑149 (MQ=255)
   cGACGGCATAGTTGTGCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGc                                                                                                                                        <  4:307022/149‑1 (MQ=255)
     aCGGCATAGTTGTGCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGTGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGa                                                                                                                                      <  2:96564/149‑1 (MQ=255)
              ttGTGCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCgg                                                                                                                             <  4:38387/149‑1 (MQ=255)
                gtgCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCt                                                                                                                           >  5:235450/1‑149 (MQ=255)
                 tgCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTg                                                                                                                          >  3:326339/1‑149 (MQ=255)
                 tgCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTg                                                                                                                          <  2:279934/149‑1 (MQ=255)
                            gcGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAgg                                                                                                                >  1:229805/1‑148 (MQ=255)
                            gcGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAgg                                                                                                                <  2:229805/148‑1 (MQ=255)
                             cGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGGTCCAGCGCTGGGGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGaa                                                                                                              <  5:284326/149‑1 (MQ=255)
                              ggTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAAc                                                                                                             <  1:227727/149‑1 (MQ=255)
                                      aGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCTTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGAt                                                                                                     <  4:195390/149‑1 (MQ=255)
                                         ccGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTg                                                                                                  <  1:319397/149‑1 (MQ=255)
                                          cGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTgc                                                                                                 <  6:2918/149‑1 (MQ=255)
                                           gAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTgcg                                                                                                >  2:141182/1‑149 (MQ=255)
                                             ggTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGcc                                                                                              <  7:5798/149‑1 (MQ=255)
                                                          aCCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGc                                                                                 >  2:65605/1‑149 (MQ=255)
                                                              gCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTc                                                                             <  1:65605/149‑1 (MQ=255)
                                                              cccTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGACCATCTCGATCTCCGCCTTCCATGGGTTGCGGTc                                                                             >  8:256900/2‑149 (MQ=255)
                                                                          ttGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCAt                                                                                        <  6:196569/126‑1 (MQ=255)
                                                                          ttGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCAt                                                                                        >  5:196569/1‑126 (MQ=255)
                                                                          ttGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTgg                                                                  >  2:323109/1‑148 (MQ=255)
                                                                            gAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTgg                                                                  <  1:323109/146‑1 (MQ=255)
                                                                            gAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCg                                                                >  4:280797/1‑148 (MQ=255)
                                                                             aaCGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTa                                                              >  8:143375/1‑149 (MQ=255)
                                                                              aCGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTAc                                                             >  5:132115/1‑149 (MQ=255)
                                                                                gCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTg                                                           <  6:12067/149‑1 (MQ=255)
                                                                                 ctgcctgcGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTgg                                                          <  7:143375/149‑1 (MQ=255)
                                                                                  tgcctgcGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGa                                                                                                                  <  2:123655/92‑1 (MQ=255)
                                                                                  tgcctgcGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTggg                                                         >  5:308051/1‑149 (MQ=255)
                                                                                     ctgcGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCg                                                      <  1:141182/149‑1 (MQ=255)
                                                                                     ctgcGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCg                                                      >  7:171353/1‑149 (MQ=255)
                                                                                        ccggTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCggg                                                    <  8:132382/146‑1 (MQ=255)
                                                                                         gggTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGaa                                                  >  5:115634/1‑149 (MQ=255)
                                                                                          ggTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGACGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAAc                                                 >  7:299351/1‑149 (MQ=255)
                                                                                           gTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACg                                                <  4:326339/149‑1 (MQ=255)
                                                                                           gTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACg                                                <  6:235450/149‑1 (MQ=255)
                                                                                              gAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGTGTCGGGAACGGGt                                             <  6:308051/149‑1 (MQ=255)
                                                                                               aaCTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGACCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCg                                                                                <  5:43728/113‑1 (MQ=255)
                                                                                               aaCTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGACCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCg                                                                                <  5:49249/113‑1 (MQ=255)
                                                                                               aaCTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGACCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCg                                                                                >  6:49249/1‑113 (MQ=255)
                                                                                               aaCTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGACCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCg                                                                                >  6:43728/1‑113 (MQ=255)
                                                                                                aCTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTAc                                           <  3:280797/149‑1 (MQ=255)
                                                                                                          gcgcTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGtt                                 <  6:132115/149‑1 (MQ=255)
                                                                                                          gcgcTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGtt                                 >  5:273631/1‑149 (MQ=255)
                                                                                                             cTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCt                              <  8:171353/149‑1 (MQ=255)
                                                                                                                gTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGTGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGa                           <  7:256900/149‑1 (MQ=255)
                                                                                                                gTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCggg                                   <  6:336226/141‑1 (MQ=255)
                                                                                                                gTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCggg                                   >  5:336226/1‑141 (MQ=255)
                                                                                                                gTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTg                            <  1:62668/148‑1 (MQ=255)
                                                                                                                gTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGa                           <  6:115634/149‑1 (MQ=255)
                                                                                                                gTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGa                           <  2:60899/149‑1 (MQ=255)
                                                                                                                      acAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGCCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGAGTACTGGGTCGGGAACGGGTACAGGCCGGATTGCTTGAGGAACg                     >  2:106663/2‑149 (MQ=255)
                                                                                                                       cAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGc                    >  4:152310/1‑149 (MQ=255)
                                                                                                                        aGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGTGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGCt                   >  1:231375/1‑149 (MQ=255)
                                                                                                                         gCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGCtt                  >  6:51603/1‑149 (MQ=255)
                                                                                                                                  tggtggCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGCTTCCAAAgcgc         >  3:191497/1‑149 (MQ=255)
                                                                                                                                    gtggCCGGCCTTCAGCGGGCGAGCCGACCGGCTGACGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGCTTCCAAAGCGCTg       <  8:299351/149‑1 (MQ=255)
                                                                                                                                         cggccTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGCTTCCAAAGCGCTGCCACt  >  2:96471/1‑149 (MQ=255)
                                                                                                                                         cggccTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGCTTCCAAAGCGCTGCCACt  <  8:166430/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CACCGACGGCATAGTTGTGCAGTTGCTTGCGGTCGCGCAGGCCGAGGTTGATGCGCGCACCGGCCTGGGCGGTCTTGAACGCTGCCTGCGGGTAGAACTCGTTCCAGCGCTGGTGCGACCAGCCTTTTCCTGGTGGCCGGCCTTCAGCGGGCGAGCCGACCGGCTGGCGGTTGAGGAACATCTCGATCTGCGCCTTCCATGGGTTGCGGTCCAGCACGTTGGCGTACTGGGTCGGGAACGGGTACAGGCCGGGTTGCTTGAGGAACGCTTCCAAAGCGCTGCCACT  >  NC_002947/3612900‑3613185

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