Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 484,688 C→T 6.9% intergenic (‑250/+33) 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,6880CT 6.9% 253.6 / 6.0 86intergenic (‑250/+33)yeaG/glpEprotein kinase/thiosulfate sulfurtransferase GlpE
Reads supporting (aligned to +/- strand):  ref base C (36/44);  new base T (3/3);  total (39/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 = 6.16e-01

ACGCTACCGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCA  >  NC_002947/484543‑484835
                                                                                                                                                 |                                                                                                                                                    
aCGCTACCGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCcgcgc                                                                                                                                                    <  4:349605/148‑1 (MQ=255)
    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                                                                                                                 >  2:39956/1‑147 (MQ=255)
                                    tcCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAattatt                                                                                                                 >  8:345660/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                                                                                                                   >  4:322826/1‑126 (MQ=255)
                                                       gAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAatta                                                                                                                   <  3:322826/126‑1 (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:307415/149‑1 (MQ=255)
                                                                                                  aaGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCAc                                                <  8:278581/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)
                                                                                                                                                 |                                                                                                                                                    
ACGCTACCGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCA  >  NC_002947/484543‑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: