Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 484,699 A→G 7.5% intergenic (‑258/+25) 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,6990AG 7.5% 640.7 / 26.1 211intergenic (‑258/+25)yeaG/glpEprotein kinase/thiosulfate sulfurtransferase GlpE
Reads supporting (aligned to +/- strand):  ref base A (93/102);  new base G (9/7);  total (102/109)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 215 total aligned reads displayed.

TACCGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGG  >  NC_002947/484550‑484848
                                                                                                                                                     |                                                                                                                                                     
   cGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTGTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACgg                                                                                                                                                    >  2:1919192/1‑150 (MQ=255)
      tGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCt                                                                                                                                                 <  2:2925005/150‑1 (MQ=255)
      tGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCt                                                                                                                                                 >  2:118185/1‑150 (MQ=255)
      tGACCCGCACGCCGGTGTACCGGCTCTCCCCGTTTTGGATGGCCTGAGCTTAAGCGTAGGTGGGTATCTGCACGGTAAAGGGGTGAGAAAGGCAGACGTTCCGAACGACCGGTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGcg                                                                                                                                                 >  2:3260165/1‑149 (MQ=255)
        aCCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttg                                                                                                                                               <  2:1337717/150‑3 (MQ=255)
        aCCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttg                                                                                                                                               <  1:423207/150‑3 (MQ=255)
        aCCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGc                                                                                                                                               <  1:5313279/150‑1 (MQ=255)
          ccGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTACGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttgcg                                                                                                                                             <  2:4365162/150‑5 (MQ=255)
          ccGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAg                                                                                                                                             <  1:4536898/150‑1 (MQ=255)
           ggagagCCGGTACACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAgg                                                                                                                                            >  2:2356812/6‑150 (MQ=255)
            gCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAggg                                                                                                                                           >  2:4937644/1‑150 (MQ=255)
            gCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAggg                                                                                                                                           <  1:2177471/150‑1 (MQ=255)
            gCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAggg                                                                                                                                           >  2:2171703/1‑150 (MQ=255)
                gCCGGTGTACCGGCTCTCCCCTTGTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGaaa                                                                                                                                       >  2:4989170/1‑150 (MQ=255)
                 ccGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGGAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGaaaa                                                                                                                                      <  2:2158193/150‑1 (MQ=255)
                  cGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAAt                                                                                                                                     >  2:69113/1‑150 (MQ=255)
                      gtACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATaaaa                                                                                                                                 <  2:768092/150‑1 (MQ=255)
                           ggCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAatta                                                                                                                            <  1:2400714/150‑1 (MQ=255)
                                tcCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAg                                                                                                                       >  1:1325273/1‑150 (MQ=255)
                                  cccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttgc                                                                                                                                              <  2:1263977/125‑4 (MQ=255)
                                  cccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCg                                                                                                                     >  1:3182394/1‑150 (MQ=255)
                                  cccTTTTTGGATGGCCTGAGCTAAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCg                                                                                                                     >  1:1080781/1‑150 (MQ=255)
                                   ccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTg                                                                                                                                                >  1:3110964/1‑122 (MQ=255)
                                   ccTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTg                                                                                                                                                <  2:3110964/122‑1 (MQ=255)
                                     tttttGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTg                                                                                                                  <  2:3379027/150‑1 (MQ=255)
                                     tttttGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTg                                                                                                                  >  2:2879382/1‑150 (MQ=255)
                                     cttttGGAGGGCCTGAGCTTAAGCGTAGTTGGGTATCGGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGAGAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTg                                                                                                                  <  2:1736777/149‑1 (MQ=255)
                                                ccTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGcc                                                                                                                <  2:3819992/141‑1 (MQ=255)
                                                ccTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGt                                                                                                       <  1:2208321/150‑1 (MQ=255)
                                                 cTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGtc                                                                                                      >  1:5238606/1‑150 (MQ=255)
                                                   gAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGtctc                                                                                                    <  1:3456156/150‑1 (MQ=255)
                                                    aGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCg                                                                                                   >  2:1726010/1‑150 (MQ=255)
                                                       ttAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGTGCGCCGCGCGGGCttgc                                                                                                                                              <  2:5248233/104‑4 (MQ=255)
                                                       ttAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGCAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGcc                                                                                                >  1:2712364/1‑150 (MQ=255)
                                                       ttAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGcc                                                                                                >  1:5234598/1‑150 (MQ=255)
                                                       ttAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGcc                                                                                                >  1:1253756/1‑150 (MQ=255)
                                                        tAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCg                                                                                               >  1:787590/1‑150 (MQ=255)
                                                         aaGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCgg                                                                                              <  1:4517621/150‑1 (MQ=255)
                                                              tAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAgg                                                                                         <  2:1500301/150‑1 (MQ=255)
                                                               aGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGt                                                                                        <  1:3015670/150‑1 (MQ=255)
                                                                gTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTc                                                                                       <  1:4722647/150‑1 (MQ=255)
                                                                 ttGTTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCCCAAGCCCGCGCGGCGCACGGGATGCAGCGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCTGGTAGGTcg                                                                                      >  2:3696586/1‑150 (MQ=255)
                                                                   ggTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGGTACGAACGACCGGTCGTCAGCTCAGCCGCACGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCGGGGCGGTCTCGGCCGGGTAGGTcgcg                                                                                    <  2:4770306/150‑1 (MQ=255)
                                                                    gTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTcgcgc                                                                                   <  1:3245655/150‑1 (MQ=255)
                                                                       aTCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGcc                                                                                <  1:3516801/150‑1 (MQ=255)
                                                                         cTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCcac                                                                              <  2:4104401/150‑1 (MQ=255)
                                                                          tGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCcaca                                                                             >  2:2325933/1‑150 (MQ=255)
                                                                          tGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCcaca                                                                             <  2:5157675/150‑1 (MQ=255)
                                                                          tGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCcaca                                                                             >  1:2537211/1‑150 (MQ=255)
                                                                           gCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTgg                                                                                                            <  2:860904/118‑1 (MQ=255)
                                                                           gCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTgg                                                                                                            >  1:860904/1‑118 (MQ=255)
                                                                            cAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGt                                                                           >  2:4309770/1‑150 (MQ=255)
                                                                            cAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGt                                                                           >  2:3706280/1‑150 (MQ=255)
                                                                               ggTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCg                                                                        >  1:3539433/1‑150 (MQ=255)
                                                                               ggTAAACGGGTGAGAAAGGCAGAAGTAACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAAGTATACAGCGATGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGATCg                                                                        >  2:842499/1‑150 (MQ=255)
                                                                                 tAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                      <  1:2754352/150‑1 (MQ=255)
                                                                                 tAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                      <  2:1948598/150‑1 (MQ=255)
                                                                                 tAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                      <  1:3788311/150‑1 (MQ=255)
                                                                                 tAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCAGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGaa                                                                      <  2:3939294/150‑1 (MQ=255)
                                                                                   aaGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGc                                                                    <  1:1028243/150‑1 (MQ=255)
                                                                                   aaGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGc                                                                    >  1:4702474/1‑150 (MQ=255)
                                                                                    aGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGcc                                                                   >  2:2686714/1‑150 (MQ=255)
                                                                                          gagaAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTc                                                             <  2:5099034/150‑1 (MQ=255)
                                                                                             aaaGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTCGGTCGCGCGCCCCAGTTCGAAGCCTCCGTcccc                                                          >  2:1074166/1‑148 (MQ=255)
                                                                                                  cAGAAGTTACGAACGACCGGTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAACATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGt                                                     <  2:1559626/150‑1 (MQ=255)
                                                                                                        ttACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGt                                               >  2:4605812/1‑150 (MQ=255)
                                                                                                            gaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGga                                           <  1:4332369/150‑1 (MQ=255)
                                                                                                            gaacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGga                                           <  2:486730/150‑1 (MQ=255)
                                                                                                             aacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGgag                                          <  1:3919461/150‑1 (MQ=255)
                                                                                                             aacgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGgag                                          >  1:3314285/1‑150 (MQ=255)
                                                                                                              acgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAgg                                                                                                                                            >  2:1150458/1‑51 (MQ=255)
                                                                                                              acgaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAgg                                                                                                                                            <  1:1150458/51‑1 (MQ=255)
                                                                                                                gaCCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGaaa                                       >  1:4057173/1‑150 (MQ=255)
                                                                                                                  ccGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAAcc                                     <  2:3048697/150‑1 (MQ=255)
                                                                                                                     ttCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTg                                  <  2:2598904/150‑1 (MQ=255)
                                                                                                                      tCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTgg                                 <  1:2879382/150‑1 (MQ=255)
                                                                                                                         taaGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCt                              >  2:4493374/3‑150 (MQ=255)
                                                                                                                         tCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGcca                              >  1:4366430/1‑149 (MQ=255)
                                                                                                                             aTAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCGGGGGAAAATTATAACTTATTCTGCGTTGCCGTGGGCGGTCTCGGCCGGGTAGGTCGGGCGCGGCGGTTCGAAGCCTCCGTCCACGCGGTACGCGTCGGAGATACCCTGGCGTAccg                          >  2:1751378/1‑149 (MQ=255)
                                                                                                                             aTAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCa                          <  1:5115900/150‑1 (MQ=255)
                                                                                                                                 cctgCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTa                      >  1:1897979/4‑150 (MQ=255)
                                                                                                                                  ccGCAAGCCCGCGCGGCGCACGGGCTGCAGGTGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAg                     >  2:2466538/1‑150 (MQ=255)
                                                                                                                                  ccGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAg                     >  2:3590907/1‑150 (MQ=255)
                                                                                                                                  ccGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAg                     >  2:2899922/1‑150 (MQ=255)
                                                                                                                                    gCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAggc                   >  1:1061388/1‑150 (MQ=255)
                                                                                                                                        gCCCGTGCGCCGCGCGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGGTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGgca               >  2:299918/1‑150 (MQ=255)
                                                                                                                                        gCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGgca               <  1:4937644/150‑1 (MQ=255)
                                                                                                                                         cccGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGgcag              >  1:1423531/1‑150 (MQ=255)
                                                                                                                                          ccGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGgcatc             >  1:2023444/1‑148 (MQ=255)
                                                                                                                                          ccGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGgcagc             >  2:5327808/1‑150 (MQ=255)
                                                                                                                                           cgcgcgGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGgcagca            >  2:3423667/1‑150 (MQ=255)
                                                                                                                                           cgagcgccgCGCGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGACCTGGGCGGTCTCGGCAGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGACCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGTTAGGAGgcagca            >  2:103137/4‑150 (MQ=255)
                                                                                                                                             cgcgGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACt          >  1:3820422/1‑150 (MQ=255)
                                                                                                                                                ggcgCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTg       <  1:4105349/150‑1 (MQ=255)
                                                                                                                                                  cgcACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGc     <  1:2155662/150‑1 (MQ=255)
                                                                                                                                                  cgcACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGc     <  2:3552245/150‑1 (MQ=255)
                                                                                                                                                   gcACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCt    >  1:2114909/1‑150 (MQ=255)
                                                                                                                                                   gcACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCt    <  2:4049631/150‑1 (MQ=255)
                                                                                                                                                     aCGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTgg  <  1:69113/150‑1 (MQ=255)
                                                                                                                                                     aCGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTgg  >  1:4783708/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACCGATGACCCGCACGCCGGTGTACCGGCTCTCCCCTTTTTGGATGGCCTGAGCTTAAGCGTAGTTGGTTATCTGCAAGGTAAAGGGGTGAGAAAGGCAGAAGTTACGAACGACCGTTCGTCAGATAAGCCGCAAGCCCGCGCGGCGCACGGGCTGCAGGGGAAAATAAAAAATTATTCAGCGTTGCCCTGGGCGGTCTCGGCCGGGTAGGTCGCGCGCCACAGTTCGAAGCCTCCGTCCACGCTGTACACGTCGGAGAAACCCTGGCCTACCAGGTAGGCGGCAGCACTTTGGCTGG  >  NC_002947/484550‑484848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: