Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 644,042 C→T 9.5% intergenic (‑75/+155) rnk ← / ← rna regulator of nucleoside diphosphate kinase/ribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913644,0420CT 9.5% 162.2 / 4.9 58intergenic (‑75/+155)rnk/rnaregulator of nucleoside diphosphate kinase/ribonuclease I
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base T (4/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCT  >  NC_000913/643743‑644342
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                            
cgaagtgcgtgtgcGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGCAAACCCTCACCGTCAGgcggcg                                                                                                                                                                                                                                                                                                             >  1:491868/14‑301 (MQ=255)
    gTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAggg                                                                                                                                                                                                                                                                                                         <  2:602149/301‑1 (MQ=255)
            acGCACTTCGCAATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCGTCGTGCGAACACATTTGGGCGGGATCCAACGCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTc                                                                                                                                                                                                                                                                                                 <  1:551357/301‑1 (MQ=255)
               cACTTCGCCATCGCTAATATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTcaca                                                                                                                                                                                                                                                                                              <  2:385865/301‑1 (MQ=255)
                             tAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTATTGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATactgact                                                                                                                                                                                                                                                                                >  2:566970/1‑301 (MQ=255)
                                   ggctgAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCt                                                                                                                                                                                                                                                                          <  1:431225/297‑1 (MQ=255)
                                    gCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCta                                                                                                                                                                                                                                                                         >  2:431789/1‑300 (MQ=255)
                                      ggAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGAAAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCGGCGAGTGGCGCAGATAAACATACCCAAAAGAAAACCCTCACCGCCAGGCGGCGAGGGTTTAACTCACATAATGATACTGACTGTTGCTCAc                                                                                                                                                                                                                                                                       >  1:518523/1‑301 (MQ=255)
                                               cccGGTTGTTTATTGGTACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATAGCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg                                                                                                                                                                                                                                                              <  1:298489/301‑1 (MQ=255)
                                                      gTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTACAATGATGATCCTGACTGTTGCTCACTCATTAAAGTGAtttt                                                                                                                                                                                                                                                       >  1:436686/1‑300 (MQ=255)
                                                        tcgttGTCACCTGGTGGTGTGGTATCGTTTCTGGCGAACACATTTTGGCGTGATCCAACGTTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCg                                                                                                                                                                                                                                                     <  1:70253/298‑1 (MQ=255)
                                                        tCATTGTCACCACTTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCg                                                                                                                                                                                                                                                     <  2:489134/301‑1 (MQ=255)
                                                         cATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCCAGGGTTTAACTCACATGAAGAAACTGACTGTTTCCCACTCTTTCAAATTCTTTGCGt                                                                                                                                                                                                                                                    >  1:639977/1‑301 (MQ=255)
                                                          ttccgatctacGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTc                                                                                                                                                                                                                                                   <  2:548210/292‑1 (MQ=255)
                                                                accacGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTc                                                                                                                                                                                                                                             >  2:295714/1‑301 (MQ=255)
                                                                   acGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCAGAtcggaa                                                                                                                                                                                                                                          >  1:547537/1‑296 (MQ=255)
                                                                                            aaCTCATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGTACATTTTTACTCCATGTCGTCGGTGCTTCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATg                                                                                                                                                                                                                 <  2:437250/301‑1 (MQ=255)
                                                                                              cacaTTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGTCATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGta                                                                                                                                                                                                               <  2:612219/301‑1 (MQ=255)
                                                                                                                   gctGCGGTTTAGGCGTCGGCGATTGGCAAACCAGCTTAGGCGGGTTGGTCCAGTGGAAGATCGATGCGTTCGGCATCCAGGTAGTTAATGATGATAGTTGGTCTGGACATTGTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCATCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCGGCGCATTATGTATAAATCTTAATCGCCTTGGtt                                                                                                                                                                                          <  1:566256/300‑1 (MQ=255)
                                                                                                                        ttttAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATgg                                                                                                                                                                                     <  2:330293/300‑1 (MQ=255)
                                                                                                                           tAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGaaa                                                                                                                                                                                  >  2:311649/1‑300 (MQ=255)
                                                                                                                             aCGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAACAGAAAACCCTCACCGTCAGGAGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGACGCGCATTATGTAACAATCTTAATCACCCTGGTTTATGAAAAAc                                                                                                                                                                                >  1:484430/1‑301 (MQ=255)
                                                                                                                                    ggCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACACTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGGAAAAACCTTAACCACCTATGTTTATGGAAGACGAAttta                                                                                                                                                                         >  1:15311/1‑299 (MQ=255)
                                                                                                                                        aTTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTAGGGAAGAAGAATTAGCgtg                                                                                                                                                                     >  2:381017/1‑301 (MQ=255)
                                                                                                                                         tggagcaaCCCGCCTATGCTGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAACGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATCTGCGTCACATTCAGGGAATTCCTCAATGCAGCTCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAATACAAATTACCgtgt                                                                                                                                                                    >  1:411495/7‑301 (MQ=255)
                                                                                                                                           gagcaaCCCGCCTATGCTGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATCTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGttt                                                                                                                                                                  <  2:412058/297‑1 (MQ=255)
                                                                                                                                                         ggtgggTTGTTCCAGCGGAATTTCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAgatga                                                                                                                                                    <  1:295310/298‑1 (MQ=255)
                                                                                                                                                                 gCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATtaataaa                                                                                                                                            >  2:389396/1‑300 (MQ=255)
                                                                                                                                                                 gCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAAc                                                                                                                                            <  2:631915/301‑1 (MQ=255)
                                                                                                                                                                  cTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATATTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAAcc                                                                                                                                           >  1:502976/1‑301 (MQ=255)
                                                                                                                                                                         gggagtggggtgggTTCGGGGTCGAGGTGGTTAGTGTTGGTATGTGGTCTGGACAGTTTTGCTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCGAAAGAAAGCCCTCACCGTCAGGCGGCGAGGGTTTAACGCACATGAGGATACTGACTGTTGCTGACTCTTTGGAGGGATGTGCGTCAGATGCAGGGAATTCCTCAGTGCAGGGCATTATGTATGAGTCTTAGTCGGCTTGGGTTATGGAGGGGGAGTGGGCGTGTGTGGTAGAGCGGATGGTTGATAGCCCgggggg                                                                                                                                    <  1:311245/288‑6 (MQ=255)
                                                                                                                                                                            atatCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTatca                                                                                                                                 <  2:492514/301‑1 (MQ=255)
                                                                                                                                                                                   tGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACaa                                                                                                                          <  2:503648/301‑1 (MQ=255)
                                                                                                                                                                                           gCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGttttt                                                                                                                  >  1:176288/1‑300 (MQ=255)
                                                                                                                                                                                                        cactctttccctacacgtcgctcttccgatcTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTg                                                                                                     <  1:61286/271‑1 (MQ=255)
                                                                                                                                                                                                                     tGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTgg                                                                                        >  2:121907/1‑301 (MQ=255)
                                                                                                                                                                                                                           ggATATTTTTACTCCATGTCGTTGGTGCTGCGAGTGTCGCATATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGaa                                                                                  <  2:15312/301‑1 (MQ=255)
                                                                                                                                                                                                                             gagttcagacgtgtgctcttccgatctgcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATg                                                                                <  2:159622/276‑1 (MQ=255)
                                                                                                                                                                                                                                       tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCAGCTTTATCAATCACAAAGGTTTTGCCCACGTTACCTGGGTGAAGTTGTGGCAAGAAAGAGTTTGCAGa                                                                      >  1:542238/1‑301 (MQ=255)
                                                                                                                                                                                                                                        ccATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACCGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGaa                                                                     >  2:81227/1‑301 (MQ=255)
                                                                                                                                                                                                                                                  ggTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGt                                                           >  1:441682/1‑301 (MQ=255)
                                                                                                                                                                                                                                                       tgcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGGGAGGTTGTGGCAAGACTGAATTTGCAGAAAGCAGAACGGTGAtt                                                      >  1:379123/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        gcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAgatcggaagagcgtcgtgtagggaaa                                                     >  1:159471/1‑276 (MQ=255)
                                                                                                                                                                                                                                                              ttCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAg                                               <  2:519195/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                cgccgttAAACATACCCAATAGAAAACCCTCACCGTCAGTCTGCGAGTGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCt                                             <  1:380521/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aaaCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCg                                      <  1:121783/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aaCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCga                                     >  2:307067/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                         aCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCgag                                    <  2:640854/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aaCCCAAAAGAAAACCCTCACCGTCAGGCGGTGAGGGTTTAATTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAAGTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATc                                 <  1:388882/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ccctctcttccgttctccttCAGGCTGCGAGGTTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGaa                             <  1:384142/282‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGaa                             >  2:274904/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                    aGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTc                         >  2:290635/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                        aaCCCTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTCATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAACGCGGAGCGTTGATGACGTCAGCTTTGATCGACACCTGAATTTCAATc                     >  1:629507/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                         tcccTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCa                    <  2:360161/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          cccTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGCTTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGACTGAGTTTGCAGAAAGCGGAGCGTTGCTGGCGACAGCTTTGCTAGAGCTCTGAATTtcaatcaa                   >  1:490689/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          cccTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTAAGTCaa                   >  1:359677/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          cccTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAAATCTGAATTTCAGTCaa                   >  1:169511/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                           ccTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGTATTAGCGTGTTTTGTAAATCAGATGATTATTAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGACAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCaaa                  >  2:152597/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                ccGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGaaagatcggaagagcaca             >  2:384638/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                                                           cgAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGTGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCt  <  2:485031/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                            
CAGCGTGCGCACACGCACTTCGCCATCGCTAAGATTGCGGAATTTAACCCGGCTGTTCATTGTCACCACGTCGTGTGGCATCTCTTCTGGCGAACACATTTGGGCGCGATCCAACTCTGCGTTTAACGCGTCGGCGATTGGCAAACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTAAATCAGATGATTAATAACCCGCTTTATCAATCACAAAGGTTTTGCCACAGTTACCTGGGTGAGGTTGTGGCAAGAATGAGTTTGCAGAAAGCGGAGCGTTGATGGCGTCAGCTTTGATCGAGATCTGAATTTCAGTCAAATACGCAGGGTTACCCT  >  NC_000913/643743‑644342

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