Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,437 T→C 9.7% intergenic (+19/+16) 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,4370TC 9.7% 373.3 / 25.2 135intergenic (+19/+16)ydfZ/ydfIselenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (78/43);  new base C (9/4);  total (88/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 = 3.89e-02

Only 100 of 136 total aligned reads displayed.

GGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCA  >  NC_000913/1629290‑1629586
                                                                                                                                                   |                                                                                                                                                     
 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)
       tCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAtt                                                                                                                                              <  2:1138543/150‑1 (MQ=255)
        ccTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACAGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttc                                                                                                                                             <  1:317724/150‑2 (MQ=255)
           gTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCTCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtt                                                                                                                                          <  2:309773/150‑1 (MQ=255)
             cGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcacacg                                                                                                                                        <  2:1933017/150‑7 (MQ=255)
               cttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGGTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtttttt                                                                                                                                      <  2:1050978/149‑1 (MQ=255)
                ttGATTCTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                                                     <  1:1658609/150‑1 (MQ=255)
                ttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAACCCCGTAGTTGTTGAAGGTTGTGAAGAGCAACTGGCACCACTGGACCTGAGTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                                                     >  2:492637/1‑150 (MQ=255)
                   aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                                  <  1:1670654/150‑1 (MQ=255)
                   aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                                  >  1:1737848/1‑150 (MQ=255)
                   aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGAGTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                                  >  1:1745532/1‑150 (MQ=255)
                    tACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTc                                                                                                                                 >  2:1996939/1‑150 (MQ=255)
                       tGAAGGTCTGACCGCTGAGCAAATCCGCCGCTGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACg                                                                                                                              <  2:1371679/150‑1 (MQ=255)
                            gacTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCGCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAata                                                                                                                         >  1:1094376/3‑150 (MQ=255)
                            gTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAata                                                                                                                         >  2:1817764/1‑150 (MQ=255)
                                 aCCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTg                                                                                                                    >  1:549405/1‑150 (MQ=255)
                                 aCCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTg                                                                                                                    <  2:511593/150‑1 (MQ=255)
                                    gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                                 <  2:1453351/150‑1 (MQ=255)
                                         gCAAATCCGCCTCGGAAAAACCGTAGTTGGTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAAGGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTc                                                                                                            <  1:1570608/150‑1 (MQ=255)
                                         gCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGCGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTc                                                                                                            >  2:118582/1‑150 (MQ=255)
                                               ccgccgCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTc                                                                                                      >  2:626360/1‑150 (MQ=255)
                                                cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcc                                                                                                     >  2:547140/1‑149 (MQ=255)
                                                  ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                   <  1:1748433/150‑1 (MQ=255)
                                                  ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                   <  1:445503/150‑1 (MQ=255)
                                                   cgcgGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCgcgc                                                                                                  >  1:365185/1‑150 (MQ=255)
                                                          aaaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATgcgc                                                                                           >  1:893606/1‑150 (MQ=255)
                                                          aaaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATgcgc                                                                                           >  1:1391256/1‑150 (MQ=255)
                                                           aaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAATCGTGTGAATGCCGCCGATGGCGGCAATGCTTTTTTACTTCAAGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCt                                                                                          >  2:15022/1‑150 (MQ=255)
                                                              cGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaat                                                                                       >  1:1396916/1‑150 (MQ=255)
                                                              cGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaat                                                                                       >  2:832390/1‑150 (MQ=255)
                                                                tAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCtaataa                                                                                     >  1:5482/1‑150 (MQ=255)
                                                                 aGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATaaa                                                                                    >  2:1350810/1‑150 (MQ=255)
                                                                  gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGGGCTAATAAAg                                                                                   >  2:1874449/1‑150 (MQ=255)
                                                                  gttgttGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAg                                                                                   <  2:1541904/150‑1 (MQ=255)
                                                                     gttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTGTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACa                                                                                >  2:1000843/1‑150 (MQ=255)
                                                                       tgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAg                                                                              <  2:249569/150‑1 (MQ=255)
                                                                        gAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCGCGGTATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGAGAAGt                                                                             >  2:789814/1‑150 (MQ=255)
                                                                         aaGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTa                                                                            <  2:477855/150‑1 (MQ=255)
                                                                           ggTTGTGAAGAGACACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTGTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACg                                                                          >  1:17482/1‑150 (MQ=255)
                                                                            gTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGt                                                                         >  1:473872/1‑150 (MQ=255)
                                                                                tgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                                                                     <  1:576345/150‑1 (MQ=255)
                                                                                tgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCCCCATGCGCTAAAAAATACAAGTACGTTTca                                                                     >  2:1726079/1‑149 (MQ=255)
                                                                                tgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTGTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                                                                     >  1:1418411/1‑150 (MQ=255)
                                                                                 gAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATTAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACTGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCg                                                                    <  1:1364180/150‑1 (MQ=255)
                                                                                  aaGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGt                                                                   <  1:1571810/150‑1 (MQ=255)
                                                                                   agagAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                                                  >  1:1528670/1‑150 (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)
                                                                                              caccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGt                                                       >  2:872487/1‑150 (MQ=255)
                                                                                              caccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGt                                                       >  1:1167467/1‑150 (MQ=255)
                                                                                               accacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCACCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTa                                                      >  1:1125980/1‑150 (MQ=255)
                                                                                                ccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCCCTCTTGCAGTaa                                                     >  1:788099/1‑150 (MQ=255)
                                                                                                ccacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaa                                                     <  2:1102026/150‑1 (MQ=255)
                                                                                                  acTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAAGAACGACAAGTACGTTGCCGTCACTCTTGCAGTaacc                                                   >  2:1672516/1‑148 (MQ=255)
                                                                                                     ggACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCAGGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaa                                                <  1:1331465/150‑1 (MQ=255)
                                                                                                        ccTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTAAGTTTCCGTCACTCTTGCAGTAAACAAACTa                                             >  1:204266/1‑150 (MQ=255)
                                                                                                        ccTGAATCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTa                                             >  2:1152875/1‑150 (MQ=255)
                                                                                                           gATTCGTCTCGTCATGAACTAAGCGTGTGAATGCGGCCGATGGCGGGATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCTCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTg                                          <  2:1427621/150‑1 (MQ=255)
                                                                                                              tCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCt                                       <  1:1720417/150‑1 (MQ=255)
                                                                                                                gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                     >  1:1731328/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)
                                                                                                                      gCATGAACTAAGCGTGTGAATGCCGCCGATCGCGGCATTGCTTTTTTACTTCCCGGTAAATTATGCCACGGTCGATTTCGCTCCATGCGCTAATAAAGACAAGTACGTGTCCGTCACTCGTGCAGTAAACAAACTATTGTCTGGCaataa                               >  2:126806/1‑147 (MQ=255)
                                                                                                                         tGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAATCTATTGTCTGGCAAatcatc                            >  2:44097/1‑150 (MQ=255)
                                                                                                                         tGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatcatc                            <  2:893606/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)
                                                                                                                              tAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCaa                       <  2:1525401/150‑1 (MQ=255)
                                                                                                                                aGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAg                     <  1:626360/150‑1 (MQ=255)
                                                                                                                                  cGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGAt                   >  1:229055/1‑150 (MQ=255)
                                                                                                                                    aagcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATGATCACGAAAGATCg                 >  1:440977/5‑150 (MQ=255)
                                                                                                                                    aagcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAATAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCg                 >  2:1643828/5‑150 (MQ=255)
                                                                                                                                     gtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGATCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc                <  1:677150/150‑1 (MQ=255)
                                                                                                                                     gtgAATGCCGCCGATGGCGGCATTGCTTATTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAAACATCACCAAAGATAGc                >  2:689897/1‑150 (MQ=255)
                                                                                                                                     agcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc                >  2:1763629/4‑150 (MQ=255)
                                                                                                                                       caaTGCCGCCATCGGCGGCGTTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGCCAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACGATTGTCTGGCAAATCATCCCGAAAGATCGccg              >  1:830023/2‑149 (MQ=255)
                                                                                                                                       caaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCCTGCGCTAATAAAGACAAGTACGATTCCGTCACTCTTGCAGTAAACACACTATTGTCTGGCAAATCATCACCAAAGATCCCCt              >  1:89981/2‑150 (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       >  1:1871580/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)
                                                                                                                                                 gATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGGCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAg    <  1:492637/150‑1 (MQ=255)
                                                                                                                                                   tGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCa  >  1:1360253/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGCA  >  NC_000913/1629290‑1629586

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: