Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,436 A→T 6.1% 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 6.1% 404.8 / 11.9 142intergenic (+18/+17)ydfZ/ydfIselenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (79/53);  new base T (3/6);  total (83/59)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

Only 100 of 144 total aligned reads displayed.

ACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGC  >  NC_000913/1629287‑1629585
                                                                                                                                                     |                                                                                                                                                     
acTGGCAAGACCCTGTCGATTCATACTGAAGGTGTGACGGGTGAGCAAATCCGCCGCGGAAAAACCGTTGTTGTTGAAGGTTGTGAAGAGAAAATGTCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGGGTGTGAATGCCGCCGa                                                                                                                                                       <  2:1491124/150‑1 (MQ=255)
 cTGGCAAGATCCTGCCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAt                                                                                                                                                      >  2:33497/1‑150 (MQ=255)
    gCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggc                                                                                                                                                   >  2:1271256/1‑150 (MQ=255)
     cAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcg                                                                                                                                                  >  1:518731/1‑150 (MQ=255)
     cAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcg                                                                                                                                                  >  2:1421730/1‑150 (MQ=255)
          tCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCAtt                                                                                                                                             >  2:2015130/1‑150 (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)
               ccGATGGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCAGGAACTACGCGTGTGAATGCCGCCGATGGCGGCATTGCttt                                                                                                                                        <  1:466929/149‑1 (MQ=255)
                cGATTGATACTGAAGGGCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAATGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATAGGCGGCAttcacacg                                                                                                                                       <  2:414068/150‑7 (MQ=255)
                 gATTGATACTGAAGGTCTGACCGCTGAGCAAATCCTCCGCTGAAAAAACGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGTACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttttt                                                                                                                                      <  2:840867/150‑1 (MQ=255)
                 gATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCttttt                                                                                                                                      >  1:1074858/1‑150 (MQ=255)
                  aTTGATACTGAAGGTCAGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtttttt                                                                                                                                     <  1:971792/150‑1 (MQ=255)
                   ttGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTa                                                                                                                                    <  1:1139215/150‑1 (MQ=255)
                      aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAAGTCGGCATTGCTTTTATACtt                                                                                                                                 >  1:657901/1‑150 (MQ=255)
                      aTACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGGCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACtt                                                                                                                                 <  2:501755/150‑1 (MQ=255)
                        aCTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCa                                                                                                                               <  2:1609881/150‑1 (MQ=255)
                           gAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACgg                                                                                                                            >  2:1745193/1‑150 (MQ=255)
                              ggTCTGACCGCTGAGCAAATCCGCCGCGGTAAACCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTCGCACCACTGGCCCTGATTCGTCTCGGCAGGAACTAAGCGTGTGAATGCCGCCGCTGGCGGCAGTGCTTTGTTACTTCACGGAat                                                                                                                         >  1:316738/1‑150 (MQ=255)
                              ggTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGCCCTGATTCGTCGCGGCAGGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACTGCACGGAat                                                                                                                         >  1:1263710/1‑150 (MQ=255)
                                tCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAatat                                                                                                                       <  1:169294/150‑1 (MQ=255)
                                   gACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGGCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATAtttt                                                                                                                    >  1:1412724/1‑150 (MQ=255)
                                    aCCGCTGATCAAATCCGCCGCGGAAAACCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTg                                                                                                                   <  1:1001636/150‑1 (MQ=255)
                                     ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGAATGAACTAATCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTAACGGAATATTTTGc                                                                                                                  >  2:382190/1‑150 (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:1480907/150‑1 (MQ=255)
                                            gCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTc                                                                                                           >  1:302879/1‑150 (MQ=255)
                                             cAACTCCGCCGCTGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGTCCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCg                                                                                                          <  2:216002/150‑1 (MQ=255)
                                              aaaTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGCAATATTTTGCCACGGTCGc                                                                                                         >  2:1659849/1‑150 (MQ=255)
                                               aaTCCTCCGCGTAAATACAGTAGTTGTTGAAGTTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCt                                                                                                        <  2:362371/150‑1 (MQ=255)
                                               aaTCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCt                                                                                                        >  1:571693/1‑150 (MQ=255)
                                                   cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                    >  2:947027/1‑150 (MQ=255)
                                                   cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                    >  1:1757082/1‑150 (MQ=255)
                                                   cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                    >  2:2044819/1‑150 (MQ=255)
                                                   cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttca                                                                                                                                           >  1:380308/1‑109 (MQ=255)
                                                     ccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcgcg                                                                                                  >  2:1997061/1‑150 (MQ=255)
                                                      cgcgGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTTCTTCACGGAATATTTTGCCACGGTCGCTTTCgcgc                                                                                                 >  1:590822/1‑150 (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)
                                                                           gAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGt                                                                            >  1:1432896/1‑150 (MQ=255)
                                                                              ggTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACg                                                                         <  2:477533/150‑1 (MQ=255)
                                                                                ttGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGtt                                                                       >  1:329658/1‑150 (MQ=255)
                                                                                    taaGAGAAACTGGCAACAGTGGACCTGATTCGTCTAGTCATGAACTAAGCGTGTGAATGCCGCCGATGTAGGCATTGCTTTTTTACTTCACGTAATATTTTGACACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCg                                                                   <  1:559688/149‑1 (MQ=255)
                                                                                     aaGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGt                                                                  <  1:1976189/150‑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)
                                                                                         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)
                                                                                                       tGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaa                                                <  1:759145/150‑1 (MQ=255)
                                                                                                        ggACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaa                                               >  2:1041568/1‑150 (MQ=255)
                                                                                                         gACCTGATTCGTCTCGGCATGAACTAAGCGTCTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTATCGCGCCATGCGCTAATAAAGACAAGTAGGTTTCCGTCACTCTTGCAGTaaacaaac                                              >  1:19420/1‑150 (MQ=255)
                                                                                                          aCCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCCCGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACt                                             >  2:1659793/1‑150 (MQ=255)
                                                                                                           ccTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATGAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTa                                            >  2:143434/1‑150 (MQ=255)
                                                                                                             tGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTAtt                                          <  2:16968/150‑1 (MQ=255)
                                                                                                                 tCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCt                                      <  2:309969/150‑1 (MQ=255)
                                                                                                                   gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    <  1:2006668/150‑1 (MQ=255)
                                                                                                                   gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  1:197458/1‑150 (MQ=255)
                                                                                                                   gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  2:26083/1‑150 (MQ=255)
                                                                                                                   gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                    >  1:936550/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)
                                                                                                                     ctcGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCa                                  >  1:801822/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)
                                                                                                                          cATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGAGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAatca                             >  2:1137661/1‑150 (MQ=255)
                                                                                                                            tGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAACCAAACTATTGTCTGGCAAatcatc                           >  1:602242/1‑150 (MQ=255)
                                                                                                                             gAACTAAGCGTGTGAACGCCGCCGATGGCGGCATTGCTTTTTTACTTCACTGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAAtcatca                          <  2:371678/150‑1 (MQ=255)
                                                                                                                                cTAAGCGTGTGAGTGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCa                       <  1:4838/150‑1 (MQ=255)
                                                                                                                                   aGCGTGTGAATGCCGGCGATGGCGGCATTGCTTTTTTACTTCACGGATTATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTACGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAg                    <  2:1381748/150‑1 (MQ=255)
                                                                                                                                   aGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGTTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAg                    <  2:302879/150‑1 (MQ=255)
                                                                                                                                    gggtgtGACTGCCGCCGATGGAGGCATTGCTTTTTTACTTCACGGAATATTTTGCCAAGGTCGCTTTCGCGCCATGCGCTAATAAAGACATGTACGTTTCCGTCACTCTTGCAGTAAACAACCTATTGTCTGGCAAATCATCACCAAAGa                   <  2:1263709/148‑1 (MQ=255)
                                                                                                                                      ttgtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATc                 >  1:329615/2‑150 (MQ=255)
                                                                                                                                       aagcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCg                >  1:1906548/5‑150 (MQ=255)
                                                                                                                                        gtgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc               <  1:990949/150‑1 (MQ=255)
                                                                                                                                        agcaaTGCCGCCATCGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGc               >  2:1173329/4‑150 (MQ=255)
                                                                                                                                         tgAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGGCACTCTTGCAGTAAACACACTATTGTCTGGCAAATCATCACCAAAGATCGcc              >  1:1793662/1‑150 (MQ=255)
                                                                                                                                          gAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGGTTCCGTCACTATTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCt             >  1:1395875/1‑150 (MQ=255)
                                                                                                                                                  ccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGcc     <  1:850992/150‑1 (MQ=255)
                                                                                                                                                    gATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAg   <  2:1050538/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAGC  >  NC_000913/1629287‑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: