Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,436 A→T 9.4% intergenic (+18/+17) ydfZ → / ← ydfI selenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,629,4360AT 9.4% 348.8 / 18.6 134intergenic (+18/+17)ydfZ/ydfIselenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (79/42);  new base T (9/4);  total (88/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 = 9.97e-01

Only 100 of 135 total aligned reads displayed.

GGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGC  >  NC_000913/1629290‑1629585
                                                                                                                                                  |                                                                                                                                                     
 gctaGATCCTTTCGATTGATACTTAAGGTCTGACCGCTTAGAATATCCGCCTCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                                                   <  1:42640/147‑1 (MQ=255)
 gCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGCGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                                                   >  2:1776327/1‑150 (MQ=255)
 gCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                                                   >  1:1000209/1‑150 (MQ=255)
 gCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCCCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                                                   >  1:1531657/1‑150 (MQ=255)
 gCAAGATCCTGTCGATTGATACTGAAGGTCGGACCGCGGCGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGGTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGGGTGGGAATGCCGACGATtgc                                                                                                                                                   >  2:1198949/1‑150 (MQ=255)
   aaGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcgg                                                                                                                                                 >  1:1623162/1‑150 (MQ=255)
    aGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcggc                                                                                                                                                >  1:1541904/1‑150 (MQ=255)
       tCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAtt                                                                                                                                             <  2:1138543/150‑1 (MQ=255)
          tGTCGATTGATACTGAAGGTCTGCCCGCTGCGCAAATCCGCCGCGGACAAACCCTAGTTGTTGAAGGTTGTGAAGCGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcac                                                                                                                                          <  1:977249/150‑4 (MQ=255)
             cGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcacacg                                                                                                                                       <  2:1933017/150‑7 (MQ=255)
                ttGATTCTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                                                    <  1:1658609/150‑1 (MQ=255)
                ttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAACCCCGTAGTTGTTGAAGGTTGTGAAGAGCAACTGGCACCACTGGACCTGAGTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                                                    >  2:492637/1‑150 (MQ=255)
                   aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                                 >  1:1737848/1‑150 (MQ=255)
                   aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                                 >  1:1515013/1‑150 (MQ=255)
                   aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGAGTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                                 >  1:1745532/1‑150 (MQ=255)
                     aCTGAAGTTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTCAAGGTTGTGAAGAGAAACTGGCACCACTGGAAATGAGTCGTCTCGGCATGAACTATGCGTGTGAATTACGCCGAGGGCGGCATTGCTTTTTTACttga                                                                                                                               >  2:1282861/1‑148 (MQ=255)
                       tGAAGGTCTGACCGCTGAGCAAATCCGCCGCTGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACg                                                                                                                             <  2:1371679/150‑1 (MQ=255)
                            gacTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCGCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAata                                                                                                                        >  1:1094376/3‑150 (MQ=255)
                            gTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAata                                                                                                                        >  2:1817764/1‑150 (MQ=255)
                               tGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAttt                                                                                                                     >  1:1367988/1‑150 (MQ=255)
                                 tccGCTGAGGAATTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTTAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCTTGTACTTAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTg                                                                                                                   <  1:1728641/149‑1 (MQ=255)
                                 aCCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTg                                                                                                                   <  2:511593/150‑1 (MQ=255)
                                   cGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCACCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGcc                                                                                                                 >  1:1762542/1‑150 (MQ=255)
                                    gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                                <  2:1453351/150‑1 (MQ=255)
                                    gCTGAGCAAATCCGCCGCGGAAAAACCGTAGATGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                                >  2:441021/1‑150 (MQ=255)
                                     cTGAGCAAATCCGCCGCGGAAAAACCCTAGATGTTGAAGGTTGTGAAGAGAAACTGGCACCACTCGACCTGATTCGTCTAGGCATGAACTAAGCGTGTGAATGCCGCCGATGCCGGCATAGCTTTTTTACATCACGGAATATTTTGCCAc                                                                                                               >  1:542388/1‑150 (MQ=255)
                                         gCAAATCCGCCTCGGAAAAACCGTAGTTGGTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAAGGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTc                                                                                                           <  1:1570608/150‑1 (MQ=255)
                                         gCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGCGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTc                                                                                                           >  2:118582/1‑150 (MQ=255)
                                           aaaTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCga                                                                                                         >  1:453927/1‑149 (MQ=255)
                                               ccgccgCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTc                                                                                                     >  2:626360/1‑150 (MQ=255)
                                                cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcc                                                                                                    >  2:547140/1‑149 (MQ=255)
                                                  ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                  <  1:445503/150‑1 (MQ=255)
                                                  ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                  <  1:1748433/150‑1 (MQ=255)
                                                  ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                  <  1:68478/150‑1 (MQ=255)
                                                   cgcgGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCgcgc                                                                                                 >  1:365185/1‑150 (MQ=255)
                                                          aaaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATgcgc                                                                                          >  1:893606/1‑150 (MQ=255)
                                                           aaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAATCGTGTGAATGCCGCCGATGGCGGCAATGCTTTTTTACTTCAAGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCt                                                                                         >  2:15022/1‑150 (MQ=255)
                                                              cGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaat                                                                                      <  2:1645711/150‑1 (MQ=255)
                                                              cGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaat                                                                                      >  2:832390/1‑150 (MQ=255)
                                                              cGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaat                                                                                      >  1:1396916/1‑150 (MQ=255)
                                                                tAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaataa                                                                                    >  1:5482/1‑150 (MQ=255)
                                                                  gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGGGCTAATAAAg                                                                                  >  2:1874449/1‑150 (MQ=255)
                                                                  gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAg                                                                                  <  1:1968552/150‑1 (MQ=255)
                                                                  gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAg                                                                                  <  2:1541904/150‑1 (MQ=255)
                                                                     gttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTGTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACa                                                                               >  2:1000843/1‑150 (MQ=255)
                                                                       tgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAg                                                                             >  2:1434233/1‑150 (MQ=255)
                                                                       tgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAg                                                                             <  2:249569/150‑1 (MQ=255)
                                                                        gAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCGCGGTATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGAGAAGt                                                                            >  2:789814/1‑150 (MQ=255)
                                                                        gAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGt                                                                            >  1:406691/1‑150 (MQ=255)
                                                                         aaGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTa                                                                           <  2:477855/150‑1 (MQ=255)
                                                                            gTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGt                                                                        >  1:473872/1‑150 (MQ=255)
                                                                                tgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                                                                    <  1:576345/150‑1 (MQ=255)
                                                                                tgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTGTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                                                                    >  1:1418411/1‑150 (MQ=255)
                                                                                 gAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATTAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACTGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCg                                                                   <  1:1364180/150‑1 (MQ=255)
                                                                                 gAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATTAAGACAAGTACGTTTCCg                                                                   >  2:394111/1‑150 (MQ=255)
                                                                                  aaGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATCCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCATTCGCGCCATGCGCTAATAAAGCCAAGTACGTTTccct                                                                  >  1:56525/1‑148 (MQ=255)
                                                                                          cttgCACCACTGGAACTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTTCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACTAGTACGTTTCCGTCACTCTTg                                                          <  1:1915912/147‑1 (MQ=255)
                                                                                            ggCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCa                                                        <  1:1996939/150‑1 (MQ=255)
                                                                                            ggCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAAAATTATGCCACGGTCGCTTTCGCGCCATGCGGTAATAAAGACAAGTACGTTTCCGTCACTCTTGCa                                                        >  2:415793/1‑150 (MQ=255)
                                                                                             gCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAg                                                       <  2:131942/150‑1 (MQ=255)
                                                                                             gCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTATCGCCATGCGCTAAAAAAGACAATTACGTTTCCGTCACTCTTGCat                                                       >  1:1022379/1‑149 (MQ=255)
                                                                                             gCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCACCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAg                                                       <  2:1762542/150‑1 (MQ=255)
                                                                                              caccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGt                                                      >  2:872487/1‑150 (MQ=255)
                                                                                              caccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGt                                                      >  1:1167467/1‑150 (MQ=255)
                                                                                               accacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCACCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTa                                                     >  1:1125980/1‑150 (MQ=255)
                                                                                                ccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaa                                                    <  2:1102026/150‑1 (MQ=255)
                                                                                                  acTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAATCCAAGTAAGTTTACGTCACTCTTGCAGTaaac                                                  >  1:1746828/1‑150 (MQ=255)
                                                                                                  acTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAAGAACGACAAGTACGTTGCCGTCACTCTTGCAGTaacc                                                  >  2:1672516/1‑148 (MQ=255)
                                                                                                     ggACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCAGGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaa                                               <  1:1331465/150‑1 (MQ=255)
                                                                                                          tGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTAtt                                          >  2:1385069/1‑150 (MQ=255)
                                                                                                              tCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCt                                      <  1:1720417/150‑1 (MQ=255)
                                                                                                               cGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAACTACGTGTCCGTCACTCGTGCAGTAAACGAACTAGTGTCTg                                     >  2:1064055/1‑150 (MQ=255)
                                                                                                                gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  2:406760/1‑150 (MQ=255)
                                                                                                                gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    <  1:560185/150‑1 (MQ=255)
                                                                                                                  ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCa                                  >  2:1941146/1‑150 (MQ=255)
                                                                                                                     ggCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAat                               >  2:170471/1‑150 (MQ=255)
                                                                                                                         tGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatcatc                           <  2:893606/150‑1 (MQ=255)
                                                                                                                            aCTAAGCGTGTGAATGCCGCCGATGTCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCAcc                        <  1:843483/150‑1 (MQ=255)
                                                                                                                            aCTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCCCTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCAcc                        >  2:1423681/1‑150 (MQ=255)
                                                                                                                            aCTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCAcc                        >  1:1106271/1‑150 (MQ=255)
                                                                                                                             cTAAGCGTGTGAATGCCGCCGATGGCGGCAATGCTTTTTTACTTCACGGAATAATTTGCCACGGTCGCTTTCGCGCCATGCGCTAAGAAAGACAAGTACCTTTCCGTCACTGTTGCAGTAAACAAACTATTGTCTGGCATATCATCTCCa                       >  2:153510/1‑150 (MQ=255)
                                                                                                                              tAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCCTGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTATACAAACTATTGTCTGGCAAATCATCACCaa                      >  1:266658/1‑150 (MQ=255)
                                                                                                                                  cGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGAt                  >  1:229055/1‑150 (MQ=255)
                                                                                                                                   gtgtGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATc                 <  1:1637818/150‑1 (MQ=255)
                                                                                                                                     gtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGATCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc               <  1:677150/150‑1 (MQ=255)
                                                                                                                                     gtgAATGCCGCCGATGGCGGCATTGCTTATTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAAACATCACCAAAGATAGc               >  2:689897/1‑150 (MQ=255)
                                                                                                                                     agcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc               >  2:1763629/4‑150 (MQ=255)
                                                                                                                                      gcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCAGTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGcc              >  2:24447/3‑150 (MQ=255)
                                                                                                                                       caaTGCCGCCATCGGCGGCGTTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGCCAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACGATTGTCTGGCAAATCATCCCGAAAGATCGccg             >  1:830023/2‑149 (MQ=255)
                                                                                                                                        aaTGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCCCTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCtt            >  2:837125/1‑150 (MQ=255)
                                                                                                                                        aaTGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCtt            >  2:294791/1‑150 (MQ=255)
                                                                                                                                        aaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTACACAAACGATTGGCTGGCAAATCATCACCAAAGATCGCCtt            >  2:601799/1‑150 (MQ=255)
                                                                                                                                        aaTGCCGCCATCGGCGGCATTGCTTTTTTACTACACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGGCACTCTTGCAGTAAACAAACTATTGGCTGGCAAATCATAACCAAAGATCGCCtt            >  2:1919390/1‑150 (MQ=255)
                                                                                                                                           gccgccGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAAt         >  1:714193/1‑150 (MQ=255)
                                                                                                                                           gccgccGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAAt         >  2:561811/1‑150 (MQ=255)
                                                                                                                                              gccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGc      >  2:1182080/1‑150 (MQ=255)
                                                                                                                                                cgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCa    <  1:1817764/150‑1 (MQ=255)
                                                                                                                                                cgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCa    >  2:1021998/1‑150 (MQ=255)
                                                                                                                                                 gATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCATTAATAGCCAg   >  1:296349/1‑150 (MQ=255)
                                                                                                                                                  aTGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGc  <  2:1528670/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGC  >  NC_000913/1629290‑1629585

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