Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 298,120 C→G 5.3% intergenic (‑117/+37) 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,1200CG 5.3% 309.9 / 5.4 94intergenic (‑117/+37)gshA/PP_0244glutamate‑‑cysteine ligase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/49);  new base G (2/3);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

GGCTCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTGCA  >  NC_002947/297973‑298267
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                                             >  4:136002/1‑116 (MQ=255)
                                        cTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCGCTGGCActgct                                                                                                                                             <  3:136002/116‑1 (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                                                     >  2:257162/1‑149 (MQ=255)
                                                                                               gaagaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTc                                                     >  8:395463/1‑149 (MQ=255)
                                                                                                aagaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCg                                                    >  6:400403/1‑149 (MQ=255)
                                                                                                aagaaAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCg                                                    >  6:406792/1‑149 (MQ=255)
                                                                                                    aaGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                               >  8:26066/1‑134 (MQ=255)
                                                                                                    aaGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCt                                                               <  7:26066/134‑1 (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                                                                  >  3:409359/1‑120 (MQ=255)
                                                                                                               gTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGTCACCAGCTTGTCGc                                                                  <  4:409359/120‑1 (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       >  7:188677/1‑149 (MQ=255)
                                                                                                                                             tGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCg       <  4:334473/149‑1 (MQ=255)
                                                                                                                                              gCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCgg      >  5:407462/1‑149 (MQ=255)
                                                                                                                                               ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt     <  7:314953/149‑1 (MQ=255)
                                                                                                                                               ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt     <  8:258598/149‑1 (MQ=255)
                                                                                                                                               ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt     >  6:140212/1‑149 (MQ=255)
                                                                                                                                               ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt     <  6:108467/149‑1 (MQ=255)
                                                                                                                                               ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt     <  4:238779/149‑1 (MQ=255)
                                                                                                                                               ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt     >  8:378529/1‑149 (MQ=255)
                                                                                                                                               ccAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGt     >  1:123714/1‑149 (MQ=255)
                                                                                                                                                cAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTg    <  6:407462/149‑1 (MQ=255)
                                                                                                                                                 aGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTGc   <  2:8326/149‑1 (MQ=255)
                                                                                                                                                   cgcACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGCGCGGTGca  >  7:322700/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCTCAGGCGGCGGTTGAGGAGTTCGCTCAAGATGTATTCCTTCACGCGTCAGTCGCCCCAATATGGGGGTAGAGATAACGGTCTACAAGGGTAGGAAGAAAGGAACTGGCGTTGTCGCCTGGTTTACGCGGTTCCAGCAGTGCCAGCGCACTGCTGAAAATACCGCAGATACTGCTTGGTGAGCTAGAGAACCGCGAAGGTTCCTTGCGCTTTGGCCACCAGCTTGTCGCCCTGAACTACGTCGGCGTCGACCACCAGGGTGCGCCGCCCGGCGTGCAGCACACGAGCGGTGCA  >  NC_002947/297973‑298267

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: