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

CCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCC  >  NC_000913/2877833‑2878322
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ccGCTAACGCGGGGAACTCGTAGTCCATCATTACACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCACGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCCGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCTGGGGGCGCCCGGTTTATCCTCGCTgg                                                                                                                                                                                                                                                 >  2:400040/1‑251 (MQ=255)
       cgcgGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCgggg                                                                                                                                                                                                                                          >  1:111411/1‑251 (MQ=255)
           gggAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAAct                                                                                                                                                                                                                                      <  1:360618/251‑1 (MQ=255)
                  cGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGCACTCAAGCAGGCTGGCAATCTCTTTCGGGGGGAGTCCGGTTTATCCTCGCTGGCGCGGGTAACTCTAGttt                                                                                                                                                                                                                               >  2:647092/1‑251 (MQ=255)
                  cGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGt                                                                                                                                                                                                                                 <  1:713141/249‑1 (MQ=255)
                      gTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGAGGGGAACTCAAGCTGGCTGGCCATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGt                                                                                                                                                                                                                           >  2:567306/1‑251 (MQ=255)
                                ccacctatggcgtatcTCCCGGTTTAGCCCCGCTGGCTCGGGGAACTCCCGGTGGACAGTGTTTACGGTCAGTCGCACCCCGGTTTATCCCCGCTGGCGCGTGGATCGCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCCATCTCTTTCGGGGTGAGGCCGGTTTATCCTCGCTGTCGCGGTGAACTCTAGTTTCCGGATCTCCGGAt                                                                                                                                                                                                                 <  2:554075/236‑1 (MQ=255)
                                ccACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGa                                                                                                                                                                                                                  <  1:567138/250‑1 (MQ=255)
                                  aCCTATGTCAGAACTCCCGGTTTATCCCCGCTGGAGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATACGCCCCCCGGCTTATCCCCGCTGGGGCGGGGAACTATGGGCGGCTTGCCTTGCAGCCAGCTCCCGCAGCGGTTTATCCCCGC‑GGGCGCGGGGAACTCAAGCTGGCTGCCAATATCTTTCGGGGTGAATCCGGTAGATCCACGATGGCGCGAGGAACGCTAGTTTACGGAtctcagtattt                                                                                                                                                                                                               >  2:117034/1‑244 (MQ=255)
                                      atggcgtaacggccgggttacacccGCTGGAGCGGGGAACCCCCGGTGGATAATGGTGACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGCGCGGCTGGCCGTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTTGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCTCGGGGAACTCTAGTTTCCGTATCTCCGGATTtata                                                                                                                                                                                                            <  2:58216/228‑1 (MQ=255)
                                        gTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATaaa                                                                                                                                                                                                          <  1:611817/250‑1 (MQ=255)
                                                  ccGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTTGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGttt                                                                                                                                                                                               <  1:567306/251‑1 (MQ=255)
                                                      tttATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATcc                                                                                                                                                                                           >  1:234509/1‑251 (MQ=255)
                                                           ccccGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTg                                                                                                                                                                                                                                                  <  1:741279/191‑1 (MQ=255)
                                                            cccacgcgcgcTGGGTACTCCGGGGGGATAATTTTGACGGTCAGGCCCCCCCCGTTTTATCCCCGCTGGCGCGGGGAACTCGGGTCGGCTTGCCTTGCAGCCAGCGCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATTTCCTGATTTATAAAGCTGACGGTTTATCCCCGCTg                                                                                                                                                                                     <  2:526514/244‑1 (MQ=255)
                                                                  ggCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCgg                                                                                                                                                                               >  1:41839/1‑251 (MQ=255)
                                                                   gcgcGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCggg                                                                                                                                                                              >  1:371367/1‑251 (MQ=255)
                                                                        gggaacgaccGGGGGACAAAGTTTACGGTCATGCGCACCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGAAGCGGTGTATCCCCGC‑TGGCGCGGGGAACGCAAGATGGCTGGCAAGCTTTTTCGGGGGTAGCCCGGTGTATCCGCGCTGGCGCTGGGGACTATAGTGTCGGGATCTACGGGTTTAAAATGCTGACGGTTTATcaccaattgcgcggggc                                                                                                                                                                            >  2:719232/9‑232 (MQ=255)
                                                                         ggtaCTCCCGGTGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTc                                                                                                                                                                        <  2:241858/248‑1 (MQ=255)
                                                                           aaCTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGc                                                                                                                                                                      >  1:369213/1‑251 (MQ=255)
                                                                                ccGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAggcg                                                                                                                                                                 <  1:689480/251‑1 (MQ=255)
                                                                                  gggggATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCCGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGGTTATCCTCGCTGGCGCGGGGAACTCTAGATTCCGTATCTCCGGATTTATAAAGCGGACGGTCGATCCCCGCTGGCCCGGCGACCTCGCatcggcc                                                                                                                                                               >  2:185645/1‑245 (MQ=255)
                                                                                     ggATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGTGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TTGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGAc                                                                                                                                                            <  1:326117/251‑1 (MQ=255)
                                                                                                 cGGTCATGCCCCCCCCGGTTGATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCCTGCAGCCAGCTCCAGCTGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATg                                                                                                                                                <  2:257901/251‑1 (MQ=255)
                                                                                                        gcgcCCCCCGGTTTATCCACGCTGGCGAGGGGAACTCTGGGCGGCTTGCCTCGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCTATGGCGCGGGGAACTCTCGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGTCACTCGCAGGCGGCGACCCACAGGGTATGCGCGAtt                                                                                                                                         >  2:338563/1‑251 (MQ=255)
                                                                                                                cGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGCCGCGCAGGGTATGCGCGATTAGCGGttt                                                                                                                                 >  1:469704/1‑251 (MQ=255)
                                                                                                                  gTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTAt                                                                                                                               >  1:56631/1‑251 (MQ=255)
                                                                                                                     tATCCCCGCTGGTGCGGGGAACTCTTGGCGGCGTGCCTTGCAGCCATCTCCCTCATCGGTTTATCCCCGC‑TGGCTCTGGGAACTTAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGTGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCTGGGATCTCTCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATcc                                                                                                                             <  2:717971/250‑1 (MQ=255)
                                                                                                                     tATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATcc                                                                                                                             <  1:400040/250‑1 (MQ=255)
                                                                                                                      aTCCCCGCCGGCGCGGGGAACTCTGGGCGGCTTGCCTGGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATcccc                                                                                                                           <  2:444248/251‑1 (MQ=255)
                                                                                                                          ccGCTGGCTCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTg                                                                                                                       <  1:691174/251‑1 (MQ=255)
                                                                                                                                  gcgGGGAACTCTGGGATGCTTGCCTGGCAGGCAGCTCCAGCATCGGTTTATCCCCGCGTGG‑GCGGGGAACTCATGCTGGCTGGCAACCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAATTCTAGTTTCCGTATCTCCGGATTTATAAAGCGGACTGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACTCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCggg                                                                                                                <  1:485996/250‑1 (MQ=255)
                                                                                                                                       gAACTCTGGGCGGCTTGCCTTGCAGCCCGCGCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTGTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTc                                                                                                          <  1:5412/251‑1 (MQ=255)
                                                                                                                                           cctgtgtggattgccgttcatccagcTCCAGCAGCGGTTTATCCCCGC‑TTACTCGGGGAACTCAAGCGGTCTGGCAATCGCTTGCGGGGTGAGCTCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGCATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCTGTTTATCCCCGCTGGCGCGGGTAACTcgcg                                                                                                       <  2:469704/229‑1 (MQ=255)
                                                                                                                                               ggcggcTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGAcctc                                                                                                  >  1:162673/1‑249 (MQ=255)
                                                                                                                                                                                gtgaTCCCCGC‑TGGCGTGGGGAACTCAACCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTAGGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTa                                                                  <  2:371367/247‑1 (MQ=255)
                                                                                                                                                                                        cGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTg                                                         >  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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGC‑TGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCC  >  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: ‑