Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,878,077 C→T intergenic (+459/+492) iap → / ← cas2 alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,878,0770CT97.6% 102.2 / NA 43intergenic (+459/+492)iap/cas2alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (16/24);  minor base G (1/0);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCC  >  NC_000913/2877833‑2878322
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ccGCTAACGCGGGGAACTCGTAGTCCATCATTACACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCACGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCCGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCTGGGGGCGCCCGGTTTATCCTCGCTgg                                                                                                                                                                                                                                                 >  2:400040/1‑251 (MQ=255)
       cgcgGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCgggg                                                                                                                                                                                                                                          >  1:111411/1‑251 (MQ=255)
           gggAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAAct                                                                                                                                                                                                                                      <  1:360618/251‑1 (MQ=255)
                  cGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGCACTCAAGCAGGCTGGCAATCTCTTTCGGGGGGAGTCCGGTTTATCCTCGCTGGCGCGGGTAACTCTAGttt                                                                                                                                                                                                                               >  2:647092/1‑251 (MQ=255)
                  cGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGt                                                                                                                                                                                                                                 <  1:713141/249‑1 (MQ=255)
                      gTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGAGGGGAACTCAAGCTGGCTGGCCATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGt                                                                                                                                                                                                                           >  2:567306/1‑251 (MQ=255)
                                ccacctatggcgtatcTCCCGGTTTAGCCCCGCTGGCTCGGGGAACTCCCGGTGGACAGTGTTTACGGTCAGTCGCACCCCGGTTTATCCCCGCTGGCGCGTGGATCGCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCCATCTCTTTCGGGGTGAGGCCGGTTTATCCTCGCTGTCGCGGTGAACTCTAGTTTCCGGATCTCCGGAt                                                                                                                                                                                                                 <  2:554075/236‑1 (MQ=255)
                                ccACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGa                                                                                                                                                                                                                  <  1:567138/250‑1 (MQ=255)
                                  aCCTATGTCAGAACTCCCGGTTTATCCCCGCTGGAGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATACGCCCCCCGGCTTATCCCCGCTGGGGCGGGGAACTATGGGCGGCTTGCCTTGCAGCCAGCTCCCGCAGCGGTTTATCCCCGCGGGCGCGGGGAACTCAAGCTGGCTGCCAATATCTTTCGGGGTGAATCCGGTAGATCCACGATGGCGCGAGGAACGCTAGTTTACGGAtctcagtattt                                                                                                                                                                                                               >  2:117034/1‑244 (MQ=255)
                                      atggcgtaacggccgggttacacccGCTGGAGCGGGGAACCCCCGGTGGATAATGGTGACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGCGGCTGGCCGTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTTGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCTCGGGGAACTCTAGTTTCCGTATCTCCGGATTtata                                                                                                                                                                                                            <  2:58216/228‑1 (MQ=255)
                                        gTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATaaa                                                                                                                                                                                                          <  1:611817/250‑1 (MQ=255)
                                                  ccGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTTGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGttt                                                                                                                                                                                               <  1:567306/251‑1 (MQ=255)
                                                      tttATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATcc                                                                                                                                                                                           >  1:234509/1‑251 (MQ=255)
                                                           ccccGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTg                                                                                                                                                                                                                                                  <  1:741279/191‑1 (MQ=255)
                                                            cccacgcgcgcTGGGTACTCCGGGGGGATAATTTTGACGGTCAGGCCCCCCCCGTTTTATCCCCGCTGGCGCGGGGAACTCGGGTCGGCTTGCCTTGCAGCCAGCGCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATTTCCTGATTTATAAAGCTGACGGTTTATCCCCGCTg                                                                                                                                                                                     <  2:526514/244‑1 (MQ=255)
                                                                  ggCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCgg                                                                                                                                                                               >  1:41839/1‑251 (MQ=255)
                                                                   gcgcGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCggg                                                                                                                                                                              >  1:371367/1‑251 (MQ=255)
                                                                        gggaacgaccGGGGGACAAAGTTTACGGTCATGCGCACCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGAAGCGGTGTATCCCCGCTGGCGCGGGGAACGCAAGATGGCTGGCAAGCTTTTTCGGGGGTAGCCCGGTGTATCCGCGCTGGCGCTGGGGACTATAGTGTCGGGATCTACGGGTTTAAAATGCTGACGGTTTATcaccaattgcgcggggc                                                                                                                                                                            >  2:719232/9‑232 (MQ=255)
                                                                         ggtaCTCCCGGTGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTc                                                                                                                                                                        <  2:241858/248‑1 (MQ=255)
                                                                           aaCTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGc                                                                                                                                                                      >  1:369213/1‑251 (MQ=255)
                                                                                ccGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAggcg                                                                                                                                                                 <  1:689480/251‑1 (MQ=255)
                                                                                  gggggATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCCGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGGTTATCCTCGCTGGCGCGGGGAACTCTAGATTCCGTATCTCCGGATTTATAAAGCGGACGGTCGATCCCCGCTGGCCCGGCGACCTCGCatcggcc                                                                                                                                                               >  2:185645/1‑245 (MQ=255)
                                                                                     ggATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGTGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTTGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGAc                                                                                                                                                            <  1:326117/251‑1 (MQ=255)
                                                                                                 cGGTCATGCCCCCCCCGGTTGATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCCTGCAGCCAGCTCCAGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATg                                                                                                                                                <  2:257901/251‑1 (MQ=255)
                                                                                                        gcgcCCCCCGGTTTATCCACGCTGGCGAGGGGAACTCTGGGCGGCTTGCCTCGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCTATGGCGCGGGGAACTCTCGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGTCACTCGCAGGCGGCGACCCACAGGGTATGCGCGAtt                                                                                                                                         >  2:338563/1‑251 (MQ=255)
                                                                                                                cGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGCCGCGCAGGGTATGCGCGATTAGCGGttt                                                                                                                                 >  1:469704/1‑251 (MQ=255)
                                                                                                                  gTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTAt                                                                                                                               >  1:56631/1‑251 (MQ=255)
                                                                                                                     tATCCCCGCTGGTGCGGGGAACTCTTGGCGGCGTGCCTTGCAGCCATCTCCCTCATCGGTTTATCCCCGCTGGCTCTGGGAACTTAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGTGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCTGGGATCTCTCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATcc                                                                                                                             <  2:717971/250‑1 (MQ=255)
                                                                                                                     tATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATcc                                                                                                                             <  1:400040/250‑1 (MQ=255)
                                                                                                                      aTCCCCGCCGGCGCGGGGAACTCTGGGCGGCTTGCCTGGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATcccc                                                                                                                           <  2:444248/251‑1 (MQ=255)
                                                                                                                          ccGCTGGCTCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTg                                                                                                                       <  1:691174/251‑1 (MQ=255)
                                                                                                                                  gcgGGGAACTCTGGGATGCTTGCCTGGCAGGCAGCTCCAGCATCGGTTTATCCCCGCGTGGGCGGGGAACTCATGCTGGCTGGCAACCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAATTCTAGTTTCCGTATCTCCGGATTTATAAAGCGGACTGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACTCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCggg                                                                                                                <  1:485996/250‑1 (MQ=255)
                                                                                                                                       gAACTCTGGGCGGCTTGCCTTGCAGCCCGCGCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTGTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTc                                                                                                          <  1:5412/251‑1 (MQ=255)
                                                                                                                                           cctgtgtggattgccgttcatccagcTCCAGCAGCGGTTTATCCCCGCTTACTCGGGGAACTCAAGCGGTCTGGCAATCGCTTGCGGGGTGAGCTCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGCATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCTGTTTATCCCCGCTGGCGCGGGTAACTcgcg                                                                                                       <  2:469704/229‑1 (MQ=255)
                                                                                                                                               ggcggcTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGAcctc                                                                                                  >  1:162673/1‑249 (MQ=255)
                                                                                                                                                                                gtgaTCCCCGCTGGCGTGGGGAACTCAACCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTAGGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTa                                                                  <  2:371367/247‑1 (MQ=255)
                                                                                                                                                                                        cGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTg                                                         >  1:110815/1‑251 (MQ=255)
                                                                                                                                                                                              gcgcTGGGAACTCAACCTGGCTGGCAAGCTCCTTCGGGGTGACGCCGGTGTATCCTCGCTGGCGCGGGGAACTCTAGGTTCCGTCTCTCCGAATTTATAAAGCTGACGGTTTCTCCCCGCTGGCGCGGGTAACTCGCAGGCGGCGGCGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGcgcg                                                    <  2:279872/251‑1 (MQ=255)
                                                                                                                                                                                               cgcgGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGCCGCGGGGAACTCTAGTTTCCGTATCTCCGGGTTTATAACGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCACGATTCGCCGTGTTTACCCGCCGGCGAGGGGTACTCGACAACGCTAAGCAATTCCAGACCCCATCCAGACGGTTTATcaccactgcatctg                                                   >  2:747847/1‑238 (MQ=255)
                                                                                                                                                                                                               gggctggcagtcgctcgcGGGGAGAGTCAGGTTTACTCTCGCTGGCGCGGGGAAATCTAGTTTCCGGACCTGCGGATTCAGACGGCTAGCGGTTTAGCCGCGCTTGCGCGGGTGACTCGCAGGCGGCGACGCGCAGGGTAAGCTCGATTCGCGGTTTATCCCCGCTGGCGCGGGGACCTCGCGCCCGCGCAGAAATTCCAGACTCGATCCAAACGGTTTAGCCCTGCTGGCGCGGGTAACTCTCAACAtca                                   <  2:160861/234‑3 (MQ=255)
                                                                                                                                                                                                                ggctggcAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTAt                                  <  1:506994/251‑1 (MQ=255)
                                                                                                                                                                                                                         tcgcttgcgcggTGAGTCCGGTTTAGCCTCGTTGCCGTGGGGAACTCTAGGTTCCGTATCTCCGGATCTATAAAGCTGTGGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTAGCCCCGCTGGCGCTGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCCCAACATTATCAATTACa                          <  2:347476/240‑1 (MQ=255)
                                                                                                                                                                                                                                                tATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGCTTATCCCCGCTGGCGCGGCTAACCCGCGCCCGATACTAAATTCCAGACCCGCTCACAACGGCTTATCCCCACAGGCGCGGGGAACCCTCAACAATATCAATAACAACCGACacacatcccgggttcacc  >  2:498139/1‑234 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCC  >  NC_000913/2877833‑2878322

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