Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,585,673 C→T W363* (TGG→TGA 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,6730CT98.6% 240.3 / NA 73W363* (TGG→TGAhsdRtype I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/39);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCG  >  NC_000913/4585428‑4585916
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGCAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCt                                                                                                                                                                                                                                                >  1:392512/1‑251 (MQ=255)
    ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCGCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTTGTCTTCATTGTTGTGCGGTATTCGCGCGCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGtt                                                                                                                                                                                                                                            <  1:420112/251‑1 (MQ=255)
    ccAGCTCGCTCATGCCAGGGTTACCGGCAAACCACTGATTCTGTTTTTGCGTTTCGCTGCCGAGCATTTCCAGCAGCTCGTCCGGGCGGTGCCACTCGTGTAAGGCTTGCGACATATTGCGGGGATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGGTGGTCGAGTAGCAGAAGGGGATTGTAAACCGTTGTTTTCCGCTTGTGTTCTGTCAGCTGGtt                                                                                                                                                                                                                                            <  2:181017/251‑1 (MQ=255)
     cAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGttt                                                                                                                                                                                                                                           <  1:27596/251‑1 (MQ=255)
     cAGCTCGCCCATGCCAGGGTTATCGGCAAACCACTTATTCTGTTTTTGCGGTTCGCTGCCGCGCATTTCAAGCAGCTCTTCCGGGCGGTGCCACTCGGGAAAGGCTTTCCACCTATGGAGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTTGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGTGTCGCTGACGTGGCTTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGtat                                                                                                                                                                                                                                           >  2:120465/1‑249 (MQ=255)
              cATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCAGTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTTCGGGTATCACGCACGTCGCGATACCAGATGCCGCGTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATActc                                                                                                                                                                                                                                    <  1:429601/249‑1 (MQ=255)
                  ccAGGGTTATCGGCAAACCACTGATGCTGTTGTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTTTCCTGTCAGCTGGTTTCATACTCTGGCAc                                                                                                                                                                                                                              <  1:146260/251‑1 (MQ=255)
                         tATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTc                                                                                                                                                                                                                        >  2:113796/1‑250 (MQ=255)
                                    cACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATc                                                                                                                                                                                                            >  2:14348/1‑251 (MQ=255)
                                      cTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCg                                                                                                                                                                                                          >  2:123457/1‑251 (MQ=255)
                                                     cGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTcaa                                                                                                                                                                                           <  1:332412/251‑1 (MQ=255)
                                                     cGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTACTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGGTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTcaa                                                                                                                                                                                           <  1:228870/251‑1 (MQ=255)
                                                     cGGTTCGCTGCCGAGCATTTCCAGCAGCTATTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTca                                                                                                                                                                                            >  1:330067/1‑250 (MQ=255)
                                                              gCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGt                                                                                                                                                                                    <  1:110491/249‑1 (MQ=255)
                                                                cGAGCATTTCCATCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTGTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGTGGATTTTAAACCGTTGTTTGCCGCCGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGTTTTcc                                                                                                                                                                                <  1:164629/251‑1 (MQ=255)
                                                                  aacattgacaccatctctTCCGCGCGGTGCAACTAGGGTAAGGCTTCAGACCTTTTGCGGGTATCACGCACGTCGCGATACCAGATGCCACTTTTGGTCTTCATTGTTGCGCGGTAGTCGCGCCCGTTGGTCGAGTAGCAGAAGGTGATTTTAAACCGTTGTGTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTccca                                                                                                                                                                              <  2:15861/237‑2 (MQ=255)
                                                                        tccagcagtgtgcccGGGGGGGGCCACCCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCCCGTCGCGATCCCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGTTTTTCGCGCCCGTTGTTGGAGTAGCAGAAGGGAATTTTAAACCGTTGTGTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAATGTTGCCCGCAGGa                                                                                                                                                                         <  2:435881/237‑1 (MQ=255)
                                                                             cagcTCTTCCGGGCGGTGCCCCTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGTCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCTTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTGTCATACTCTGGCACTGCTTCATGCCCTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCAt                                                                                                                                                                   <  1:371618/251‑1 (MQ=255)
                                                                               gcTCTTCCGGGCGGTGCCACTCGGGTAAGGCTGTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCTCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCGTCATGCACTTCATCCGGGGAGTATTGCTCAAGCAAGGTTGCCCGCAGGAAGCCATTa                                                                                                                                                                 <  2:392512/251‑1 (MQ=255)
                                                                                     ccGGGTGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGaaa                                                                                                                                                           <  2:134509/251‑1 (MQ=255)
                                                                                      cGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAAc                                                                                                                                                          <  1:167416/251‑1 (MQ=255)
                                                                                       gggCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCCCGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAAc                                                                                                                                                          >  1:263341/1‑250 (MQ=255)
                                                                                       gggCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAAc                                                                                                                                                          >  1:261662/1‑250 (MQ=255)
                                                                                         gCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAtt                                                                                                                                                       <  2:89903/251‑1 (MQ=255)
                                                                                                cACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGt                                                                                                                                                 >  2:356536/1‑250 (MQ=255)
                                                                                                          aGGCTTTCGACATATTGCGGGTATAACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGCATTCGCGCCCGTTGGTCGAGTAGCAGACGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCACGGTTTCCCGCAGGAAGCCATTATCGAAACATTTCCTGTAGCGCTACGAc                                                                                                                                      >  2:408741/1‑251 (MQ=255)
                                                                                                                 cGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGa                                                                                                                               <  2:366124/251‑1 (MQ=255)
                                                                                                                     atatTGCGGGTATCACGCACGTCGAGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTAGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCATGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCt                                                                                                                           >  2:230497/1‑251 (MQ=255)
                                                                                                                      tatTGCGGGAATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGCGTAGCCGCAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCCGCTGGTTTCATCCTCAGGCTCTTCTTCATGCACTTCATCAGGTGAGTAGTGCTCCAGCAAGGTTTCCCGCCGGAAGCCCTTATCGACACATTTACTGTAGCGTTACGACTCATTGAGCCTg                                                                                                                          >  2:277945/1‑251 (MQ=255)
                                                                                                                       atTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTAGTTGTTGCGCGGTATTCGAGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAAACGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATACGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTgg                                                                                                                         >  2:384232/1‑251 (MQ=255)
                                                                                                                       atTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTGTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGTGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTGTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATGGAGCCTg                                                                                                                          <  2:263341/250‑1 (MQ=255)
                                                                                                                         tGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTAGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACCTTTACTGTAGCGATACGACTCATTGAGCCTGGCg                                                                                                                       >  2:113868/1‑251 (MQ=255)
                                                                                                                           cGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTCGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCCTACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTCATCCTGGCggg                                                                                                                     >  2:201705/1‑251 (MQ=255)
                                                                                                                           cGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATCCTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCAtt                                                                                                                                                                  >  1:407401/1‑206 (MQ=255)
                                                                                                                           cGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCAAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCggg                                                                                                                     <  1:339067/251‑1 (MQ=255)
                                                                                                                           cGGGTATCACGCACGTCGCGAGACCAGATGCCGCTTTGGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCAtt                                                                                                                                                                  <  2:407401/206‑1 (MQ=255)
                                                                                                                              gTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAAc                                                                                                                  <  1:405634/251‑1 (MQ=255)
                                                                                                                                  cacgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCg                                                                                                              >  2:52654/1‑251 (MQ=255)
                                                                                                                                     gcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGata                                                                                                           <  1:99688/251‑1 (MQ=255)
                                                                                                                                                     gATGGCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCTCTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCct                                                                                           <  2:244770/251‑1 (MQ=255)
                                                                                                                                                         ccGCTTTTGGTCTTCAATGTTGCGCGGTATTCGCGCCCGTTGGTCGCGTAGCAGAAGGGGATTTTAAACCGTTGTTTGGCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATTCACTTCATCCGGTGAGTAGTGCTCAATCAAGGTTTCCCGCAGGAAGACATTATCGAAACATTTACTTTACCGATACGACTCATTGAGCCTGGCGGGAACGTCGATAGTGTTACGTTTCGCCTCTa                                                                                        <  2:248484/250‑1 (MQ=255)
                                                                                                                                                                         ttgttgCGCGGTATTCGCGCCCGTTGGTCGAGTAGAAGAAGGGGATTTTAAACAGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGGTTg                                                                       >  2:464634/1‑251 (MQ=255)
                                                                                                                                                                           gttgCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGa                                                                      <  1:464634/250‑1 (MQ=255)
                                                                                                                                                                           gttgCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAg                                                                     <  1:113796/251‑1 (MQ=255)
                                                                                                                                                                             tgCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGc                                                                   >  2:130577/1‑251 (MQ=255)
                                                                                                                                                                                 gggTATTCGCGGCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTGTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTGCGTTTCGCCTCTGCCATCGCGATGGGTTTGAGGCCGAc                                                               <  2:70150/250‑1 (MQ=255)
                                                                                                                                                                                  ggTATTCGCGCCCGGTGGTAGAGTAGCCGAAGGGGGTTTTAAACCGGTGTTTGCCGCTGGTGTCCTGTCAGCGGGTTTCCTACTCTGGCACTGCTTCATGCACTTCATCAGGGGAGTAGTGCTCAAGCACGGTTTCCCGCAGGGAGCCATTATCGCAGCATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGGTATTGTTGCGTTTAGCCTCTACCACCGCGACGGGGTTGAGGCCGAc                                                               >  2:185722/1‑250 (MQ=255)
                                                                                                                                                                                        cgcgCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAg                                                         <  1:356536/250‑1 (MQ=255)
                                                                                                                                                                                        cgcgCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAGCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCTACAAACAg                                                         <  2:231339/250‑1 (MQ=255)
                                                                                                                                                                                         gcgcCCGTTGGTAGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGca                                                       <  1:14348/251‑1 (MQ=255)
                                                                                                                                                                                              cGTTGGTCGAGTCGCAGAAGGGGATTTTCAACTGTTGTTTGCCGCTGTCGTACTGAAAGCTGGTTTTATACTCTGGCCCTGCTTTATGTACTGCATCCGGTGAGCAGTGCTCAAGTACGGTTTCTCGTAGGAAGCCATTATAGAACAATTTAATGTAGAGATATGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGTCTATACTACCGCGATGGGTTTGAGGCCGACTAACAGCACAt                                                    <  1:399237/249‑1 (MQ=255)
                                                                                                                                                                                                 tGGTCGAGTAGCAGAAGGGGATTTTAATCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCAGTATCGAAACATTTACTGGAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACTTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATc                                                <  1:73525/250‑1 (MQ=255)
                                                                                                                                                                                                 tGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCCTCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGACGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAACCAGCACATAATcc                                               >  2:186836/1‑251 (MQ=255)
                                                                                                                                                                                                      gAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCgg                                                                                                                      >  1:341576/1‑175 (MQ=255)
                                                                                                                                                                                                      gAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCgg                                                                                                                      <  2:341576/175‑1 (MQ=255)
                                                                                                                                                                                                          agcagcAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGCAAAACAGCACATAATCCGCAAAGccc                                      >  1:357826/1‑251 (MQ=255)
                                                                                                                                                                                                          aGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGccc                                      >  1:330919/1‑251 (MQ=255)
                                                                                                                                                                                                           gcagacgaggctttaaaaCCGTTATTTGCCGCTGGGGCCCTGCCAGATGGTTTCATACTCTGGCACTGCTTGATGCACTCCATCCGGTGAGAAGTGCTCAAGAAAGGTTTCCCGCAGGAAGCCATTATCGAAACCTTTAATGTAGAGGTACGACTCATTGAGCCTGGCGGGCAAGCCGATATTGTTACGTTTCGCCTATCCCACCGCGATGGGCTTGAGGCCGGCAAACCCAACATAATCCGCAAATcccc                                     >  2:191128/16‑250 (MQ=255)
                                                                                                                                                                                                           gCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCt                                     <  1:176187/251‑1 (MQ=255)
                                                                                                                                                                                                           gCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCt                                     >  2:134055/1‑251 (MQ=255)
                                                                                                                                                                                                            cAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGTAAAGCCCTg                                    <  1:332389/251‑1 (MQ=255)
                                                                                                                                                                                                               aaGGGGAATTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGGTTGAGGCCGACAAACAGCCCATAATCCGCCAAACCCTGAtt                                 >  2:329370/1‑251 (MQ=255)
                                                                                                                                                                                                                       tttAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCAATGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTAGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGCGGCCGCCAAACCGCACATAATCCGCAAAGCCCTGATTACCCGttt                         >  1:39470/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTAGCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATTGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTTCCCGttt                         <  1:445475/250‑1 (MQ=255)
                                                                                                                                                                                                                              cGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTAATACTCTGGCACTGCTTCATGCACTTCATCCGGTGCGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGtttattcttt                  >  2:290366/1‑244 (MQ=255)
                                                                                                                                                                                                                                  gttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCTCCTCTACCACCGCGATGGTTTTGAGTCCGACATACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCtttgcc               <  2:69834/250‑4 (MQ=255)
                                                                                                                                                                                                                                  gttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCg              <  1:197468/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGCTACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGGTGGGTTTGCGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGATTCATCttttacgttcg          >  2:223929/1‑244 (MQ=255)
                                                                                                                                                                                                                                        cGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGc        >  1:322546/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGc        >  1:434354/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGccc      >  1:376688/1‑250 (MQ=255)
                                                                                                                                                                                                                                             gtgtCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCGCAAGCAAGGGTTCCCGCAGGAAGCCAGTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCTGTCGGCCATTc   <  2:330067/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gtCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATAGTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCTCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCg  <  1:384232/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
AGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCG  >  NC_000913/4585428‑4585916

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