Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,146,110 C→T E162E (GAG→GAA yceF ← m(7)GTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,146,1100CT98.5% 193.6 / NA 67E162E (GAG→GAAyceFm(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (31/34);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTT  >  NC_000913/1145864‑1146355
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
aaGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCa                                                                                                                                                                                                                                                    >  1:230369/1‑250 (MQ=255)
        tttGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGGTTAGCGGTAAACCAACCAGCGTGTTAGGATCCCGCCCCTCTCAACGTTCCAACAGCGt                                                                                                                                                                                                                                            >  2:272415/1‑250 (MQ=255)
         ttGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTa                                                                                                                                                                                                                                           >  1:29154/1‑250 (MQ=255)
              aaaGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGcc                                                                                                                                                                                                                                      >  1:674143/1‑250 (MQ=255)
              aaaGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGcc                                                                                                                                                                                                                                      >  1:513579/1‑250 (MQ=255)
                        aaCGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCa                                                                                                                                                                                                                            <  1:552299/250‑1 (MQ=255)
                          cgcgcgCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCActc                                                                                                                                                                                                                         >  1:238072/1‑251 (MQ=255)
                            cgcgagaaTACAACCATAGAAACCTGATACTGGGCGGTTGGCCTGATAGAGCGTTTACGCCGCATCAGGCCGCCAGCCCCTATGGCCGGATGCGACGGAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCGGACGTAACATCGGACACAGCGTGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCGCTAAACGTTCAAACAGCGTAATGCCAAATCCttttctctt                                                                                                                                                                                                                       <  2:239800/245‑8 (MQ=255)
                             gcgcTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTa                                                                                                                                                                                                                      <  1:450404/251‑1 (MQ=255)
                                      cAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTAccc                                                                                                                                                                                                             <  1:58688/251‑1 (MQ=255)
                                      cAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTAccc                                                                                                                                                                                                             >  1:538189/1‑251 (MQ=255)
                                                       tAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCAGCGCAAAAAACTATCCCATCAGCGGGATTTTCCCTTCCCGACGTAAAATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCCAACAGCGTAAAGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGGGCCGGGGAtt                                                                                                                                                                                            >  2:361649/1‑250 (MQ=255)
                                                           tggCGGGTGGCCTGATACGGCGTTTACGCCGCACCAGGCCGACAGCACCGACTGCCGGAGGCGACGTAACCCTATCCGGCGAAAAAAACTACCCCAGCAGCTGGTTTTCCCAATCCCGACGTAACAGCTGACACAGCTCGAGTAGCGGTAAACCAGCCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAAGACTTCACTCTGAAATCTTCCCGCGCAGTGCAGGGGATGctc                                                                                                                                                                                         <  2:219836/249‑1 (MQ=255)
                                                           gggCGTTTGGCCGAATAAGGCGTTTCCGCCGCATCAGGACGCAAGCACTGATTGCCTGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGTTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGctct                                                                                                                                                                                        <  2:242075/251‑1 (MQ=255)
                                                            ggCGGTTGGCCTGAAAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAATCGACCCCATCAGCGGGTTTTTCCCTCCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCtt                                                                                                                                                                                       <  2:674143/251‑1 (MQ=255)
                                                              cgtttGGCCTGATAAGGCGTTGACGCCGCCTCAGGCCGCCAGCATCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCTCGCCTCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTa                                                                                                                                                                                     <  1:655212/248‑1 (MQ=255)
                                                                     cctttaaatgtgTTTCCGCCGCATCAGGCCGCCAGCACCGATTTCCATATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCTACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCGTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTGCGCACGt                                                                                                                                                                              <  2:567155/240‑1 (MQ=255)
                                                                      gatataggggattaccgccgcACCCGGCCGCCAGCACCGATTGCAGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCTATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCCTTACGCACGTa                                                                                                                                                                             <  2:257922/236‑1 (MQ=255)
                                                                           aaGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGACAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTAt                                                                                                                                                                        >  2:642301/1‑251 (MQ=255)
                                                                            aGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCGAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATc                                                                                                                                                                       <  2:244476/251‑1 (MQ=255)
                                                                                 tttACGCCGCATCAGTCAGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGGTTTTCCCTTCCCGACGTAACCACTGACACAGCGCGATTAGCGGAACACCAACCAGCGTGTTAGGATCACGCCCCTATAAACGTTCAAACAACGTAATGCCAATACCATCACTCTTCAAGCTACCCGCGCATAGCAGGGAGTGCTCTTTACGCACGTAATCATAAAtcc                                                                                                                                                                  >  2:103339/1‑250 (MQ=255)
                                                                                    tgcccgaatcaggaggccactatcGATTGCCGGAAGCGACGTAACCACATCCGGCGCAAAAAACTAATCCATCAGCGGGTTTTTCCCTTCCCGTCGTAACAACTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGAGCGCGCCCCCCTAAACGGTCAAACAGGGTAATGCCAAATCCTGCACCCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCATGTAATAATCAATCTcc                                                                                                                                                                <  2:423524/227‑1 (MQ=255)
                                                                                             cAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAACGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCCCGTAATTATCAATCTCCGCCTCGCTc                                                                                                                                                       >  2:41258/1‑250 (MQ=255)
                                                                                             cAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCa                                                                                                                                                      <  1:272415/251‑1 (MQ=255)
                                                                                              agcccgccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCTCTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAg                                                                                                                                                     <  2:228872/248‑1 (MQ=255)
                                                                                                   gccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCa                                                                                                                                                                                                                                     >  1:403763/1‑166 (MQ=255)
                                                                                                   gccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGAc                                                                                                                                                <  1:485738/251‑1 (MQ=255)
                                                                                                                gCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGt                                                                                                                                    <  1:42866/250‑1 (MQ=255)
                                                                                                                 ccGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACCTCTGACACAGAGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATAACGCCCCTCTAAACGTTCAAACAACGTAATGCCAAATCATTCACTCTTAAAGCTACCCGCGAAGTGAAGGGGATGATCTTTACGACCGTAATTATCAATCTCAGCCTAGCTAAGGTGACGGAAATGGACGTc                                                                                                                                   >  2:663868/1‑250 (MQ=255)
                                                                                                                     aTGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAg                                                                                                                              <  1:65028/251‑1 (MQ=255)
                                                                                                                         gACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTc                                                                                                                          <  1:594122/251‑1 (MQ=255)
                                                                                                                                ccACAGCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCt                                                                                                                   <  1:461838/251‑1 (MQ=255)
                                                                                                                                       cGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCCCTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCAc                                                                                                                                                                                >  2:204647/1‑183 (MQ=255)
                                                                                                                                       cGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCAc                                                                                                                                                                                <  1:204647/183‑1 (MQ=255)
                                                                                                                                             aaaaaaCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCCGAGCAGGGGATGCTCTTTACg                                                                                                                                                                                   >  2:365245/1‑174 (MQ=255)
                                                                                                                                             aaaaaaCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACg                                                                                                                                                                                   <  1:365245/174‑1 (MQ=255)
                                                                                                                                                  aCTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTc                                                                                                 <  1:260765/251‑1 (MQ=255)
                                                                                                                                                    tATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATg                                                                                                                                         <  2:187589/209‑1 (MQ=255)
                                                                                                                                                    tATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATg                                                                                                                                         >  1:187589/1‑209 (MQ=255)
                                                                                                                                                     aTCCAATCAGCGGGTTTACCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTTCAGATGCCCATTCGc                                                                                               <  2:230369/250‑1 (MQ=255)
                                                                                                                                                       cccATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGCTTAGCGGTAAACCCACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGGGACGGAAATGGGCGTCAAAAGGCTCCACTTCTGTTTGCCGATGCCCATTCGCCg                                                                                             >  2:212671/1‑250 (MQ=255)
                                                                                                                                                             aGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTCATGCCAAATCCTTCACTATTAAAGCTACCCGCGCAGTGCAGGGGATGATCTTTACGCCCGCAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAACGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTaa                                                                                       >  1:102623/1‑250 (MQ=255)
                                                                                                                                                             aGCGGGTTTTTCCCTTCACGACGTAACATCTGACACAGCGCGCTTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAACCGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCAGCGAAGAGCAGGGGATGCTCTTTACGCACGTAATAATCACTCTCCGCCTCGCCCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTATTTGCAGAAGCCCATTCGCCGACTTaa                                                                                       >  2:327061/1‑250 (MQ=255)
                                                                                                                                                                      ttcccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAg                                                                                                                              <  1:99296/202‑1 (MQ=255)
                                                                                                                                                                      ttcccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGAGCCGGAAATGGACGTCAAAAg                                                                                                                              >  2:99296/1‑202 (MQ=255)
                                                                                                                                                                       tcccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGctct                                                                                                                                                                                        >  2:483969/1‑143 (MQ=255)
                                                                                                                                                                       tcccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGctct                                                                                                                                                                                        <  1:483969/143‑1 (MQ=255)
                                                                                                                                                                         ccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTaa                                                                                                                                                                            <  2:416543/153‑1 (MQ=255)
                                                                                                                                                                         ccttcccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTaa                                                                                                                                                                            >  1:416543/1‑153 (MQ=255)
                                                                                                                                                                            ccccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATGCGCCGAATTAAACAACGccagtccagt                                                                       <  1:273006/250‑1 (MQ=255)
                                                                                                                                                                              ccGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAAGCCTTCACTCGTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTGCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGtat                                                                     <  2:296791/251‑1 (MQ=255)
                                                                                                                                                                               cGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCa                                                                                                            <  1:79865/211‑1 (MQ=255)
                                                                                                                                                                               cGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCa                                                                                                            >  2:79865/1‑211 (MQ=255)
                                                                                                                                                                               cGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGtat                                                                     >  2:138584/1‑250 (MQ=255)
                                                                                                                                                                                gACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACAAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGAGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGtata                                                                    >  2:540028/1‑250 (MQ=255)
                                                                                                                                                                                            gACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGAc                                                           <  1:642301/247‑1 (MQ=255)
                                                                                                                                                                                               acaGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTACACGTTCAAACAGCGTAATGCCAACTCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGTTATTg                                                    >  1:332606/1‑251 (MQ=255)
                                                                                                                                                                                               acaGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTg                                                    >  1:233214/1‑251 (MQ=255)
                                                                                                                                                                                                      gATTAGCGTTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCTTTCGCCGAATTAaaca                                                                                    <  2:492577/212‑1 (MQ=255)
                                                                                                                                                                                                      gATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAaaca                                                                                    >  1:492577/1‑212 (MQ=255)
                                                                                                                                                                                                          aggggtaaccaaaccaGCGGGTTATGTTCAAGCCCCTCTAACCGTTCAAACAGCGTAATGCCAAAACCTTCTCTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACTTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCttt                                         <  2:496124/240‑1 (MQ=255)
                                                                                                                                                                                                          aGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCGTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAAGGGACGTCAAAAGGCTCCACTTCTGTTTGCATATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCtt                                          <  2:29154/250‑1 (MQ=255)
                                                                                                                                                                                                                  accaaccaGCGTGTTAGGATCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGCAGGTGACGATATTGCCGCTGGCTTTGCGTAAtt                                 >  1:341682/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aTCACGCCCCTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTttctt               <  1:212671/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cccTCTAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGttc         >  2:440109/1‑249 (MQ=255)
                                                                                                                                                                                                                                             ctctAAACGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACg      >  1:650803/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aaaCGTTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGGGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGttt  >  2:530179/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
AAGCGTTGTTTGCAAAAGCTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCACTCTTAAAGCTACCCGCGCAGTGCAGGGGATGCTCTTTACGCACGTAATTATCAATCTCCGCCTCGCTCAGGTGACGGAAATGGACGTCAAAAGGCTCCACTTCTGTTTGCAGATGCCCATTCGCCGAATTAAACAACGCCAGTCCAGTATAGAAGGTGACGATATTGCCGCTGGCTTTGCGTAATTGCAGACGGGCATTTTCTTCCGTTAACGGTTT  >  NC_000913/1145864‑1146355

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