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.3% 185.7 / ‑6.3 60W363* (TGG→TGAhsdRtype I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/31);  minor base G (0/1);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCA  >  NC_000913/4585378‑4585968
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                       
tGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTNATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGc                                                                                                                                                                                                                                                                                                     <  2:553262/300‑1 (MQ=255)
        gCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCa                                                                                                                                                                                                                                                                                             >  1:121584/1‑300 (MQ=255)
          cTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATa                                                                                                                                                                                                                                                                                           >  1:336700/1‑300 (MQ=255)
          cTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATa                                                                                                                                                                                                                                                                                           >  2:649563/1‑300 (MQ=255)
                    gcgcgGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCAc                                                                                                                                                                                                                                                                                  <  1:681162/299‑1 (MQ=255)
                    gcacgGACTGCATCTTCCTAATAATAACGCAGGCTCAACTCGCTCATGCAAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGTCGGTGCCACTCGGTTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGATTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGTTTTCATACTCTGacac                                                                                                                                                                                                                                                                                  <  2:208709/296‑1 (MQ=255)
                      gcgGACGGCATCTTTCTGATAATACCGCGGGCCCGGCTCGCTCCTGCCAGGGTTATCGTCAAACCATTGATTTTGTTTTTTCGGTTCTCTGCTGAGCGTTTCCAGCAGCTCTTCCGGTCGGTGCCATTCGGGTAAGGCTTTCGACATATTGCGGGTATCATNCATGTCGCGAGACCAGATGCCGCGTTTGGTCTTCATTGTTGTGCGGTATTCGCGCCCGTTGGTCGAGTAGCATAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCAc                                                                                                                                                                                                                                                                                  <  2:526267/297‑1 (MQ=255)
                            ggCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCAt                                                                                                                                                                                                                                                                          <  2:642248/299‑1 (MQ=255)
                             gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGTTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGc                                                                                                                                                                                                                                                                        <  2:272714/300‑1 (MQ=255)
                             gCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGc                                                                                                                                                                                                                                                                        <  1:576897/300‑1 (MQ=255)
                                  ttCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTTTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTc                                                                                                                                                                                                                                                                   <  2:667215/300‑1 (MQ=255)
                                    ccTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCAt                                                                                                                                                                                                                                                                 >  1:569961/1‑300 (MQ=255)
                                      tGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGANGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATc                                                                                                                                                                                                                                                                <  2:510033/299‑1 (MQ=255)
                                      tGATAATAACACAGGCCCAACTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTTTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATc                                                                                                                                                                                                                                                                <  2:137605/299‑1 (MQ=255)
                                         taataccgcattcccAGCTCGCCCATACCAGGGTGATCGGCAAACCACTAATGCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTATTCCGGGCGGTGCCACTCGGGTAATGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTAGTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTTGCACTGCTTCATGCACTTCATc                                                                                                                                                                                                                                                                <  2:137504/284‑1 (MQ=255)
                                          aataaCACAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCTCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTg                                                                                                                                                                                                                                                           <  2:703688/300‑1 (MQ=255)
                                            gcgagatgtgtataagagacagtGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGa                                                                                                                                                                                                                                                          >  1:653966/23‑299 (MQ=255)
                                               cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTACGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCAGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTACTGCACTTCATCCGGTGagta                                                                                                                                                                                                                                                       >  2:340675/1‑299 (MQ=255)
                                                gCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGagtag                                                                                                                                                                                                                                                      <  2:673947/299‑1 (MQ=255)
                                                   ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCt                                                                                                                                                                                                                                                  >  2:432187/1‑300 (MQ=255)
                                                   ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTAATGCACTTCATCCGGTGAGTAGTGCt                                                                                                                                                                                                                                                  >  2:436710/1‑300 (MQ=255)
                                                    gCCCAGCTCGCTCATGCCAGGTTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTc                                                                                                                                                                                                                                                 <  2:91146/300‑1 (MQ=255)
                                                            cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTcaagcaa                                                                                                                                                                                                                                           >  1:446872/1‑298 (MQ=255)
                                                             actcatgccagaagtacctgCAAACCACTAATGCTGTTTTTGCGGGTCGCTGCCGAGCATTTCCAGCCACTCGCCCAGGCGGTGCCACTCGGGGAAGGCTTTCGACATATTGCGGGTATCACGCACTTCGCGATACCAAATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGGAGGGGATGTTAAACCGTTGTTTGCCGCTAGTGTCCTTTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGagtag                                                                                                                                                                                                                                                      <  2:175274/267‑1 (MQ=255)
                                                                 aTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTcc                                                                                                                                                                                                                                    >  1:165637/1‑300 (MQ=255)
                                                                 aTGCAAGGGTTATCGGCAAACCACTGATTCTGTTTTTGAGGTTCGCTGCCGAGCATTTCCAGCAGCTATTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGAATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTACTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTcc                                                                                                                                                                                                                                    >  2:79425/1‑300 (MQ=255)
                                                                  tGCCAGGGTTATCGGCNAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCNCTTCATCCGGTGAGTAGTNCTCAAGCAAGGTTTccc                                                                                                                                                                                                                                   >  1:603176/1‑300 (MQ=255)
                                                                  tGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTccc                                                                                                                                                                                                                                   >  1:227298/1‑300 (MQ=255)
                                                                                                                       aTTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCTTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCt                                                                                                                                                                               >  1:425020/1‑299 (MQ=255)
                                                                                                                                    ctTCCGGGCGTTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGc                                                                                                                                                                                 <  2:701529/284‑1 (MQ=255)
                                                                                                                                    ctTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGc                                                                                                                                                                                 >  1:701529/1‑284 (MQ=255)
                                                                                                                                    ctTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGa                                                                                                                                                                 >  1:562858/1‑300 (MQ=255)
                                                                                                                                    ctTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCATGGCGGGAACGTCGa                                                                                                                                                                 >  1:561415/1‑300 (MQ=255)
                                                                                                                                              gTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTAc                                                                                                                                                        >  1:463273/1‑299 (MQ=255)
                                                                                                                                                          taatgCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATTCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCANTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTac                                                                                                                                           <  2:523835/296‑1 (MQ=255)
                                                                                                                                                             ggCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTaccac                                                                                                                                        <  2:463273/300‑1 (MQ=255)
                                                                                                                                                             ggCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTaccac                                                                                                                                        <  1:182514/300‑1 (MQ=255)
                                                                                                                                                                       atatTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCTGGAACGTCGATATTGTTACGTTTCGCCTCTACTACCGCGATggg                                                                                                                               <  2:485734/299‑1 (MQ=255)
                                                                                                                                                                              gggTACCACGCACTTCCCGATACCAGATGCCTCGTTTGTTCTTCATTTTGGCGCGGTGTTCGCGCCCGTTGTTCGAGGAGCAGAAGGGGATGTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTTATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGCGCTCAAGCAAGGTTTCCCGCAGTAAGCCATTACCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGc                                                                                                                       <  2:258088/300‑1 (MQ=255)
                                                                                                                                                                               ggTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGc                                                                                                                       >  2:440213/1‑299 (MQ=255)
                                                                                                                                                                                  aTCATGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGTTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGATTAGTGCTCAAGCAAGGTTTCCCGCAGTAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCTTGAACGTCGATATTGTTCCGTTTCGCCTCTACCACCGCGATGGGTTTGAGGc                                                                                                                       <  2:258476/296‑1 (MQ=255)
                                                                                                                                                                                  aTCACGCACGTCGCGANACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACNTTTCGCCTCTACCACCGCNATGGGTTTGAGGCCGAc                                                                                                                   >  1:593933/1‑300 (MQ=255)
                                                                                                                                                                                          cgTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGca                                                                                                           >  1:624340/1‑300 (MQ=255)
                                                                                                                                                                                             cgcgATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACAt                                                                                                        >  1:558771/1‑300 (MQ=255)
                                                                                                                                                                                                      aGATGCCGCTTTTGGCCTTCATTGTTGTGCGGCATTCGCGCCCGTTGTTCTATTAGCAGAATTGAATTTTAAACCGTTGTTTGCCGTTGGTGTCCTTTCAGCTGGTTTCATACTCTGTCACTGCTCCATGCACTTCATCCGGTGAGTAGTACTTAAGCAAGGTTTCCCGCAGGAAGCCATTAGCGAAACATTTACTGTAGCGATACGACTCATTGAGTCTGTCGGGAACGTCGATATTGTTACGTTTCGCTTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCTGCaa                                                                                               <  2:401419/300‑1 (MQ=255)
                                                                                                                                                                                                                tttGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGAtt                                                                                     <  2:603176/300‑1 (MQ=255)
                                                                                                                                                                                                                ttgggTCTTCATGGTTGCGCGGTATTCGCGCCCGTTGTTCGAGTAGCAGAAGGGGATTCTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTNAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCTGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGAtt                                                                                     <  2:227298/297‑1 (MQ=255)
                                                                                                                                                                                                                 ttcgtattCATTGTGGCGTGGTATTTACGCCCGTTGGTCGAGTAGCAGAGGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGGCAAATGGTTTCAGACTCTGGTACTGCTTTAGGCACTTCATCCGGTGAGTAGTGCTCAAGTAAGGTTTCCCGCAGGAAGCCATTAGCGAAACATCTACTGTAGCGATACGACTCATTGAGCCCGGCGGGAACGGCGACATTGTTACGTTTCGCCTCTATCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGAtt                                                                                     <  2:18084/293‑1 (MQ=255)
                                                                                                                                                                                                                   ggTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTAc                                                                                   >  2:332781/1‑299 (MQ=255)
                                                                                                                                                                                                                      cTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGt                                                                               >  1:445043/1‑300 (MQ=255)
                                                                                                                                                                                                                                     gTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCATGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGACATGATTACCCGTTTCATCTTCTcc                                                                   >  2:353466/1‑297 (MQ=255)
                                                                                                                                                                                                                                       aTTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTc                                                               >  1:620471/1‑299 (MQ=255)
                                                                                                                                                                                                                                                cGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTc                                                       >  1:532144/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                             aaCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTTGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGNTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGa                        <  2:562858/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                             aaCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGNTTTCCCGCAGGAAGCCATTATCTAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGa                        <  2:561415/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                             aaCCGTTGTTTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGa                        <  2:314412/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    gttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCNCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTgcgc                 <  2:569961/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    gttTGCCGCTGGTGTCCTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTgcgc                 <  1:564204/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             tGGTGTCCTGTCAGCTGTTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGgag         <  1:679366/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   ccTGTCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCa  >  2:430354/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                       
TGACGATTGCCTTTTCAACCGCGCGGACGGCATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGGTATCACGCACGTCGCGATACCAGATGCCGCTTTTGGTCTTCATTGTTGCGCGGTATTCGCGCCCGTTGGTCGAGTAGCAGAAGGGGATTTTAAACCGTTGTTTGCCGCTGGTGTCCTGCCAGCTGGTTTCATACTCTGGCACTGCTTCATGCACTTCATCCGGTGAGTAGTGCTCAAGCAAGGTTTCCCGCAGGAAGCCATTATCGAAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCA  >  NC_000913/4585378‑4585968

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