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

GTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATA  >  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)
                                                                                                           gcaagccgggattacccGCCCATTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCAGATAGACGTCGCGGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGat                                                                                                                            <  2:136637/236‑1 (MQ=255)
                                                                                                             cAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCTACGCGGGTTTGAGCCTGAATGgtat                                                                                                                         >  2:290290/1‑251 (MQ=255)
                                                                                                                  gggCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTGCGCCTGAATGGTTTTGgtg                                                                                                                    >  2:123835/1‑251 (MQ=255)
                                                                                                                              ccAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTTCGTGCCGCTTCt                                                                                                        <  1:339102/251‑1 (MQ=255)
                                                                                                                               cAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCta                                                                                                       <  1:25339/251‑1 (MQ=255)
                                                                                                                                  tttCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGCATGATATTGGTGCGTGCCGCTTCtata                                                                                                     >  2:200389/1‑250 (MQ=255)
                                                                                                                                    gcAGCGGGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGc                                                                                                   <  1:194007/249‑1 (MQ=255)
                                                                                                                                                ccTGTTCGGTTGGCAGGAAGAAAGTCATATATACTTCGCTGAGATCGACCATATTCAGCACCCGCCCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTTCCGCTTCTATAGCCGCTTTAGAAGcc                                                                                      <  2:382190/251‑1 (MQ=255)
                                                                                                                                                  tGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCg                                                                                     >  1:403824/1‑250 (MQ=255)
                                                                                                                                                   cgacctttGGCAGGAAGAAAGTCATATAGACGTCGCTGAGACCGACCATATTCACCCCCCGACCGCCTGCCGCCAGCACTCCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGAt                                                                                   <  2:102936/245‑1 (MQ=255)
                                                                                                                                                         ttGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGAt                                                                                   >  1:102936/1‑245 (MQ=255)
                                                                                                                                                             cAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCAAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTACCGCGAGGCTTTCAGTTCGACGTCATCAATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCCTCGACGCGGGCTTGCGCCGGAAAGCGTTTTGGGAGTGCCGCTTCTATATGCGATTTAGAAggagctacccgagctt                                                                         >  2:369368/1‑236 (MQ=255)
                                                                                                                                                              aGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGGCGCTTCTATAGCCGCCTTAGAAGCCGGTACCTGAGCttt                                                                        >  1:10400/1‑251 (MQ=255)
                                                                                                                                                              aGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACGGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTTTGCCGCTTCGACGCGGGGTTGCGCCTGAATGAGAGTTGGGCGTGCCGCTTCTATAGCCACTTTATAGGCCGATACCTGAGCttt                                                                        >  2:318144/1‑251 (MQ=255)
                                                                                                                                                                  agaaAGTCATATAGACGTCGCTGAGATCGCCCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCCCTTCGCCGCGGGTTTGCGCCTGCATGATATTGGGGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGcc                                                                    >  2:295344/1‑251 (MQ=255)
                                                                                                                                                                   gaaAGTCATACAGACGTCGCTGAGATCGAACATATTCAGCACCAGACCGCCTGCCGCCAGCCCTTCGCCTGGCTTAGCAACCCGATACTCCACGCGTCCGTCACGCGGGGCTTTCCGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTCTCCCGCCTAGACGCGGGTTTGCGCCTGAAAGATATTGGTGCGGCCCGCTTCTCTAGCCGCTTTAGAACCCGACACATCAGCCTTCGCCg                                                                   >  2:184836/1‑251 (MQ=255)
                                                                                                                                                                     aaGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTTGGCCGCTTCGACGCGGGTTTGCGCCTGAATGGTATTGGTGCGGGCCGCTTTTATTACCGCTTTAGAAGCCGATACCTGAGCTTTCGcccg                                                                  >  2:361421/1‑248 (MQ=255)
                                                                                                                                                                     aaGTCATATAGACGTCGCTGAGACCGACCTGATTCAGCCCCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGGGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGAt                                                                 <  1:200389/251‑1 (MQ=255)
                                                                                                                                                                         cATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTcc                                                              <  1:416817/250‑1 (MQ=255)
                                                                                                                                                                              agtcGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCa                                                        <  1:214880/248‑1 (MQ=255)
                                                                                                                                                                                aCGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGc                                                      >  1:161245/1‑251 (MQ=255)
                                                                                                                                                                                    cGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCACTGCGCCGTTCAGTGGCTTGTGCCGCTTCGCCGCGGGTTTGCGACTGAATGCTATTAGTGCGTGCCGCTTCTATAGCCCCTTTTGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTccg                                                  >  1:441674/1‑249 (MQ=255)
                                                                                                                                                                                       tGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCCCTTCGACGCGGGTTTTCGCCTGAATGATATTTGTGCGTGCCGCTTCTATAGCCGCTTTAGAATCCGTTACCTTCACTTTCGCCGTTTCCAGCGCAGCTCGGgcc                                               >  1:39675/1‑250 (MQ=255)
                                                                                                                                                                                         agaTCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGGTACCTTAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGct                                             >  1:33311/1‑251 (MQ=255)
                                                                                                                                                                                              gACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCtctc                                          >  2:113431/1‑249 (MQ=255)
                                                                                                                                                                                               cccATATTCAGCCCCCGCCCCCCTGCCGCCAGCCCTTCGCCGGGCTTGGCAACCCGATCCTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAgc                                       <  2:105007/250‑1 (MQ=255)
                                                                                                                                                                                                   aattcaggccccGACCGCCTTCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCGCGTTTCCGTCACGCGGGGCTTTCAGTTCTCTTTCATCGATATCTGCGGCAATGTGCCGTTCTGTGGCTTGTGTCGCTTCGACGCGGGTTTGCGCTTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGAGACCGGAGCTTTCGCCGATCCCAGCGCAGCTCGGGCGCTCTCAgcggcg                                   <  1:266142/242‑1 (MQ=255)
                                                                                                                                                                                                    atTCAGCCCCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAgcggcg                                   <  1:339315/250‑1 (MQ=255)
                                                                                                                                                                                                               gACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTTGTGCGTGCCGCCTCTATAGCCGCTTTTGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCCGCTCGGGCGCTCTCCGCGGCGGCGCGATCggcc                       >  2:388648/1‑248 (MQ=255)
                                                                                                                                                                                                               gACCGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGgtgc                                                                                                                   >  2:371695/1‑159 (MQ=255)
                                                                                                                                                                                                                 ccGCCTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATTATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGCTACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATc                     >  2:308641/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGCCGCGGGTTTGCGCCTGAATGGTATTGGTTCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGGGCTTTCTCCGGTTCCCGCGCAGCTCGGGCGCTCTCACCGGCGGCGCGGTCGTCATcccg                  >  2:176020/1‑249 (MQ=255)
                                                                                                                                                                                                                    ccTGCCGCCAGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTGCGCCTGAATGATATTGGTGCGTGCCGCCTCTATAACCGCTTTAGAAGCCGGTACCTGGGCTTTCGCCGATTCCAGCGCAGCACGGGCGCTCTCAGCGGCGGCGCGATCGtcctcccc                  >  1:144956/1‑249 (MQ=255)
                                                                                                                                                                                                                       gccgccAGAACTTCGCCGGGCTTGGTAACCCGATACTGCACGCGTCCGTCGCGCGGGGTTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCGTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATTCCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCt                <  2:130118/250‑1 (MQ=255)
                                                                                                                                                                                                                          gccAGCACTTCGCCTGGCTTGGCAACCCGATCCTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGGGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTg            >  1:94016/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGCACTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTgcgc         <  2:11001/251‑1 (MQ=255)
                                                                                                                                                                                                                                aCTTCGCCTGGCTTGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTgcgc         >  1:154047/1‑248 (MQ=255)
                                                                                                                                                                                                                                     gCCTGGCTCGTCAACCCGATACTGCACCCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGCTTGCGCCTGGATGGTATTGGTTCGGGCCGCTTCTATAGCCGCTTTAGAAACCGATAACTGCGACTTCGCCCGTTCCAGCGCCGCTCCGGCTCTCTCCGCCGCGGCGCGATCGTCCTCCCGCTGTTgcgccagaata  >  2:3997/1‑243 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
GTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCAACCCGATACTGCACGCGTCCGTCACGCGGGGCTTTCAGTTCGCTGTCATCGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATA  >  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: