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% 219.1 / NA 70W363* (TGG→TGAhsdRtype I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/35);  minor base G (0/1);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCG  >  NC_000913/4585427‑4585916
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGCGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGCTGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGc                                                                                                                                                                                                                                                 >  1:585814/1‑251 (MQ=255)
 aGGCCCAGCTCGCTCATGCCTGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCt                                                                                                                                                                                                                                                <  1:429180/251‑1 (MQ=255)
 aGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCt                                                                                                                                                                                                                                                >  1:198685/1‑251 (MQ=255)
  ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTg                                                                                                                                                                                                                                               <  1:209749/251‑1 (MQ=255)
      cAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAGTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGttt                                                                                                                                                                                                                                           <  1:254612/251‑1 (MQ=255)
           cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCAt                                                                                                                                                                                                                                        <  1:175458/249‑1 (MQ=255)
             ctcATGCCAGGGTTATCGGCAAACCAATGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGCTGGTCGGGTCGCAGAAGGGGGTTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGGTTCATActc                                                                                                                                                                                                                                    >  2:330149/1‑251 (MQ=255)
                     aGGGTTATCGGCAAACAACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTACGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTg                                                                                                                                                                                                                            >  2:278695/1‑251 (MQ=255)
                      gggTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGc                                                                                                                                                                                                                           >  2:428551/1‑251 (MQ=255)
                             cGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGc                                                                                                                                                                                                                    >  2:271878/1‑251 (MQ=255)
                               gCAAAACACTGATTCTGTTTTTGAGGTTAGATGCAGCGCATTACCAAAAGCTATTCCGGGCGGTGCCACTCGGGAAAGGCTTACGACCCATTGCGGGTATCACGCACGTCCCAATACCAGATGCCGCTTTTGGTATTCATTGTTGCGCGGCATTCGCGCCCGTTGGTCGAGTCGCAGAAGGGAATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGATGGTTTCATACTCTGGTACTGCTTCATGCa                                                                                                                                                                                                                   >  2:375851/1‑250 (MQ=255)
                                   aCCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCa                                                                                                                                                                                                              <  2:13163/251‑1 (MQ=255)
                                      aCTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATcc                                                                                                                                                                                                           >  1:5161/1‑251 (MQ=255)
                                       cTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCTCTTTTGGTCTTCATTGTTTCGCGGGATTCGCGCCCGTTGGTCGAGGAGCAGAAGGGGATTTTCAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCCGGCACTGCTTCATGCACTTCATCCg                                                                                                                                                                                                          >  2:403504/1‑251 (MQ=255)
                                       cTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCg                                                                                                                                                                                                          >  1:211703/1‑251 (MQ=255)
                                            tCTGTTTTTGCGGGTCGCTGCCGAGCTTGTCCAGCAGCGCTACCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATGCCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAATGGTATTTTAAACCGTTGTTTGCCGCTTGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGa                                                                                                                                                                                                      <  2:131977/250‑1 (MQ=255)
                                                  tttGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGCCATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGc                                                                                                                                                                                               >  1:249654/1‑251 (MQ=255)
                                                    tGCGGTTCGCTGCCGAGCATTTCCAGTAGCTCTTCCGGGCTGTGCCACTCGAGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTGGGTCGAGGAGCAGAAGGGGATTTTAAACCGTTGTTTGGCGCTGGTGTCCTGTCAGCTGGTGTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCt                                                                                                                                                                                              <  2:313414/250‑1 (MQ=255)
                                                     gCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGCAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTc                                                                                                                                                                                             >  1:84029/1‑250 (MQ=255)
                                                      cGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTcaa                                                                                                                                                                                           <  2:10522/251‑1 (MQ=255)
                                                      cGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGAATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTca                                                                                                                                                                                            >  1:589101/1‑250 (MQ=255)
                                                        gTTCGCTGCCGAGCATTTCCCGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGGGAGTCGTGCTcacgc                                                                                                                                                                                         >  2:461188/1‑248 (MQ=255)
                                                              tGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGtt                                                                                                                                                                                   >  2:133331/1‑251 (MQ=255)
                                                              tGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGtt                                                                                                                                                                                   >  2:231224/1‑251 (MQ=255)
                                                                ccGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTc                                                                                                                                                                                 <  1:334882/251‑1 (MQ=255)
                                                                  gAGCATTTCCAGCAGCTCTTCCCGGCGGTGCCACTCGGGTAAGGCTTTCGCCATATTGCGGGAATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGACGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCGCAAGCAAGGTTTccc                                                                                                                                                                               >  2:521977/1‑251 (MQ=255)
                                                                     cattgaccgcatctctGCCGGGCGGTGCCACTCTGTTAAGGCTTGCGACTTATTGCGATTATCACGCACGTCGCGATACCTGATTCCGCTTTTGGTCTTCATTGTGGCGCTGTATTCGCGCGCGTTTGTCGAGTAGCATAAGGGGATTTTAAATCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGGAGGGCTCAAGCAAGTTTTCCCGCa                                                                                                                                                                            <  2:135452/239‑1 (MQ=255)
                                                                               agcTCGTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCAGTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCGGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCAtt                                                                                                                                                                  <  1:575965/251‑1 (MQ=255)
                                                                                  tctTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAACCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGGTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGACGCCATTATc                                                                                                                                                               >  2:461472/1‑251 (MQ=255)
                                                                                      ccGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGaaa                                                                                                                                                           <  2:235324/251‑1 (MQ=255)
                                                                                      ccGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCGGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGaa                                                                                                                                                            <  1:68223/250‑1 (MQ=255)
                                                                                        gcgCGGTGCCACTCGGGTAAGGCTTTTGACATATTGCGGGTATCACGCTCGTCGCGATACCAGATGCTGCGTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGTTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACa                                                                                                                                                         <  2:561402/249‑1 (MQ=255)
                                                                                           gggTGCCACTCGGGTAAGGCCGTCCACATATTGCGGGTATCACGCGCGGCGCGATACCAGATGCCGCTTTTTGGCGTCATTGTTTCGCGGTAGTCGCGCCCGTTGGTCGGGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTTGTTTCATACTCTGGCACTGCTTCATGCACGTCATTCGGGGAGTTGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACAttt                                                                                                                                                      <  1:97479/250‑1 (MQ=255)
                                                                                                ccACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGt                                                                                                                                                 >  1:448550/1‑251 (MQ=255)
                                                                                                      gggaaaGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGAt                                                                                                                                           <  2:393080/247‑1 (MQ=255)
                                                                                                       ggTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATa                                                                                                                                          <  1:338473/251‑1 (MQ=255)
                                                                                                        gTAAGGCTTGCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATAc                                                                                                                                         <  1:323045/251‑1 (MQ=255)
                                                                                                           aGGCTTTCTGCAGATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTTTCACTGCTTCATGCACTTCATCCGGGGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGGAGCCATTATCGAAACATTTACTGGAGTGATACGAc                                                                                                                                      <  2:514917/251‑1 (MQ=255)
                                                                                                           aGGCTTTCGTCAGATTGCGGTTATCACGCACGTCGCGATACCAGATGCCGCTTGTGGTCGTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCTAGTAGCAGAAGGGGAGTTTAAACCGTTTTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAGCCGGGGAGTAGTGCTCAAGCAAGGTTTCCCTCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGAc                                                                                                                                      <  2:226038/251‑1 (MQ=255)
                                                                                                                      atatTGCGGGTATCAAGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCt                                                                                                                           >  1:96190/1‑251 (MQ=255)
                                                                                                                         ttGCGGGTAGCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTgg                                                                                                                         <  2:37382/250‑1 (MQ=255)
                                                                                                                              ggTATCACGCACGTCGCGATACCAGAAGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTAAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGa                                                                                                                    >  2:419505/1‑250 (MQ=255)
                                                                                                                               gTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAATCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAAc                                                                                                                  <  1:515802/251‑1 (MQ=255)
                                                                                                                                  tCACGCACATCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCCGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGCGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTc                                                                                                               >  2:278014/1‑251 (MQ=255)
                                                                                                                                   cacgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTc                                                                                                               <  1:196039/250‑1 (MQ=255)
                                                                                                                                     cgcacgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGat                                                                                                            >  2:12804/1‑251 (MQ=255)
                                                                                                                                         cgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTg                                                                                                        >  1:459806/1‑251 (MQ=255)
                                                                                                                                                aTACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGtt                                                                                                 <  1:464054/251‑1 (MQ=255)
                                                                                                                                                                     cTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGCGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGACACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGCTATTGTTCCGTTTCGCCTCTACCACCGCGATgg                                                                            >  2:251693/1‑251 (MQ=255)
                                                                                                                                                                            gttgCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTg                                                                       >  1:130440/1‑249 (MQ=255)
                                                                                                                                                                                    gTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGATTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATACGGTGAGTAGTGCTCAAGCACGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATCTTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGCTGCCGAcc                                                              >  2:597096/1‑249 (MQ=255)
                                                                                                                                                                                      aTTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGGAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAAc                                                           <  1:398717/251‑1 (MQ=255)
                                                                                                                                                                                         cgcgCTCGTTTGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAGCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCAGTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAATCAGc                                                        <  1:449101/251‑1 (MQ=255)
                                                                                                                                                                                          gcgcCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGc                                                        >  1:225772/1‑250 (MQ=255)
                                                                                                                                                                                            gcCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGcac                                                      >  1:284562/1‑250 (MQ=255)
                                                                                                                                                                                             cccGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACAt                                                    >  2:112945/1‑251 (MQ=255)
                                                                                                                                                                                                  tGGTCGAGTCTCAGAAGGTGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTTGTTTCAGACTCTGGCACTGCTTCATGCACTGCATCCGGTTAGTAGTGCTCAAGCAAGGTTTCCCGCAGTAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTGGAGGCCGACAAACAGCACATAATc                                                <  2:175910/250‑1 (MQ=255)
                                                                                                                                                                                                       gAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTTGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAATCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGTGAACGTCGATATTGTTACGTTTCTCCTCGTCCATCGCGATGGGTTTGATTCCGACAAACAGCACATAATCCGCaaa                                          <  2:513963/251‑1 (MQ=255)
                                                                                                                                                                                                              aGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGa                                   <  1:313045/251‑1 (MQ=255)
                                                                                                                                                                                                                        tttAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGtt                          <  1:428551/250‑1 (MQ=255)
                                                                                                                                                                                                                             aCCGTTGTTTGCCGCCGGTGACCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGGCGCTATTGTTAAGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAAAAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCt                    >  2:501752/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggtgttGGCCGCTTGTGGCCTGTCAGCTGGTTTCATTCTCTTGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCGCAAGCAATGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTTCCACCTCGATTTGTGTGATGCCGACATACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCtatt                 <  2:145326/249‑4 (MQ=255)
                                                                                                                                                                                                                                  tgttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTcc               >  1:29320/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGTTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTcc               <  1:597096/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTcc               <  1:12804/250‑1 (MQ=255)
                                                                                                                                                                                                                                    tttGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATCTTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAACGCCCTGATTACCCGTTTCCTCTTTTACgg             >  2:383457/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGcc       <  1:133331/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tgtgtCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTCGTGCTCAAGCAAGGTTTCCCGCAGCATGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGTTATTGTTCATTTTCTCATCTCCCACCGCGATGGGTTTCAGCCCGCCAAAAAGCAACTAATCTGCAAACCCCTACTAACCCGATTCATCTTTTCCGGTCCGCCAtc    >  2:4748/2‑250 (MQ=255)
                                                                                                                                                                                                                                              gtgtCTTGGCAGCTGGTTTCATACTCTGGCACTTCTTCCTGCTCTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGATGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTTAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACCGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCCTATCCGGTCGGCCactc   <  2:692/251‑4 (MQ=255)
                                                                                                                                                                                                                                                gtCCTGTCAGCTGGTTTCATTCTCTGGCACTGCTTCATGCACGTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATAGTGTGACGTTGCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTAATGACCCGTTTCATCTTTTCCGGGCTGCCATTCg  <  2:228890/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCG  >  NC_000913/4585427‑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: