Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,437 T→C 6.6% 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 6.6% 412.3 / 13.4 142intergenic (+19/+16)ydfZ/ydfIselenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (74/54);  new base C (4/5);  total (82/60)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 143 total aligned reads displayed.

CTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGC  >  NC_000913/1629288‑1629585
                                                                                                                                                     |                                                                                                                                                    
cTGGCAAGATCCTGCCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAt                                                                                                                                                      >  2:33497/1‑150 (MQ=255)
 tGGCACGATCCTGTCGATTGATACTGAAGGGCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGGGAAGAGAAACTGGCACCACTGGACCTGCTTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATg                                                                                                                                                     <  2:1784780/150‑1 (MQ=255)
         tCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAtt                                                                                                                                             >  2:2015130/1‑150 (MQ=255)
          ccTGTCGATTGATACTGAATGTCTGACCGCGGAGCAAATCCGCCGCGGAAAAACCGTAGTTGATGAAGGTTGTGATGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttc                                                                                                                                            <  1:498576/150‑2 (MQ=255)
          ccTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttc                                                                                                                                            <  2:1843259/150‑2 (MQ=255)
            tGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCt                                                                                                                                          <  1:1905390/150‑1 (MQ=255)
             gTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcaca                                                                                                                                         <  1:1782730/150‑5 (MQ=255)
             gTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGAAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtt                                                                                                                                         >  1:2001238/1‑150 (MQ=255)
               cGATTGATACTGAAGGGCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAATGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATAGGCGGCAttcacacg                                                                                                                                       <  2:414068/150‑7 (MQ=255)
                gATTGATACTGAAGGTCTGACCGCTGAGCAAATCCTCCGCTGAAAAAACGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGTACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttttt                                                                                                                                      <  2:840867/150‑1 (MQ=255)
                gATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttttt                                                                                                                                      >  1:1074858/1‑150 (MQ=255)
                  ttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                                                    <  1:1139215/150‑1 (MQ=255)
                     aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGGCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                                 <  2:501755/150‑1 (MQ=255)
                       aCTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCa                                                                                                                               <  2:1609881/150‑1 (MQ=255)
                             ggTCTGACCGCTGAGCAAATCCGCCGCGGTAAACCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTCGCACCACTGGCCCTGATTCGTCTCGGCAGGAACTAAGCGTGTGAATGCCGCCGCTGGCGGCAGTGCTTTGTTACTTCACGGAat                                                                                                                         >  1:316738/1‑150 (MQ=255)
                             ggTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGCCCTGATTCGTCGCGGCAGGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACTGCACGGAat                                                                                                                         >  1:1263710/1‑150 (MQ=255)
                                  gACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGGCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAtttt                                                                                                                    >  1:1412724/1‑150 (MQ=255)
                                   aCCGCTGATCAAATCCGCCGCGGAAAACCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTg                                                                                                                   <  1:1001636/150‑1 (MQ=255)
                                    ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGc                                                                                                                  <  1:570263/150‑1 (MQ=255)
                                    ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAAGGCGGCATAGCTTTTTTACTTCACGGAATATTTTGc                                                                                                                  >  2:1088110/1‑150 (MQ=255)
                                    ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGAATGAACTAATCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTAACGGAATATTTTGc                                                                                                                  >  2:382190/1‑150 (MQ=255)
                                     cgcgGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGcc                                                                                                                 <  2:1465975/150‑1 (MQ=255)
                                     cGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAACCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAGGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGcc                                                                                                                 >  2:279949/1‑150 (MQ=255)
                                     cGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGCCCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGACGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGcc                                                                                                                 >  2:483716/1‑150 (MQ=255)
                                      gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                                >  2:990948/1‑150 (MQ=255)
                                      gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                                <  1:1975084/150‑1 (MQ=255)
                                      gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                                <  1:1480907/150‑1 (MQ=255)
                                          aGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCGGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGACGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGt                                                                                                            >  1:326093/1‑150 (MQ=255)
                                             aaaTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGCAATATTTTGCCACGGTCGc                                                                                                         >  2:1659849/1‑150 (MQ=255)
                                              aaTCCTCCGCGTAAATACAGTAGTTGTTGAAGTTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCt                                                                                                        <  2:362371/150‑1 (MQ=255)
                                                  cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                    >  2:947027/1‑150 (MQ=255)
                                                  cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                    >  1:1757082/1‑150 (MQ=255)
                                                  cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                    >  2:2044819/1‑150 (MQ=255)
                                                  cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttca                                                                                                                                           <  2:380308/111‑3 (MQ=255)
                                                  cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttca                                                                                                                                           >  1:380308/1‑109 (MQ=255)
                                                    ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                  >  2:1997061/1‑150 (MQ=255)
                                                    ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                  >  1:460605/1‑150 (MQ=255)
                                                     cgcgGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTTCTTCACGGAATATTTTGCCACGGTCGCTTTCgcgc                                                                                                 >  1:590822/1‑150 (MQ=255)
                                                       cgGAAAAACCGTAGTTGTTGCAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCa                                                                                               <  1:639821/150‑1 (MQ=255)
                                                        ggAAAAGCCGTAGGTGTTGAAGGTTGTCCAGAGAACCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACGTCAATGAATATTTTGCAACGGTCGCTGTCGCGcccc                                                                                              >  2:741855/1‑148 (MQ=255)
                                                         gAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATg                                                                                             >  1:1129132/1‑150 (MQ=255)
                                                         gAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCAGATTCGTCTCGGCATGAACTTAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATg                                                                                             <  1:1239489/150‑1 (MQ=255)
                                                           aaaaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATgcg                                                                                           >  2:475198/1‑150 (MQ=255)
                                                            aaaCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGAGGCATTGCTTTTTTACATCACGGAATATTTTGCAACGGTCGCTTTCGCGCCATgcgc                                                                                          >  1:1196578/1‑150 (MQ=255)
                                                             aaCCGTAGTTGTTGACTGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGTCATGAACTATGCGGGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCt                                                                                         <  2:484498/150‑1 (MQ=255)
                                                              aCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCCCCACTGGCCCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCCCGGTCGCTTTCGCGCCATGCGCta                                                                                        >  1:1246136/1‑150 (MQ=255)
                                                                     ttgttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATaaata                                                                                 >  2:1278091/1‑148 (MQ=255)
                                                                      tgttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGTATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACTGTCGCTTTCGCGCCATGCGCTAATAAAGAc                                                                                >  2:1856254/1‑150 (MQ=255)
                                                                       gttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACa                                                                               >  1:770367/1‑150 (MQ=255)
                                                                          gAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGt                                                                            >  1:1432896/1‑150 (MQ=255)
                                                                             ggTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACg                                                                         <  2:477533/150‑1 (MQ=255)
                                                                               ttGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGtt                                                                       >  1:329658/1‑150 (MQ=255)
                                                                                 gtgAAGAGAAACTGGCACCACTGGACCTGCTTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTc                                                                     <  1:167878/150‑1 (MQ=255)
                                                                                  tgAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTcc                                                                    >  1:285820/1‑150 (MQ=255)
                                                                                   taaGAGAAACTGGCAACAGTGGACCTGATTCGTCTAGTCATGAACTAAGCGTGTGAATGCCGCCGATGTAGGCATTGCTTTTTTACTTCACGTAATATTTTGACACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCg                                                                   <  1:559688/149‑1 (MQ=255)
                                                                                     agagAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGCACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTGACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                                                 >  1:1085143/1‑150 (MQ=255)
                                                                                     agagAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGTGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                                                 <  2:712712/150‑1 (MQ=255)
                                                                                     agagAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTc                                                                 >  1:277971/1‑150 (MQ=255)
                                                                                      gagaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCa                                                                >  1:866101/1‑150 (MQ=255)
                                                                                       agaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCAc                                                               <  2:1407465/150‑1 (MQ=255)
                                                                                        gaAACTGGCACCCCTGGACCTGATTAGTCTCGGGATGAACTAAGCCTGTTAAAGCCCCCGAAGGCGGGATTGCGTTTTTACATCAAGGAATATTTTGCCACGGTCGCGTAAGCGCCAAGCGGGAAAAAATAAAAGTACGTTTCCGTCAca                                                              >  2:215661/1‑149 (MQ=255)
                                                                                          aaCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCActct                                                            >  1:273641/1‑150 (MQ=255)
                                                                                          aaCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGACACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCActct                                                            >  1:1588937/1‑150 (MQ=255)
                                                                                            cTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTg                                                          <  1:1629379/150‑1 (MQ=255)
                                                                                                   cacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaa                                                   <  2:818247/150‑1 (MQ=255)
                                                                                                   cacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaa                                                   >  2:593253/1‑150 (MQ=255)
                                                                                                    acTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaac                                                  <  2:1700218/150‑1 (MQ=255)
                                                                                                      tGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaa                                                <  1:832573/150‑1 (MQ=255)
                                                                                                      tGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaa                                                <  1:759145/150‑1 (MQ=255)
                                                                                                         aCCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCCCGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACt                                             >  2:1659793/1‑150 (MQ=255)
                                                                                                          ccTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATGAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTa                                            >  2:143434/1‑150 (MQ=255)
                                                                                                            tGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTAtt                                          >  1:1915585/1‑150 (MQ=255)
                                                                                                                tCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCt                                      <  2:309969/150‑1 (MQ=255)
                                                                                                                  gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  1:936550/1‑150 (MQ=255)
                                                                                                                  gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    <  1:2006668/150‑1 (MQ=255)
                                                                                                                  gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  2:242723/1‑150 (MQ=255)
                                                                                                                  gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  2:26083/1‑150 (MQ=255)
                                                                                                                  gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  1:197458/1‑150 (MQ=255)
                                                                                                                  gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACGCTTGCAGTAAACAAACTATTGTCTgg                                    <  1:33579/150‑1 (MQ=255)
                                                                                                                  gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTATTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  1:80632/1‑150 (MQ=255)
                                                                                                                   tctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGc                                   <  2:1401909/150‑1 (MQ=255)
                                                                                                                    ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTTCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTAACTCTTGCAGTAAACAAACTATTGTCTGGCa                                  >  2:628468/1‑150 (MQ=255)
                                                                                                                    ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCa                                  >  2:1710635/1‑150 (MQ=255)
                                                                                                                        gCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCTTTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTTATAAAGACAAGTACGGTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatc                              >  1:1333129/1‑150 (MQ=255)
                                                                                                                        gCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatc                              <  1:382190/150‑1 (MQ=255)
                                                                                                                         cATGAACTCAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatca                             <  1:1813459/150‑1 (MQ=255)
                                                                                                                           tGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAACCAAACTATTGTCTGGCAAatcatc                           >  1:602242/1‑150 (MQ=255)
                                                                                                                               cTAAGCGTGTGAGTGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCa                       <  1:4838/150‑1 (MQ=255)
                                                                                                                                  aGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGTTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAg                    <  2:302879/150‑1 (MQ=255)
                                                                                                                                   gggtgtGACTGCCGCCGATGGAGGCATTGCTTTTTTACTTCACGGAATATTTTGCCAAGGTCGCTTTCGCGCCATGCGCTAATAAAGACATGTACGTTTCCGTCACTCTTGCAGTAAACAACCTATTGTCTGGCAAATCATCACCAAAGa                   <  2:1263709/148‑1 (MQ=255)
                                                                                                                                      tgtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGCTCg                >  2:1185418/1‑150 (MQ=255)
                                                                                                                                      aagcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCg                >  1:1906548/5‑150 (MQ=255)
                                                                                                                                       gtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc               <  1:990949/150‑1 (MQ=255)
                                                                                                                                       agcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc               >  2:1173329/4‑150 (MQ=255)
                                                                                                                                        tgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGcc              >  1:907162/1‑150 (MQ=255)
                                                                                                                                        tgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGGCACTCTTGCAGTAAACACACTATTGTCTGGCAAATCATCACCAAAGATCGcc              >  1:1793662/1‑150 (MQ=255)
                                                                                                                                        tgAAAGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGcc              <  2:1008780/150‑1 (MQ=255)
                                                                                                                                         gAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGGTTCCGTCACTATTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCt             >  1:1395875/1‑150 (MQ=255)
                                                                                                                                                   gATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAg   <  2:1050538/150‑1 (MQ=255)
                                                                                                                                                    aTGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGc  <  2:131653/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGC  >  NC_000913/1629288‑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: