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.1% 160.7 / NA 54Q428* (CAA→TAA) hsdRtype I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (24/28);  minor base C (1/0);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTA  >  NC_000913/4585234‑4585723
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
ccTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                                                                                                                                                                                                 <  1:303312/251‑1 (MQ=255)
  tGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATAGCCAGCAGGATATCTTGTTGCCCCTTGACGATTGCCTTTTAAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCAGCCAACCACTGATTCTGTTTTTACGGTTc                                                                                                                                                                                                                                               >  1:259125/1‑251 (MQ=255)
  ggTTCGCCAAGAGAACGGTGTTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCATCCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTTATTCTGTTTTTACGGTTc                                                                                                                                                                                                                                               <  2:177420/250‑1 (MQ=255)
   gTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGGACGGGTTTTAACGGTACCGGTCGCCATCGCCCGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGCCGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCCGGGTTATCGGCAAACCCCTGATTCTGTTTTTACGGTTCg                                                                                                                                                                                                                                              >  1:67073/1‑251 (MQ=255)
      cGCCAAGAGAACGGCGGCATACACCGGAGAGAATGCCTTTAACACGATGGGGCCGGGTCAGGCGGCACATCATGGCGATTGCCGTACGGGTTTTACCGGGACCGGTCGCCATCGCCAGCAGGATTTCTTGGTACCCCTTGACGCTTGCCTTTTCCACCGCGCCGCCCGCATTTTCCTGCTAATAACGCCGGCCCAGCTCGCCCCTGCCAGGGTTTTCTGCAACCCCCTGCTTCTGTTTTTCCGGTTCGCtc                                                                                                                                                                                                                                           >  2:164032/1‑250 (MQ=255)
        ccAAGAGAACGGCGGTCGACAAGGAAGAGAATGAGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCCGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGcc                                                                                                                                                                                                                                         >  1:374127/1‑251 (MQ=255)
        ccAAGAGAACGGCGGTCGACAAGGAAGAGAAAGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCCTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGcc                                                                                                                                                                                                                                         >  1:449701/1‑251 (MQ=255)
          aaGAGAACGGCGGTCGACAAGGGAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCCTCTTCCTGATAATAACGCAGGCCCCGCTCGCTCATGCCAGGGGTATCGGCAAACCACTGATTCTGTTTTTACGGTTCgccgcccg                                                                                                                                                                                                                                       >  2:185506/1‑249 (MQ=255)
              gaACGGCGGTCGACAAGGAAAAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGGTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAt                                                                                                                                                                                                                                   >  1:247669/1‑251 (MQ=255)
                    cggTCGACAAGGAAGATAATGCTTTTAAAACGCTGGGACTGGAGCAGGCGGAAAATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTGTTCAACCGCGCGGACGGCATCGTCCTGATAATAACGCAGGCCCAGCTCGCTCCGGCCAGGGTTATCGGCAAACCACTGATTCTGTGTTTGGGGTTCGCTGCCGAGCTTTTCcag                                                                                                                                                                                                                             <  2:229375/251‑1 (MQ=255)
                     ggTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACAACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagc                                                                                                                                                                                                                            >  2:74301/1‑251 (MQ=255)
                         gACAAGGAAGAGAATGCGTTTAAAACGCTGTGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTACCGGCAAACCACTGATTCTGGTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                                                                                                                                                                                                        <  1:391308/251‑1 (MQ=255)
                          aCAAGGAAGAGAATGCGTTTCAAACGCTGGGACTGGATCAGGCGGACCATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCCGGCCCAGCTCGCTCATGCCAGGGTTATCGGCCAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                                                                                                                                                                                                       >  2:269219/1‑251 (MQ=255)
                             aGGAAGATAATTCGTTCAAAACGCTGGGACTGGAGTAGGCGGAACATCATGGCGATTTCCGTACTGGTTTTACCGGTACCTGTCGCCATCGCCAGCAGGATTTCTTTTTTCTCCTTTACGATTTCCTTTTCATCCGCGGGGACGTCATCTACCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGTGTTATCGTAAAACCACTGATGCTGTGTTTACGGTTCGCTGCCGAGCAGGTCCAGCCGCTCttg                                                                                                                                                                                                                    <  2:356291/251‑2 (MQ=255)
                                  gagaATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCGTCCGGGc                                                                                                                                                                                                               <  2:62331/251‑1 (MQ=255)
                                     aaTGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACGGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGt                                                                                                                                                                                                            <  1:416369/251‑1 (MQ=255)
                                        gCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTGTACCGGTACCGGTCGCCATCGCCAGCAGGAGCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCGGTTGTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCGTCCGTGCGGTGcc                                                                                                                                                                                                         <  2:215596/251‑1 (MQ=255)
                                               aaaCGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                                                                                                                                                                                                  >  1:204113/1‑251 (MQ=255)
                                               aaaCGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                                                                                                                                                                                                  >  1:283909/1‑251 (MQ=255)
                                                 aCGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTa                                                                                                                                                                                                <  1:61648/251‑1 (MQ=255)
                                                        gACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCttt                                                                                                                                                                                         <  1:371388/251‑1 (MQ=255)
                                                             gATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTCCCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGAAAACAACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGGGCCACTCGGGGAAGGCTTTCgcc                                                                                                                                                                                     >  2:415833/1‑248 (MQ=255)
                                                                 aGGCGGAACATCATGGCGATAGCCGTACGGGCTTTACCGGTACCGGACGACATCGCCAGCAGGAACTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGACCAGCTCGGTCATGCCAGGGTTATCGGCAAAACACTGATTCTGTTTTTACGATTCGCTTCCGAGCATTTCCAGCAACTCTTCCGGGCGGTGCCACCAGGGTAAGGCTTTCGACATAtt                                                                                                                                                                                >  2:341551/1‑251 (MQ=255)
                                                                  ggCGGAACATCATCGCGCTTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCCGCAGGATCTCTTGTTGACCCTTGACAATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCCGCTCGCTCATGCCAGGGTTCTCGGCAAACCACTGCTTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCTCTCGGGTAAGGCTTTCGACATATTg                                                                                                                                                                               >  2:431079/1‑251 (MQ=255)
                                                                   gcgtaaCAGCATGGCGATTGTCGTACGGGTTTTACCGGTACCGGTCGCCAGCGCCTGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGTGTGATCGGCAAACCACTTATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCATCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGc                                                                                                                                                                              <  2:407682/247‑1 (MQ=255)
                                                                               ggcggtgtccGTACGGGTTTTACCGGTGCCGGTCGCCATCGCCATCAGGATCTCTTGTTGCGCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCGTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACGGAGTCGGTTGTTACTGTGCGCGGCCGAGCATTTCCAGTAGCTCTTCCTGGCGGTGTCACTCGGTTAATGCTGTCGACATATTGCGGGGAGcacgca                                                                                                                                                                  <  2:425258/243‑1 (MQ=255)
                                                                                   aTTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTc                                                                                                                                                              <  1:94414/251‑1 (MQ=255)
                                                                                    ttGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCTAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGCCATATTGCGGGTATCACGCACGTcg                                                                                                                                                             >  1:208979/1‑251 (MQ=255)
                                                                                                  ttACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGc                                                                                                                                                <  2:67073/250‑1 (MQ=255)
                                                                                                       ggtactggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCttt                                                                                                                                          <  1:380222/251‑1 (MQ=255)
                                                                                                             ggtCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCt                                                                                                                                    >  1:244911/1‑251 (MQ=255)
                                                                                                                cgccatcgccaGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTTGTCTTCa                                                                                                                                 <  1:335914/251‑1 (MQ=255)
                                                                                                                 gccatcgccaGCAGGATCGCTTGTTGCCCCTTGACGATTGCCTTTGCAACCGCGCTGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGTGTTATCGGCAAACCACTGATTCTTTGTTTACGGTTCGCTGCCTAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAt                                                                                                                                <  1:448537/251‑1 (MQ=255)
                                                                                                                    atcgccaGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAttgt                                                                                                                             <  2:192770/251‑1 (MQ=255)
                                                                                                                          agcagGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGcgcg                                                                                                                       >  1:226944/1‑251 (MQ=255)
                                                                                                                               gATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTAtt                                                                                                                  >  2:331678/1‑251 (MQ=255)
                                                                                                                                       gttgCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCg                                                                                                          <  1:104899/251‑1 (MQ=255)
                                                                                                                                                   cGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGtcg                                                                                              >  1:302794/1‑249 (MQ=255)
                                                                                                                                                            tttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgg                                                                                       <  1:185506/249‑1 (MQ=255)
                                                                                                                                                            tttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCAGGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAggg                                                                                      <  2:449701/250‑1 (MQ=255)
                                                                                                                                                             ttCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATGCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCGTCATGGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAgggg                                                                                     <  2:283381/250‑1 (MQ=255)
                                                                                                                                                              tcatccGCGCGGACGGCATCTTCCTGATAATTACGCAGGCCCAGCTCGCGCATGCCAGGGTTATCGGCAAACCGCTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGTCTTTCGATATATTGCGGGTATCACGCACGTCGCGATGCCAGATGCCGGTTTGGGTCTTCATTGTTGGGCTGTGTTCGCGCCCGTTGGTCGATTAGCAGAAGGGGa                                                                                    <  1:261932/246‑1 (MQ=255)
                                                                                                                                                                    gcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTaaa                                                                             <  1:269219/251‑1 (MQ=255)
                                                                                                                                                                            ggCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCgttgtt                                                                     >  2:76640/1‑251 (MQ=255)
                                                                                                                                                                             gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAtt                                                                                                                               <  1:355918/192‑1 (MQ=255)
                                                                                                                                                                             gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAtt                                                                                                                               >  2:355918/1‑192 (MQ=255)
                                                                                                                                                                                   tCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCCGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCt                                                              >  2:348523/1‑251 (MQ=255)
                                                                                                                                                                                                 cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGc                                                >  1:262192/1‑251 (MQ=255)
                                                                                                                                                                                                   ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTg                                              <  1:415833/251‑1 (MQ=255)
                                                                                                                                                                                                   ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCGGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGTTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATGTTAAAGCGTTTTTTGCCGCTGGTGTCCTGCCAGCTg                                              <  2:364480/251‑1 (MQ=255)
                                                                                                                                                                                                                           tATCGGCAAATCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCTTTTCCAGCAGGTCTTCCGTGCGGTGCCTCTCGGGTTAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGTTTTTGGTCTTCATGGTTGCGCGGTATGCGCGATCTTTGGTCGAGTAGCTGAAGGCGATGTTAAACCGTTGTGTGCCGCTGGTGGCCTTCCAGCTGGTTTCATACTCTGGCACTGcct                        <  2:305870/249‑3 (MQ=255)
                                                                                                                                                                                                                                       aCTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATcc          >  1:103070/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCg         >  1:133293/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggTTTTTACGGTTCGATGCCGAGCACTTCCAGCAGCTCTTCCGGGCGGTGCCAGTCGGGTAAGGCTTTCGGCATATTGCGGGTATCACGCACGTCGCGATAGCAGATGCCGCTTTTGGTGTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGTATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCGGCTGGTTTCAGACTCTGGCACTGCTACATGCACTTCATCCGGTGagta  <  2:69696/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CCTGTTCGCCAAGAGAACGGCGGTCGACAAGGAAGAGAATGCGTTTAAAACGCTGGGACTGGATCAGGCGGAACATCATGGCGATTGCCGTACGGGTTTTACCGGTACCGGTCGCCATCGCCAGCAGGATCTCTTGTTGCCCCTTGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTA  >  NC_000913/4585234‑4585723

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: