Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,585,480 G→A Q428* (CAA→TAA)  hsdR ← type I restriction enzyme EcoKI endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,585,4800GA98.5% 216.9 / ‑5.1 68Q428* (CAA→TAA) hsdRtype I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (29/37);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGT  >  NC_000913/4585234‑4585725
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ccTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGATAATGCGTTTAATACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                                                                                                                                                                                                   <  2:300416/251‑1 (MQ=255)
 cTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGCGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCCGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCCGGGGTTTCGGCAAACCCCTGATTCTGTTTTTACGGtt                                                                                                                                                                                                                                                  >  2:53371/1‑251 (MQ=255)
   gTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGtt                                                                                                                                                                                                                                                  <  2:401854/249‑1 (MQ=255)
       gCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGc                                                                                                                                                                                                                                            <  1:72233/251‑1 (MQ=255)
        ccAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGcc                                                                                                                                                                                                                                           <  1:12522/251‑1 (MQ=255)
                     ggTCGATAAGTAAGCGAATGCGTTTAAAACGCTGGAACTGGATCAGGCGGAACATCATGGCGATTGCAGTACGGGTTTTCCCGGTCCCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGCCGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTAATTCTGTTTTTACGGTTCGCTGCCGAGCCTTTCcagc                                                                                                                                                                                                                              >  2:57847/1‑251 (MQ=255)
                        cGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCagcagc                                                                                                                                                                                                                           >  1:401490/1‑251 (MQ=255)
                         gACAAGGACGAGAATGCGTTTAAAACGCTGGGACTCGATCAGGCGGAACATCATGGCGATTGCAGTACGGGTTTTAACGTAACCCGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCACCCGCGCGGACGGCATCTTCCTGATAATAACTCATGCCCAGCTCGCTCATGCCAGGGGTCTCGGCAAACCACTTCTTCTGTTTTCACGGTTCGCTGCCGAGCCTTTCcagcaacc                                                                                                                                                                                                                          >  2:249680/1‑248 (MQ=255)
                          aCAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                                                                                                                                                                                                         <  2:364681/251‑1 (MQ=255)
                          aCAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGAGTCTGTTTTTACGGTTCGCTGCCGTGCATTTCCAGCAGctc                                                                                                                                                                                                                         <  2:38840/251‑1 (MQ=255)
                              ggAAGAGAATGCGTTTACAACGCTGGGACTGCATCCGGCGGAACATCATGGCGATTGCCGTACAGGTTTTACTGGTAACGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGAAGATTGCCTTTTCAACCGCGCGGACGGCATATTCCTGATAATAACGCAGGACAAGCTCGCTCATGCCACGGGTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTACAGCAGCTCTTcc                                                                                                                                                                                                                     >  2:208951/1‑251 (MQ=255)
                                   agaATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGACGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCg                                                                                                                                                                                                                >  1:285040/1‑251 (MQ=255)
                                        gcggttAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCTATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGcc                                                                                                                                                                                                           <  2:40253/247‑1 (MQ=255)
                                                   gCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAg                                                                                                                                                                                                <  1:318009/251‑1 (MQ=255)
                                                    cTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAgg                                                                                                                                                                                               >  2:147992/1‑251 (MQ=255)
                                                               tCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGGGCCACTCGGGTAAGGCTTTCGACat                                                                                                                                                                                     >  2:323243/1‑250 (MQ=255)
                                                                 aGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTTCCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTAATAATAACGCAGACCCAGCTCGCTCATGACAGGGTTATCCGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTACCAGCAGATCTTCCGGGCGGTGCCACTCGGGTAAGACTTTCGACATAtt                                                                                                                                                                                  >  1:268703/1‑251 (MQ=255)
                                                                      gAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCggg                                                                                                                                                                             >  1:408479/1‑251 (MQ=255)
                                                                               ggCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATcacgca                                                                                                                                                                    <  1:113651/251‑1 (MQ=255)
                                                                                    ttGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTcg                                                                                                                                                               >  1:183216/1‑251 (MQ=255)
                                                                                               gttgtaccggtaccggTCGCCATCGCCAGCAGGATCCCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCTCGTCGCGATTCCAGAt                                                                                                                                                    <  2:41252/247‑1 (MQ=255)
                                                                                                ttttACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATg                                                                                                                                                   >  1:268385/1‑251 (MQ=255)
                                                                                                  ttACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGcc                                                                                                                                                 >  1:39759/1‑251 (MQ=255)
                                                                                                        gtaccggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTGCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCtttt                                                                                                                                           <  1:248399/251‑1 (MQ=255)
                                                                                                             ggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCt                                                                                                                                      >  2:147141/1‑251 (MQ=255)
                                                                                                                      cgcccGCAGGATCTCTTGTTTCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCGTTCGACAGATTGCGTGTATCACGCACGTCGCGATACCAGATGTCGCTTTTGGTCTTCAttgttg                                                                                                                             <  2:285040/251‑1 (MQ=255)
                                                                                                                       gccaGCGGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATCATAACGCAGGCCCAGCTCGCTCATGCCAGGGTAATCGGCAAGACACTGATTCTGTTTTTACGGTTCGCGGCCGGGCATTTACAGCAGCTCTTCCTGGCGGTGCCACTCGGGTAAGTCTTTCGACATATTGCGGGTATCACGCACGTCGCGATCCCAGCTGCCGCTTTTGGTCTTCAttgttgt                                                                                                                            <  2:227296/251‑2 (MQ=255)
                                                                                                                         cagcagGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTgcgc                                                                                                                          >  1:404274/1‑251 (MQ=255)
                                                                                                                            caggctatctTGTTGCCCCTTGACGATTGCCTTTTCCACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGATCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGACACACGGGTAAAGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCCGATGCCGCTTTTGGTCTTCATAGTTGCGCggc                                                                                                                       >  2:278612/8‑250 (MQ=255)
                                                                                                                             agGATCTCTTGTTGCCCCTTGACGATTGCCTGTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGGTTTTACGGTGCGCGGCCGAGCAGTTCCAGCAGCTCTTCCGGGCGGTGCCGCTCGGGTAAGGCGTTCGACATATTGCGGGTATGACGCACGTCGCGATTCCAGATGCCGCTGTTGGTCTTCATTGTTGCGCGGTa                                                                                                                      <  1:362094/251‑1 (MQ=255)
                                                                                                                                    gttgttgCGCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCAGTGATTCTGTTTTTACGGTGCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGGCATATTGCGGGTATCACTCACGTCGGGATGCCAGATGCGGCTTTTGGTCTTCATTGTTGCGCGGTATTcg                                                                                                                  <  1:249680/247‑1 (MQ=255)
                                                                                                                                    ctTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACTCAGTCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAg                                                                                                                                                                                                                                        <  2:283419/130‑1 (MQ=255)
                                                                                                                                    ctTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAg                                                                                                                                                                                                                                        >  1:283419/1‑130 (MQ=255)
                                                                                                                                    ctTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTc                                                                                                                   <  1:87397/247‑1 (MQ=255)
                                                                                                                                      tgttgCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTAATCgcgc                                                                                                               <  1:116245/249‑1 (MQ=255)
                                                                                                                                        ttgCCCCTTGACGATTGACTTTTCAACCGAGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATACCAGGGTTATAGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGGGCCACTCGGGTAAGGCTTTTGACATATTGCGGGTATCAAGCACGACGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCTGTATTCGCGCCCGt                                                                                                           >  2:95396/1‑251 (MQ=255)
                                                                                                                                            cccTTGACGATTGCCTTTTCAACCCCGCGTACGGCATCTTCCTGATAATAACGCAGGGCCAGCTCGCTCATTCCAGGTTTAGCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGTGGTGCCACTCGGGTAATGCTTTCGACATATTGCGGGTATCACGCACGGCGCGATGCCAGATGTCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGt                                                                                                       <  1:165937/251‑1 (MQ=255)
                                                                                                                                             ccTTGACGATTGCCTTTTCAACCTCGCGGACTGCATCTTCCGGATAATAACGCAGGCCCAGCTCGCTCATGCCAGTGTTATCGGCAAACCACTGATTCTGTTTTTTCGGTTCGCTGCCGAGCATGTCCAGCAGCTCTTCCGTGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATCCCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCTCCCGTTGGTc                                                                                                      <  2:78745/251‑1 (MQ=255)
                                                                                                                                               ttGACGATTGACTTTTCAACCGCGCGGCCAGCATATTCCTGATCATAACACCGGCCCTCATCGATCATGCCAGGGTTATCGGAAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGACACTCGGGAAAGGCTTTCGACATATTGCGGGTATCACGAACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGATGCGCGGTATTAGCGCCCGTTGGTCg                                                                                                     >  2:51058/1‑250 (MQ=255)
                                                                                                                                                   cGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAg                                                                                                <  1:95396/251‑1 (MQ=255)
                                                                                                                                                    gATTTCCTTTTCAACCGCGCGGACGGCGTCTTACTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGAGTTATCGGAAAAACACTGATTCTGTTTTTACGGTTCACTGCCGAGCCTTTCCAGCAGCACTTCCGGGCGGTACCACTCGGGTAAGGCTTTCGACATCTCACGGATATCACGCACGTCGCGATCACAGCTGCCTCTTTTGGTCTTAATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGc                                                                                               >  2:24033/1‑251 (MQ=255)
                                                                                                                                                        gCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGa                                                                                            <  1:276817/250‑1 (MQ=255)
                                                                                                                                                             ttCAACCGCGCGGACGGCTTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAGTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgggg                                                                                       <  2:404274/250‑1 (MQ=255)
                                                                                                                                                              tCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGAAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGAGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGAt                                                                                     >  2:20958/1‑251 (MQ=255)
                                                                                                                                                                aaCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGAttt                                                                                   <  1:98818/251‑1 (MQ=255)
                                                                                                                                                                  ccGCGCAGACGGCATCTTCCTGATAATAACGCAGACCCAGCTCCCTCATGCCAGGATTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGAAAAGCTTTCGACATATTGCGGGTATCACGCACGTCGCGCTACCAGATGCCGCTTTTGGTCTTCATTGTGGCGCGGTACTTGCTCCCGTTGGTCGAGTAGCAGAAGGGGAtttt                                                                                  >  2:378778/1‑250 (MQ=255)
                                                                                                                                                                       cgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCg                                                                            <  1:60861/251‑1 (MQ=255)
                                                                                                                                                                       cgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCg                                                                            >  1:419138/1‑251 (MQ=255)
                                                                                                                                                                       cagacggtcacTTCCTGATAATAACGCAGGCTCAGCCCGCTCATGCTAGGGTTATCGGCAAACCTCTGATTCTCTTTTTACTGTTCGCTGCCGAGCATTTCAAGCAGTTCTTCCGGGCGGTGCCAATCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCCCGTCGCGATACCAGATGACGCTTTTGGTCTTCATTGTTGCGCGTGATTCGCGCCTGTTGGTTGAGTAGCAGAAGGGGATTTTAAAcag                                                                            <  1:60860/241‑3 (MQ=255)
                                                                                                                                                                            ggCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGAATTTCCAGCAGCTCTTCCGGGCGGt                                                                                                                                                                                                              >  1:37890/1‑116 (MQ=255)
                                                                                                                                                                                    ccTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCt                                                                <  1:53371/250‑1 (MQ=255)
                                                                                                                                                                                       gATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGgtg                                                            <  1:57847/251‑1 (MQ=255)
                                                                                                                                                                                               cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCa                                                    <  1:319586/251‑1 (MQ=255)
                                                                                                                                                                                                gCAGGCCCAGCTCGCTCATGCCAGGGTTAGCGGCAAACCACTTATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACAGATTGCGGGTATCACGCACGTCGCGATTCCAGATGCGGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGCTGCAGAAGGGGATTTTAAAGCGTTGTTTGCCGCGGGTGTCCTGCCAg                                                   <  2:273874/251‑1 (MQ=255)
                                                                                                                                                                                                 cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCCAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGc                                                  >  2:271953/1‑251 (MQ=255)
                                                                                                                                                                                                          gtcgctcATGCCAGGGTTATCGGCAAACCACTTATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAGTGTTGCGCGGTATTCGCGCCCGGTGGTCGAGTATCAGAATGGCATATGAACCCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCAt                                         <  1:264361/250‑1 (MQ=255)
                                                                                                                                                                                                            cgctcATGCCAGGGTTATCGGCAAACCACGGATGCTTGTTTTACGGTTCGCTGCCGAGCACTTCCAGCAGCGCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATCCCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGTGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCTGCTGGTTTCATa                                        <  2:293451/250‑1 (MQ=255)
                                                                                                                                                                                                             gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGAAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGGTTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGATTGATAct                                      >  2:224059/1‑251 (MQ=255)
                                                                                                                                                                                                              ctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGTGTATCACGCACGTCGCGATACCAGATGACGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGGTGGTCGAGTAGCAGAAGGAGATTTTAACCCGTTGTTTGCCGCTGGTGTCCTGCCAACTGGTTTCATActc                                     >  2:93767/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccAGGGTTTTCGGCAAACCCCTGATTCTGTTTTTACGGTTCGCCGCCGATCATTTCCAGCCGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAGTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCAc                               <  1:377493/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCAc                               <  1:147992/251‑1 (MQ=255)
                                                                                                                                                                                                                     cAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGGTTCATACTCTGGCACt                              >  1:38745/1‑251 (MQ=255)
                                                                                                                                                                                                                      aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCCGCTCTTCCGGGCGGTGCCACTCGGGAAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTTGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTg                             >  2:36109/1‑251 (MQ=255)
                                                                                                                                                                                                                            aTCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCAt                       <  1:147141/251‑1 (MQ=255)
                                                                                                                                                                                                                            aTCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTACGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGCTGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCAt                       >  1:47472/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccGGCAAACCACTGATTCTGGTTTTACGGTTCGCTTCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATg                      <  2:183216/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACt                  <  1:36109/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttttACGGTTCGCTGCCTAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCTCTCGGGTAAGGCTTTCGACATATTGCGGGTTTCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTTTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAGCCGGTGagtagt  <  2:401490/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGT  >  NC_000913/4585234‑4585725

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