Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,629,435 G→A 6.4% intergenic (+17/+18) 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,4350GA 6.4% 434.7 / 11.8 140intergenic (+17/+18)ydfZ/ydfIselenoprotein, function unknown/putative NAD‑dependent D‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (79/52);  new base A (3/6);  total (82/58)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

Only 100 of 147 total aligned reads displayed.

CACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAG  >  NC_000913/1629286‑1629584
                                                                                                                                                     |                                                                                                                                                     
cacTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGccgccg                                                                                                                                                       >  1:16968/1‑150 (MQ=255)
cacTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGccgccg                                                                                                                                                       <  1:1165089/150‑1 (MQ=255)
cacTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGccgccg                                                                                                                                                       >  1:131653/1‑150 (MQ=255)
 acTGGCAAGACCCTGTCGATTCATACTGAAGGTGTGACGGGTGAGCAAATCCGCCGCGGAAAAACCGTTGTTGTTGAAGGTTGTGAAGAGAAAATGTCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGGGTGTGAATGCCGCCGa                                                                                                                                                      <  2:1491124/150‑1 (MQ=255)
   tGGCACGATCCTGTCGATTGATACTGAAGGGCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGGGAAGAGAAACTGGCACCACTGGACCTGCTTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATg                                                                                                                                                    <  2:1784780/150‑1 (MQ=255)
      cAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcg                                                                                                                                                 >  1:518731/1‑150 (MQ=255)
      cAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATggcg                                                                                                                                                 >  2:1421730/1‑150 (MQ=255)
         gATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCa                                                                                                                                              >  1:1432612/1‑150 (MQ=255)
            ccTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAGACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCCCTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTg                                                                                                                                           >  2:281347/1‑150 (MQ=255)
            ccTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttc                                                                                                                                           <  2:1843259/150‑2 (MQ=255)
            ccTGTCGATTGATACTGAAGCTCTGACCGCTGACCAAATCCGCCGCGCAAAAACCGTAGTTGTGGAAGGCTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttc                                                                                                                                           <  1:1240650/150‑2 (MQ=255)
               gTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCATCGGCGGCAttcaca                                                                                                                                        <  1:1782730/150‑5 (MQ=255)
               gTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGAAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtt                                                                                                                                        >  1:2001238/1‑150 (MQ=255)
                ccGATGGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCAGGAACTACGCGTGTGAATGCCGCCGATGGCGGCATTGCttt                                                                                                                                       <  1:466929/149‑1 (MQ=255)
                 cGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCtttt                                                                                                                                      <  1:1171365/150‑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)
                               ggTCTGACCGCTGAGCAAATCCGCCGCGGTAAACCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTCGCACCACTGGCCCTGATTCGTCTCGGCAGGAACTAAGCGTGTGAATGCCGCCGCTGGCGGCAGTGCTTTGTTACTTCACGGAat                                                                                                                        >  1:316738/1‑150 (MQ=255)
                               ggTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGCCCTGATTCGTCGCGGCAGGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACTGCACGGAat                                                                                                                        >  1:1263710/1‑150 (MQ=255)
                                 tCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAatat                                                                                                                      <  1:169294/150‑1 (MQ=255)
                                 tCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAatat                                                                                                                      >  1:1544129/1‑150 (MQ=255)
                                     cccGCTGAGCAAATCCGCAGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATGCGTCTCGGCATGAACTAAGCGGGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCCCGGAATATTTTg                                                                                                                  <  2:477622/149‑1 (MQ=255)
                                     aCCGCTGATCAAATCCGCCGCGGAAAACCCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTg                                                                                                                  <  1:1001636/150‑1 (MQ=255)
                                      ccGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAAGGCGGCATAGCTTTTTTACTTCACGGAATATTTTGc                                                                                                                 >  2:1088110/1‑150 (MQ=255)
                                       cgcgGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGcc                                                                                                                <  2:1465975/150‑1 (MQ=255)
                                       cGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAACCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGAGGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGcc                                                                                                                >  2:279949/1‑150 (MQ=255)
                                        gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                               >  2:990948/1‑150 (MQ=255)
                                        gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                               >  1:178467/1‑150 (MQ=255)
                                        gCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCa                                                                                                               <  1:1975084/150‑1 (MQ=255)
                                            aGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCGGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGACGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGt                                                                                                           >  1:326093/1‑150 (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)
                                                    cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCTTTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                   >  2:48609/1‑150 (MQ=255)
                                                    cgccgcGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTcg                                                                                                   >  2:947027/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)
                                                       cgcgGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTTCTTCACGGAATATTTTGCCACGGTCGCTTTCgcgc                                                                                                >  1:590822/1‑150 (MQ=255)
                                                         cgGAAAAACCGTAGTTGTTGCAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCa                                                                                              <  1:639821/150‑1 (MQ=255)
                                                         cgGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCa                                                                                              <  1:1386453/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)
                                                               aaCCGTAGTTGTTGACTGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGTCATGAACTATGCGGGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCt                                                                                        <  2:484498/150‑1 (MQ=255)
                                                                       ttgttgAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATaaata                                                                                >  2:1278091/1‑148 (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)
                                                                                     taaGAGAAACTGGCAACAGTGGACCTGATTCGTCTAGTCATGAACTAAGCGTGTGAATGCCGCCGATGTAGGCATTGCTTTTTTACTTCACGTAATATTTTGACACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCg                                                                  <  1:559688/149‑1 (MQ=255)
                                                                                     gAAGAGAAGCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCTTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCg                                                                  <  2:463027/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)
                                                                                         agaAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCAc                                                              <  2:1407465/150‑1 (MQ=255)
                                                                                            aaCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCActct                                                           >  1:273641/1‑150 (MQ=255)
                                                                                            aaCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGACACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCActct                                                           >  1:1588937/1‑150 (MQ=255)
                                                                                             aCTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCtt                                                          >  2:1989462/1‑150 (MQ=255)
                                                                                              cTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTg                                                         <  1:1629379/150‑1 (MQ=255)
                                                                                                     cacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaa                                                  >  2:593253/1‑150 (MQ=255)
                                                                                                     cacTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaa                                                  <  2:818247/150‑1 (MQ=255)
                                                                                                      acTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaac                                                 <  2:1700218/150‑1 (MQ=255)
                                                                                                        tGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaa                                               <  1:832573/150‑1 (MQ=255)
                                                                                                         ggACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCGTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTaaacaaa                                              >  2:1041568/1‑150 (MQ=255)
                                                                                                           aCCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCCCGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACt                                            >  2:1659793/1‑150 (MQ=255)
                                                                                                              tGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTAtt                                         >  2:2027330/1‑150 (MQ=255)
                                                                                                              tGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTAtt                                         >  1:1915585/1‑150 (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:242723/1‑150 (MQ=255)
                                                                                                                    gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTgg                                   >  1:936550/1‑150 (MQ=255)
                                                                                                                    gTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACGCTTGCAGTAAACAAACTATTGTCTgg                                   <  1:33579/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)
                                                                                                                        cGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCaaa                               >  2:1631374/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)
                                                                                                                              gAACTAAGCGTGTGAACGCCGCCGATGGCGGCATTGCTTTTTTACTTCACTGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAAtcatca                         <  2:371678/150‑1 (MQ=255)
                                                                                                                                 cTAAGCGTGTGAGTGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCa                      <  1:4838/150‑1 (MQ=255)
                                                                                                                                    aGCGTGTGAATGCCGGCGATGGCGGCATTGCTTTTTTACTTCACGGATTATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTACGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAg                   <  2:1381748/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)
                                                                                                                                          tgAAAGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGcc             <  2:1008780/150‑1 (MQ=255)
                                                                                                                                           gAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTGTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCt            <  2:326093/150‑1 (MQ=255)
                                                                                                                                           gAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGGTTCCGTCACTATTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCt            >  1:1395875/1‑150 (MQ=255)
                                                                                                                                                   ccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTTTATGGCAAATCATCACAAAAGAACGACGTAATCGcc    >  1:1904958/1‑150 (MQ=255)
                                                                                                                                                   ccgATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGcc    <  1:850992/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACTGGCAAGATCCTGTCGATTGATACTGAAGGTCTGACCGCTGAGCAAATCCGCCGCGGAAAAACCGTAGTTGTTGAAGGTTGTGAAGAGAAACTGGCACCACTGGACCTGATTCGTCTCGGCATGAACTAAGCGTGTGAATGCCGCCGATGGCGGCATTGCTTTTTTACTTCACGGAATATTTTGCCACGGTCGCTTTCGCGCCATGCGCTAATAAAGACAAGTACGTTTCCGTCACTCTTGCAGTAAACAAACTATTGTCTGGCAAATCATCACCAAAGATCGCCTTAATCGCCAG  >  NC_000913/1629286‑1629584

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: