Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 298,117 C→G 5.5% intergenic (‑114/+40) gshA ← / ← PP_0244 glutamate‑‑cysteine ligase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947298,1170CG 5.5% 281.5 / 3.3 87intergenic (‑114/+40)gshA/PP_0244glutamate‑‑cysteine ligase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/44);  new base G (2/3);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCTCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTG  >  NC_002947/297971‑298265
                                                                                                                                                  |                                                                                                                                                    
cagacaaaggaggCGGATGAGGAGTTAGCTCAAGATCTATTCCTTCACGCGTCAGTCGCCCCAATATCGGGGTAGAGATAACGGTCTACAAGGGTATGAAGAAAGGAACTGGCGTTGTCGCCAGGTTTACGCGGTTCCAGCAGTGCCAg                                                                                                                                                    >  6:424746/8‑149 (MQ=255)
  ggCTCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAgcg                                                                                                                                                  <  4:207303/149‑1 (MQ=255)
         gcggcgGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCGCTGGCActgct                                                                                                                                           <  8:289163/149‑1 (MQ=255)
              ggTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGaaaa                                                                                                                                      >  5:317816/1‑149 (MQ=255)
               gTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCGTCGTTTACGCGGTTCCAGCACTGCCAGCGCACTGCTGaaaa                                                                                                                                      >  6:24522/1‑148 (MQ=255)
               gTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCGCTGGCACTGCTGgaacc                                                                                                                                     <  1:31836/149‑3 (MQ=255)
                  gaggagTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATAc                                                                                                                                   <  6:71895/148‑1 (MQ=255)
                       gTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCg                                                                                                                                 >  7:402535/1‑145 (MQ=255)
                       gTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGa                                                                                                                             <  6:403143/149‑1 (MQ=255)
                          cGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATa                                                                                                                           <  5:133289/148‑1 (MQ=255)
                           gCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATAc                                                                                                                          <  1:380346/148‑1 (MQ=255)
                           ccgaaaGATGTATTCCTTCACGCGGCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTATACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATAc                                                                                                                          <  3:269538/144‑1 (MQ=255)
                             tCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGc                                                                                                                       >  2:6892/1‑149 (MQ=255)
                                aGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGTGTAGAGATAACGGTCTACAATGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    <  5:46704/149‑1 (MQ=255)
                                aGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    >  5:626/1‑149 (MQ=255)
                                  aTGTATTCCTTCACGCGTCAGTCGCCGCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGTAGAAAGGAACTGTCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGt                                                                                                                  <  1:169535/149‑1 (MQ=255)
                                   tGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTg                                                                                                                 <  6:317816/149‑1 (MQ=255)
                                    gTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGt                                                                                                                  >  6:200331/1‑147 (MQ=255)
                                       ttCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGc                                                                                                              <  3:221521/148‑1 (MQ=255)
                                          cTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCGCTGGCActgct                                                                                                                                           <  3:136002/116‑1 (MQ=255)
                                          cTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCGCTGGCActgct                                                                                                                                           >  4:136002/1‑116 (MQ=255)
                                          cTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTag                                                                                                           <  6:626/148‑1 (MQ=255)
                                          cTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTa                                                                                                            >  8:233322/1‑147 (MQ=255)
                                              aCGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAAc                                                                                                      >  2:409606/1‑149 (MQ=255)
                                                   tCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGa                                                                                                 <  2:179526/149‑1 (MQ=255)
                                                     aGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGaa                                                                                                <  8:231639/148‑1 (MQ=255)
                                                     aGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAg                                                                                               >  7:237125/1‑149 (MQ=255)
                                                      gTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTCCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAgg                                                                                              >  4:312156/1‑149 (MQ=255)
                                                         gCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGtt                                                                                            <  3:34948/148‑1 (MQ=255)
                                                         gCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTCCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTc                                                                                           <  3:312156/149‑1 (MQ=255)
                                                          ccccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    >  1:180553/1‑123 (MQ=255)
                                                          ccccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTg                                                                                                                    <  2:180553/123‑1 (MQ=255)
                                                          ccccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTcc                                                                                          >  3:227501/1‑149 (MQ=255)
                                                            ccAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                        >  3:84658/1‑149 (MQ=255)
                                                            ccAATATGGGGGTAGAGATAACGCTCTACAAGGGTAGGAAGAATGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCAGTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                        <  5:24522/149‑1 (MQ=255)
                                                             cAATATGGGGTTAGAGATAACGGGCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTg                                                                                       <  2:261841/149‑1 (MQ=255)
                                                                    ggggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCttt                                                                                 <  5:105952/148‑1 (MQ=255)
                                                                     gggTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTgg                                                                               >  7:328765/1‑149 (MQ=255)
                                                                          gagaTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGccac                                                                           <  1:6892/148‑1 (MQ=255)
                                                                          gagaTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGccac                                                                           <  8:402535/148‑1 (MQ=255)
                                                                            gaTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAg                                                                        >  7:207691/1‑149 (MQ=255)
                                                                             aTAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGc                                                                       <  4:161397/149‑1 (MQ=255)
                                                                                   gTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                        >  1:161769/1‑126 (MQ=255)
                                                                                   gTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCtt                                                                                        <  2:161769/126‑1 (MQ=255)
                                                                                        cAAGGGTAGGTAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGccc                                                              <  7:202221/147‑1 (MQ=255)
                                                                                        cAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGccc                                                              <  8:237125/147‑1 (MQ=255)
                                                                                        cAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGcc                                                               <  4:84658/146‑1 (MQ=255)
                                                                                        cAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGcc                                                               <  2:49380/146‑1 (MQ=255)
                                                                                            ggTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAAc                                                         >  7:258598/1‑148 (MQ=255)
                                                                                              tAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTg                                                                                                                        <  5:350235/83‑1 (MQ=255)
                                                                                              tAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTg                                                                                                                        >  6:350235/1‑83 (MQ=255)
                                                                                              tAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTAc                                                      <  7:254866/149‑1 (MQ=255)
                                                                                               aGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACg                                                     >  5:108467/1‑149 (MQ=255)
                                                                                               aGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACg                                                     >  6:203964/1‑149 (MQ=255)
                                                                                                 gaagaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                   >  8:395463/1‑149 (MQ=255)
                                                                                                 gaagaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                   >  2:257162/1‑149 (MQ=255)
                                                                                                  aagaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCg                                                  >  6:400403/1‑149 (MQ=255)
                                                                                                  aagaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCg                                                  >  6:406792/1‑149 (MQ=255)
                                                                                                      aaGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                             <  7:26066/134‑1 (MQ=255)
                                                                                                      aaGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                             >  8:26066/1‑134 (MQ=255)
                                                                                                       aGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTc                                             <  4:227501/149‑1 (MQ=255)
                                                                                                        ggAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTc                                             >  7:284063/1‑148 (MQ=255)
                                                                                                        ggAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCg                                            <  5:400403/149‑1 (MQ=255)
                                                                                                        ggAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCg                                            <  5:406792/149‑1 (MQ=255)
                                                                                                             tGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGaccac                                       >  3:404392/1‑149 (MQ=255)
                                                                                                             tGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGaccac                                       >  3:238779/1‑149 (MQ=255)
                                                                                                               gCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGaccac                                       >  3:334473/1‑147 (MQ=255)
                                                                                                                 gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGTCACCAGCTTGTCGc                                                                <  4:409359/120‑1 (MQ=255)
                                                                                                                 gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGTCACCAGCTTGTCGc                                                                >  3:409359/1‑120 (MQ=255)
                                                                                                                 gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAgg                                   <  8:328765/149‑1 (MQ=255)
                                                                                                                       gCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTgcgc                             <  1:409606/149‑1 (MQ=255)
                                                                                                                        ccTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                   <  8:329195/126‑1 (MQ=255)
                                                                                                                        ccTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                   >  7:329195/1‑126 (MQ=255)
                                                                                                                         cTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGcgccg                           >  8:314953/1‑149 (MQ=255)
                                                                                                                            gTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGcgccgc                          <  7:233322/147‑1 (MQ=255)
                                                                                                                             tttACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCg                       >  1:8326/1‑149 (MQ=255)
                                                                                                                                 cgctGTTCCAGCAGTGCCAGCGCACTTCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGt                   <  5:424746/149‑1 (MQ=255)
                                                                                                                                  gcgGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTg                  <  8:207691/149‑1 (MQ=255)
                                                                                                                                   cgGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgc                 <  4:404392/149‑1 (MQ=255)
                                                                                                                                    ggTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgca                <  6:270895/149‑1 (MQ=255)
                                                                                                                                        tcagcagTGCGCTGGCACTGCTGGAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgcagca             >  1:54280/2‑148 (MQ=255)
                                                                                                                                        ccAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTgcag               <  1:437811/146‑1 (MQ=255)
                                                                                                                                         cagcagTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGcaca           >  5:11819/1‑149 (MQ=255)
                                                                                                                                           gcagTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCacac          >  7:387794/1‑148 (MQ=255)
                                                                                                                                               tGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCg     <  4:334473/149‑1 (MQ=255)
                                                                                                                                               tGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCg     >  7:188677/1‑149 (MQ=255)
                                                                                                                                                gCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCgg    >  5:407462/1‑149 (MQ=255)
                                                                                                                                                 ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt   <  7:314953/149‑1 (MQ=255)
                                                                                                                                                 ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt   <  8:258598/149‑1 (MQ=255)
                                                                                                                                                 ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt   <  4:238779/149‑1 (MQ=255)
                                                                                                                                                 ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt   <  6:108467/149‑1 (MQ=255)
                                                                                                                                                 ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt   >  8:378529/1‑149 (MQ=255)
                                                                                                                                                 ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt   >  6:140212/1‑149 (MQ=255)
                                                                                                                                                 ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt   >  1:123714/1‑149 (MQ=255)
                                                                                                                                                  cAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTg  <  6:407462/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAGGCTCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTG  >  NC_002947/297971‑298265

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: