Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,629,921 C→T E228K (GAG→AAG)  yhiI ← putative membrane fusion protein YhiI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,9210CT98.6% 255.4 / NA 72E228K (GAG→AAG) yhiIputative membrane fusion protein YhiI
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (38/33);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATA  >  NC_000913/3629684‑3630162
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
gTTTCAGCCGTTAAGCGCTGGTTTCGACGGTTTTTGGCGTGAACGGGGCGACCCTGGCGACAAAACGGATGGTTGCAGAAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATAc                                                                                                                                                                                                                                        <  2:333051/251‑1 (MQ=255)
      gCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCAcg                                                                                                                                                                                                                                  >  1:323203/1‑251 (MQ=255)
            cATCGCTGGTTTCGACGGTTTTTGTCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGTCGCGGCATCGAGGATCAGCCGGGCTTCCCCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCCCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGcgtcc                                                                                                                                                                                                                             <  1:313640/250‑1 (MQ=255)
            cATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCAcgcg                                                                                                                                                                                                                                >  2:41224/1‑247 (MQ=255)
            cATCGCTGGTTTCGACGGTTTTTGGCGGGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCAcgcg                                                                                                                                                                                                                                <  1:41224/247‑1 (MQ=255)
            aaTCGCTGGTTTCGACGGTTTTCGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATGCGCAGATCTGGCGCGGCATCGAGTATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATAGAGACGTCGCTGGGATCGGCCATATTCAGCACCCGACCGCCTGCGGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGcgtccg                                                                                                                                                                                                                            <  1:258955/250‑1 (MQ=255)
             aTCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGcgtccgt                                                                                                                                                                                                                           >  1:386661/1‑251 (MQ=255)
                gCTGGGTTCGACGGTTTTTGGCGTGAACTGGCCGACACGGGCGCCAAAACTGATGGTTGCAGGAATCCGCAGATCGGGCGCGCCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCATCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCAc                                                                                                                                                                                                                        <  1:326643/251‑1 (MQ=255)
                   ggTTTCGACGGTTTTTGGGGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTGCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCTATACTGTACGCGTCCGTCAcg                                                                                                                                                                                                                       <  1:306894/249‑1 (MQ=255)
                                         aacggggCGACGCTGGCGACAAAACTGAGGGTTGCTGGAATACGCAGGTCTGGCGCGGCATCGGGGAACTGCCGGGCTTCACCGCCCAGTTGCAGCGTCCCCGTCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCGGAGATCGTCCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGTGATCCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATc                                                                                                                                                                                               <  1:436535/247‑1 (MQ=255)
                                          aCTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCg                                                                                                                                                                                              >  2:389788/1‑251 (MQ=255)
                                           cTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACAGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCg                                                                                                                                                                                              >  2:385601/1‑250 (MQ=255)
                                                cGACACGGGCGAAAAAACGGATGGTTGCAGGAATCCGCAGTTCTGGCTGGGCAGCGGGGATCAGCCGGGCTTCCCCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATAGAGACGTCGCTGAGATCGACCATATCCAGCACCCGACCGCCTTCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCTGGGCTTTCAGTTCGCTGTCAGCGATATCt                                                                                                                                                                                        <  2:231666/251‑1 (MQ=255)
                                                    ccgggCGACAAAACTGATGTGTGCAGGAATCCGCAGTTCGGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCTTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCCCCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCgg                                                                                                                                                                                    <  2:358701/248‑1 (MQ=255)
                                                        gCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAAt                                                                                                                                                                                >  1:382190/1‑251 (MQ=255)
                                                             aaaaCTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGcc                                                                                                                                                                           >  1:29914/1‑251 (MQ=255)
                                                              aaaCTGATGGTTGCAGGAATACGCAGATCGGGCGCGGCATCGAGGATCAGCCGGGCTTCCCCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCg                                                                                                                                                                          <  2:444268/251‑1 (MQ=255)
                                                                   gATGGTTGCAGGAATACGCAGATCTGGGGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCa                                                                                                                                                                      <  1:262730/250‑1 (MQ=255)
                                                                            aGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTtgtg                                                                                                                                                            >  1:309650/1‑251 (MQ=255)
                                                                            aGGAATACGCAGATCTGGCGCGGCATCGAGGATAAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTtgtg                                                                                                                                                            >  1:148932/1‑251 (MQ=255)
                                                                               aattcGCAGAGCGGGCGCGGCATGGAGGATCAGCCGGGCTTCACCTCCAAGGTGAAGCGTGCCCGCCTTTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCg                                                                                                                                                         <  2:238625/247‑1 (MQ=255)
                                                                                aTACGCAGATCGGTCGCGGCTTGGAGTATCAGCCGGTCTTCCCCGCCCATTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGc                                                                                                                                                        <  1:237788/251‑1 (MQ=255)
                                                                                  aCGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCACCTGGCTTGGCAACCCGATACTGCACGCCTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCACTGCGCCGTTCAGTGGCTTGGGCCGCtt                                                                                                                                                      >  2:346989/1‑251 (MQ=255)
                                                                                          cgggCGCGGCATCGAGGATCAGCCGGCCTTCCCCCCCCAGTTTCAGCGGCCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTGCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCgg                                                                                                                                              <  2:170308/249‑1 (MQ=255)
                                                                                            ggcgcggccacgtggacgagcgggccTTCCCCGCCCAGGTCCGGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCACCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGt                                                                                                                                            <  2:8392/226‑1 (MQ=255)
                                                                                            ggCGGGGCATCGAGGATCAGCCGGGCTTCCCCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCgg                                                                                                                                              <  1:20544/249‑1 (MQ=255)
                                                                                              cgcgGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGcgtc                                                                                                                                                                                                                              >  2:359716/1‑167 (MQ=255)
                                                                                              cgcgGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGttt                                                                                                                                          <  1:303266/251‑1 (MQ=255)
                                                                                              cgcgGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTGCGCCTGGCTTGGCAACCCGATACTGCACGcgtc                                                                                                                                                                                                                              <  1:359716/167‑1 (MQ=255)
                                                                                                cgGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTgc                                                                                                                                        <  1:122477/251‑1 (MQ=255)
                                                                                                 ggCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTgcg                                                                                                                                       >  1:215259/1‑251 (MQ=255)
                                                                                                      cGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTGCGCCTg                                                                                                                                   >  1:327408/1‑250 (MQ=255)
                                                                                                       gAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGa                                                                                                                                  >  1:11001/1‑250 (MQ=255)
                                                                                                          gATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGGGCCGCTTCGCCGCGGGTTTGTGCCTGAATg                                                                                                                               >  2:51784/1‑250 (MQ=255)
                                                                                                           gcaagccgggattacccGCCCATTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCAGATAGACGTCGCGGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑at                                                                                                                            <  2:136637/236‑1 (MQ=255)
                                                                                                             cAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCTACGCGGGTTTGAGCCTGAATG‑gta‑t                                                                                                                         >  2:290290/1‑251 (MQ=255)
                                                                                                                  gggCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTGCGCCTGAATG‑GTT‑TTGgtg                                                                                                                    >  2:123835/1‑251 (MQ=255)
                                                                                                                              ccAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTTCGTGCCGCTTCt                                                                                                        <  1:339102/251‑1 (MQ=255)
                                                                                                                               cAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCta                                                                                                       <  1:25339/251‑1 (MQ=255)
                                                                                                                                  tttCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGCATG‑ATA‑TTGGTGCGTGCCGCTTCtata                                                                                                     >  2:200389/1‑250 (MQ=255)
                                                                                                                                    gcAGCGGGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGc                                                                                                   <  1:194007/249‑1 (MQ=255)
                                                                                                                                                ccTGTTCGGTTGGCAGGAAGAAAGTCATATATACTTCGCTGAGATCGACCATATTCAGCACCCGCCCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTTCCGCTTCTATAGCCGCTTTAGAAGcc                                                                                      <  2:382190/251‑1 (MQ=255)
                                                                                                                                                  tGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCg                                                                                     >  1:403824/1‑250 (MQ=255)
                                                                                                                                                   cgacctttGGCAGGAAGAAAGTCATATAGACGTCGCTGAGACCGACCATATTCACCCCCCGACCGCCTGCCGCCAGCACTCCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGAt                                                                                   <  2:102936/245‑1 (MQ=255)
                                                                                                                                                         ttGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGAt                                                                                   >  1:102936/1‑245 (MQ=255)
                                                                                                                                                             cAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCAAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTACCGCGAGGCTTTCAGTTCGACGTCATCAATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCCTCGACGCGGGCTTGCGCCGGAAAGCGTT‑TTGG‑GAGTGCCGCTTCTATATGCGATTTAGAAggagctacccgagctt                                                                         >  2:369368/1‑236 (MQ=255)
                                                                                                                                                              aGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGGCGCTTCTATAGCCGCCTTAGAAGCCGGTACCTGAGCttt                                                                        >  1:10400/1‑251 (MQ=255)
                                                                                                                                                              aGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACGGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTTTGCCGCTTCGACGCGGGGTTGCGCCTGAATG‑AGAGTTGG‑GCGTGCCGCTTCTATAGCCACTTTATAGGCCGATACCTGAGCttt                                                                        >  2:318144/1‑251 (MQ=255)
                                                                                                                                                                  agaaAGTCATATAGACGTCGCTGAGATCGCCCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCCCTTCGCCGCGGGTTTGCGCCTGCATG‑ATA‑TTGGGGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGcc                                                                    >  2:295344/1‑251 (MQ=255)
                                                                                                                                                                   gaaAGTCATACAGACGTCGCTGAGATCGAACATATTCAGCACCAGACCGCCTGCCGCCAGCCCTTCGCCTGGCTTAGCAACCCGATACTCCACGCGTCCGTCACGCGGGGCTTTCCGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTCTCCCGCCTAGACGCGGGTTTGCGCCTGAAAG‑ATA‑TTGGTGCGGCCCGCTTCTCTAGCCGCTTTAGAACCCGACACATCAGCCTTCGCCg                                                                   >  2:184836/1‑251 (MQ=255)
                                                                                                                                                                     aaGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTTGGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑GTA‑TTGGTGCGGGCCGCTTTTATTACCGCTTTAGAAGCCGATACCTGAGCTTTCGcccg                                                                  >  2:361421/1‑248 (MQ=255)
                                                                                                                                                                     aaGTCATATAGACGTCGCTGAGACCGACCTGATTCAGCCCCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGGGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGAt                                                                 <  1:200389/251‑1 (MQ=255)
                                                                                                                                                                         cATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTcc                                                              <  1:416817/250‑1 (MQ=255)
                                                                                                                                                                              agtcGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCa                                                        <  1:214880/248‑1 (MQ=255)
                                                                                                                                                                                aCGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGc                                                      >  1:161245/1‑251 (MQ=255)
                                                                                                                                                                                    cGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCACTGCGCCGTTCAGTGGCTTGTGCCGCTTCGCCGCGGGTTTGCGACTGAATG‑CTA‑TTAGTGCGTGCCGCTTCTATAGCCCCTTTTGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTccg                                                  >  1:441674/1‑249 (MQ=255)
                                                                                                                                                                                       tGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCCCTTCGACGCGGGTTTTCGCCTGAATG‑ATA‑TTTGTGCGTGCCGCTTCTATAGCCGCTTTAGAATCCGTTACCTTCACTTTCGCCGTTTCCAGCGCAGCTCGGgcc                                               >  1:39675/1‑250 (MQ=255)
                                                                                                                                                                                         agaTCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGGTACCTTAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGct                                             >  1:33311/1‑251 (MQ=255)
                                                                                                                                                                                              gACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCtctc                                          >  2:113431/1‑249 (MQ=255)
                                                                                                                                                                                               cccATATTCAGCCCCCGCCCCCCTGCCGCCAGCCCTTCGCCGGGCTTGGCAACCCGATCCTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAgc                                       <  2:105007/250‑1 (MQ=255)
                                                                                                                                                                                                   aattcaggccccGACCGCCTTCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCGCGTTTCCGTCACGCGGGGCTTTCAGTTCTCTTTCATCGATATCTGCGGCAATGTGCCGTTCTGTGGCTTGTGTCGCTTCGACGCGGGTTTGCGCTTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGAGACCGGAGCTTTCGCCGATCCCAGCGCAGCTCGGGCGCTCTCAgcggcg                                   <  1:266142/242‑1 (MQ=255)
                                                                                                                                                                                                    atTCAGCCCCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAgcggcg                                   <  1:339315/250‑1 (MQ=255)
                                                                                                                                                                                                               gACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTTGTGCGTGCCGCCTCTATAGCCGCTTTTGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCCGCTCGGGCGCTCTCCGCGGCGGCGCGATCggcc                       >  2:388648/1‑248 (MQ=255)
                                                                                                                                                                                                               gACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGgtgc                                                                                                                   >  2:371695/1‑159 (MQ=255)
                                                                                                                                                                                                                 ccGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATT‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGCTACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATc                     >  2:308641/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGCCGCGGGTTTGCGCCTGAATG‑GTA‑TTGGTTCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGGGCTTTCTCCGGTTCCCGCGCAGCTCGGGCGCTCTCACCGGCGGCGCGGTCGTCATcccg                  >  2:176020/1‑249 (MQ=255)
                                                                                                                                                                                                                    ccTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCCTCTATAACCGCTTTAGAAGCCGGTACCTGGGCTTTCGCCGATTCCAGCGCAGCACGGGCGCTCTCAGCGGCGGCGCGATCGtcctcccc                  >  1:144956/1‑249 (MQ=255)
                                                                                                                                                                                                                       gccgccAGAACTTCGCCGGGCTTGGTAACCCGATACTGCACGCGTCCGTCGCGCGGGGTTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCGTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATTCCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCt                <  2:130118/250‑1 (MQ=255)
                                                                                                                                                                                                                          gccAGCACTTCGCCTGGCTTGGCAACCCGATCCTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGGGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTg            >  1:94016/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTgcgc         <  2:11001/251‑1 (MQ=255)
                                                                                                                                                                                                                                aCTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTgcgc         >  1:154047/1‑248 (MQ=255)
                                                                                                                                                                                                                                     gCCTGGCTCGTCAACCCGATACTGCACCCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGCTTGCGCCTGGATG‑GTA‑TTGGTTCGGGCCGCTTCTATAGCCGCTTTAGAAACCGATAACTGCGACTTCGCCCGTTCCAGCGCCGCTCCGGCTCTCTCCGCCGCGGCGCGATCGTCCTCCCGCTGTTgcgccagaata  >  2:3997/1‑243 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
GTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATG‑ATA‑TTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATA  >  NC_000913/3629684‑3630162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑