Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 484,692 G→C 7.1% intergenic (‑254/+29) yeaG ← / ← glpE protein kinase/thiosulfate sulfurtransferase GlpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947484,6920GC 7.1% 261.9 / 9.1 85intergenic (‑254/+29)yeaG/glpEprotein kinase/thiosulfate sulfurtransferase GlpE
Reads supporting (aligned to +/- strand):  ref base G (36/43);  new base C (3/3);  total (39/46)
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.32e-01

TACCGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCA  >  NC_002947/484547‑484835
                                                                                                                                                 |                                                                                                                                                
tACCGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCgcggcg                                                                                                                                                >  1:344787/1‑148 (MQ=255)
  ccGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGgcgc                                                                                                                                               <  8:68060/147‑1 (MQ=255)
   cGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACg                                                                                                                                            <  8:327779/149‑1 (MQ=255)
     aTGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACggg                                                                                                                                          <  2:233905/149‑1 (MQ=255)
      tGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGc                                                                                                                                         <  7:284448/149‑1 (MQ=255)
          ccGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGCACGCCCGGTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTg                                                                                                                                       >  4:252694/1‑147 (MQ=255)
          ccGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACAGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTg                                                                                                                                       >  1:297605/1‑147 (MQ=255)
          ccGCACGCCGGTGTACCGGCTCTCCACTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCa                                                                                                                                     <  6:305023/149‑1 (MQ=255)
           cGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAg                                                                                                                                    <  8:162117/149‑1 (MQ=255)
                gCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttg                                                                                                                                      >  7:351671/1‑140 (MQ=255)
                gCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttg                                                                                                                                      <  8:351671/142‑3 (MQ=255)
                 ccGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGaaa                                                                                                                              <  8:295787/149‑1 (MQ=255)
                    atgtACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAAt                                                                                                                            >  3:76738/2‑148 (MQ=255)
                       tACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATaaa                                                                                                                         <  4:76738/148‑1 (MQ=255)
                         ccGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATaaaaaa                                                                                                                      <  4:270564/149‑1 (MQ=255)
                           ggCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAatt                                                                                                                    <  1:135510/149‑1 (MQ=255)
                                tcCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAattatt                                                                                                                 >  8:345660/1‑147 (MQ=255)
                                tcCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAattatt                                                                                                                 >  2:39956/1‑147 (MQ=255)
                                tcCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAattatt                                                                                                                 <  2:234981/147‑1 (MQ=255)
                                  gccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCATGGTAAAGGGGGGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGc                                                                                                            >  4:212349/2‑149 (MQ=255)
                                  cccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGTGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGc                                                                                                            <  4:16215/149‑1 (MQ=255)
                                  cccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCa                                                                                                               <  8:190542/147‑1 (MQ=255)
                                  cccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGc                                                                                                            >  3:382924/1‑149 (MQ=255)
                                   ccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCg                                                                                                           >  2:299405/1‑149 (MQ=255)
                                    cTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGt                                                                                                          >  7:278581/1‑149 (MQ=255)
                                          ggATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGccc                                                                                                     <  2:386737/148‑1 (MQ=255)
                                            aTGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTgg                                                                                                  <  8:37818/149‑1 (MQ=255)
                                               gCCTGAGCTTAAGCGTAGTTGGTTTTCTGGAAGGTAAAGGGCTGAGAAAGGCTGAAGTTACGAACGACCGTGCGTCAGATAAGCCGCAAGGCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCg                                                                                               <  7:303311/149‑1 (MQ=255)
                                                ccTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGc                                                                                                                                      <  6:323168/110‑1 (MQ=255)
                                                ccTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGc                                                                                                                                      >  5:323168/1‑110 (MQ=255)
                                                ccTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCgg                                                                                              >  4:64218/1‑149 (MQ=255)
                                                ccTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCgc                                                                                              >  3:382177/1‑148 (MQ=255)
                                                 cTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTGCCCTGGGCGGtc                                                                                            >  3:380484/1‑149 (MQ=255)
                                                   gAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAatta                                                                                                                   <  3:322826/126‑1 (MQ=255)
                                                   gAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAatta                                                                                                                   >  4:322826/1‑126 (MQ=255)
                                                      cTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCgg                                                                                        >  6:213003/1‑149 (MQ=255)
                                                       ttAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCAGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGAGTTGCCCTGTGCGGTCTCGGc                                                                                       >  5:200866/1‑149 (MQ=255)
                                                       ttAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGACAAACCAAGAAGTTACGAACGACAGTTAGTCAGATAAGCAGCAAGCCAGCGCGGAGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGAGGTCTAGGc                                                                                       >  7:295951/1‑149 (MQ=255)
                                                                gttgttTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCGTGGGCGGTGACGGCCGGGGAGGt                                                                              >  4:250996/1‑149 (MQ=255)
                                                                gTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGt                                                                              >  6:310776/1‑149 (MQ=255)
                                                                 ttGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTc                                                                             >  6:360018/1‑149 (MQ=255)
                                                                    gTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTcgcg                                                                          >  4:78007/1‑149 (MQ=255)
                                                                     ttATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTcgcgc                                                                         <  3:64218/149‑1 (MQ=255)
                                                                          tGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttgcg                                                                                                                                    <  5:154011/86‑5 (MQ=38)
                                                                          tGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttgcg                                                                                                                                    >  6:154011/1‑82 (MQ=38)
                                                                          tGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCca                                                                     >  4:331010/1‑148 (MQ=255)
                                                                            cAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCca                                                                     <  4:72330/146‑1 (MQ=255)
                                                                            cAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAg                                                                  <  5:213003/149‑1 (MQ=255)
                                                                            cAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCcac                                                                    <  4:380484/146‑1 (MQ=255)
                                                                               ggTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTc                                                               <  2:83713/149‑1 (MQ=255)
                                                                                 tAAAGGGGTGAGAAAGGCAGAAGTTACGAACGCCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGCAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCg                                                              >  7:338339/1‑148 (MQ=255)
                                                                                   aaGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                            <  4:382924/148‑1 (MQ=255)
                                                                                     ggggTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCACGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGc                                                          <  1:299405/149‑1 (MQ=255)
                                                                                      ggggGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCt                                                        <  5:89642/149‑1 (MQ=255)
                                                                                      gggTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCt                                                        >  2:261995/1‑149 (MQ=255)
                                                                                              aaGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCAc                                                <  8:278581/149‑1 (MQ=255)
                                                                                              aaGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCAc                                                <  8:307415/149‑1 (MQ=255)
                                                                                                  cAGAAGTTACGAACGACCGTTCGTGAGATATGCCGCAAGCCCGCGCTGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGGGTTGGCCTTGGCGGTCTAGTCCGGGTAGGTAGCGCGCCACAGTTCGAAGCCTCCGTCCTCGCt                                             <  8:295951/148‑1 (MQ=255)
                                                                                                   aGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTg                                            >  3:157350/1‑148 (MQ=255)
                                                                                                          acgaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTc                                    >  6:197306/1‑149 (MQ=255)
                                                                                                          acgaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTAAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTc                                    >  3:122319/1‑149 (MQ=255)
                                                                                                           cgaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCg                                   <  1:39956/149‑1 (MQ=255)
                                                                                                            gaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTa                                          <  1:66468/141‑1 (MQ=255)
                                                                                                            gaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTa                                          >  2:66468/1‑141 (MQ=255)
                                                                                                            gaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCgg                                  <  8:290436/149‑1 (MQ=255)
                                                                                                             aacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCGCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGga                                 <  5:310776/149‑1 (MQ=255)
                                                                                                                  ccGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAAc                            >  3:374596/1‑149 (MQ=255)
                                                                                                                   cGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAAcc                           <  5:360018/149‑1 (MQ=255)
                                                                                                                        gTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTgg                       <  1:189442/148‑1 (MQ=255)
                                                                                                                        gTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGc                      >  7:10913/1‑149 (MQ=255)
                                                                                                                          cAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGcc                     >  1:139830/1‑148 (MQ=255)
                                                                                                                              tAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCa                <  4:157350/149‑1 (MQ=255)
                                                                                                                              tAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCa                <  3:252694/149‑1 (MQ=255)
                                                                                                                               aaGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCTACAGTTCGAAGCCTCCGTCCACGCTGTACACGt                                     <  1:53935/127‑1 (MQ=255)
                                                                                                                               aaGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGt                                     >  2:53935/1‑127 (MQ=255)
                                                                                                                                aGCCGCACGCCCGCGCGGCTCACGGGCTGCAGGTGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCTgg              <  8:338339/149‑1 (MQ=255)
                                                                                                                                aGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAg               <  5:347600/148‑1 (MQ=255)
                                                                                                                                 gCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGt             <  3:212349/149‑1 (MQ=255)
                                                                                                                                  ccGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAgg              >  2:392676/1‑147 (MQ=255)
                                                                                                                                   ctgcaGCCCGTGCGCCGCGCGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAg           >  8:217767/5‑149 (MQ=255)
                                                                                                                                    gCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACAAGGTAgg          >  5:111471/1‑149 (MQ=255)
                                                                                                                                     cAAGCCCGCGCGGCGAACGTTCTGCATGTGAAAATAAAAAATTATTCAGAGTTGCCCTGGGCGGTCTCGGCGGGGTAGGTCGCGAGCCACAGTTCGAAGCCTCCGTCCAGGCTGTACTCGGCGGAGAAACCCTGTCCTACCAGGTAggc         <  2:297605/149‑1 (MQ=255)
                                                                                                                                         cccGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGgca     <  2:139830/149‑1 (MQ=255)
                                                                                                                                         cccGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTAGCCTACCAGGTAGGCGgca     <  2:344787/149‑1 (MQ=255)
                                                                                                                                          ccGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACg                                               >  7:7142/1‑106 (MQ=255)
                                                                                                                                          ccGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACg                                               <  8:7142/106‑1 (MQ=255)
                                                                                                                                            gcgcgGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGgcagca  >  5:168696/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TACCGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCA  >  NC_002947/484547‑484835

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: