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,0770CT86.1% 84.6 / 4.3 37intergenic (+459/+492)iap/cas2alkaline phosphatase isozyme conversion protein/CRISPR‑associated endoribonuclease Cas2
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (13/18);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

CCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCC  >  NC_000913/2877831‑2878322
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ccccGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTAGCGCGGGGACATCAAGCTGGCTGGCAATCTCTTTCGGGGGGGGTCCGATTTTTCTCCGCt                                                                                                                                                                                                                                                   >  2:387663/1‑251 (MQ=255)
              ggAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTAGCTGGCGCGGGGAAct                                                                                                                                                                                                                                      >  1:333737/1‑250 (MQ=255)
                    cGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGt                                                                                                                                                                                                                                 <  1:246418/249‑1 (MQ=255)
                        gTCCATCATGCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGt                                                                                                                                                                                                                           <  1:179970/251‑1 (MQ=255)
                              catTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGCTAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGATGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGGGAGGCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTcc                                                                                                                                                                                                                     >  2:57402/1‑251 (MQ=255)
                                         tGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATaa                                                                                                                                                                                                           >  1:229822/1‑250 (MQ=255)
                                          gTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATaaa                                                                                                                                                                                                          <  1:57402/250‑1 (MQ=255)
                                              gAACTCCCGGTTTATCCCCGCTAGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCCTGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGAGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCGGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCTCGGGGAACTCTTATTTCCGTATCTCCGGATTtatacaactg                                                                                                                                                                                                      >  2:78112/1‑244 (MQ=255)
                                               aaCTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGCTAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGAGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGAc                                                                                                                                                                                                    >  1:31963/1‑251 (MQ=255)
                                                   cccGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGGTTCCGTATCTCCGGATTTATAAAGCTGACGGt                                                                                                                                                                                                 <  1:396331/250‑1 (MQ=255)
                                                   cccGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCGCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGtt                                                                                                                                                                                                >  2:137700/1‑251 (MQ=255)
                                                              cccGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTg                                                                                                                                                                                     <  1:137700/251‑1 (MQ=255)
                                                                  cTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGGCCGGTTTATCCTCGCTTGCCCGGGGAACACTCTTTTCCGTCTATCCGGATTTATAAAGCTAACGGTTTATCCACGCTGgcgc                                                                                                                                                                                 >  2:281943/1‑251 (MQ=255)
                                                                           ggAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTc                                                                                                                                                                        <  1:43711/251‑1 (MQ=255)
                                                                                   cGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAggcgg                                                                                                                                                                >  1:360270/1‑251 (MQ=255)
                                                                                   cGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAggcgg                                                                                                                                                                <  1:76121/251‑1 (MQ=255)
                                                                                         aTAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGGACTCTAGGTTCCGTATCTCCGGATTTATAAAGCTGACGGTCTATCCCCGCTGGCGCGGGGACCTCGCAGGCGGCGAcgc                                                                                                                                                          >  2:273899/1‑251 (MQ=255)
                                                                                                                            ccgggggCGCGGGGAACTCTGGGCGGCGTGCCTTGCAGCCAGCGGCAGCCGCGGTTTATCCCCGCTGACGCGGGGAACTCAAGCTGGCTGGCAAGCTCCTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTg                                                                                                                       <  2:129244/246‑1 (MQ=255)
                                                                                                                                 ggCGCGGGGAACTCTGGGCGGCTTGCCGTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCGGGCAATCTCGTTCGGTGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGcgcg                                                                                                                  <  2:301316/251‑1 (MQ=255)
                                                                                                                                          caCTCTGGGCGGTTTGCCTGGAAGCCGTCGCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTGTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTcg                                                                                                         <  2:21009/250‑1 (MQ=255)
                                                                                                                                              ctGGGTGGCTTGCCTTTCAGCCAGCTCCAGCTGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGTGGTGAGGCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTAGCTCCGGATGTATAAAGCTGACGGTTTAGCCCCGCTGGCGCGGGGAACTCTCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTc                                                                                                          <  2:388014/246‑1 (MQ=255)
                                                                                                                                                      cttgccttgcAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGtt                                                                                                                                  >  1:353974/1‑214 (MQ=255)
                                                                                                                                                      cttgccttgcAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGGGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGtt                                                                                                                                  <  2:353974/214‑1 (MQ=255)
                                                                                                                                                            ttgcAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTcc                                                                                       >  1:83518/1‑251 (MQ=255)
                                                                                                                                                             tgcAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTc                                                                                        <  1:403537/249‑1 (MQ=255)
                                                                                                                                                                  ccagcTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCTATTCGCGGTTTAGCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGAccc                                                                                 <  2:174220/251‑1 (MQ=255)
                                                                                                                                                                        cctgcagcGGTTTATCCCCGTTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTGCGGGGTGAGTCCGGTTTAGCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATc                                                                             <  1:281943/246‑1 (MQ=255)
                                                                                                                                                                                              tgcgcGGGTATCTCAAGCGGGCGGGCAGTCTCTTTCGGGGGGAGTTAGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTGTCCGTATCTCTGGATTTATAACGCTGACGGTTTATCCCCGCTGGCGCGGGGAGCTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGATCTCGCGACCGCTGAGAAATTCCAGACCCGAGCCAAATGTTTTAGACCCGCTGgc                                                       <  2:297192/248‑1 (MQ=255)
                                                                                                                                                                                                  cgGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGTGGAACTCGCAGGCGGCGACGCGCAGGGTATgcgc                                                                                                                                             >  1:75320/1‑159 (MQ=255)
                                                                                                                                                                                                  cgGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGTGGAACTCGCAGGCGGCGACGCGCAGGGTATgcgc                                                                                                                                             <  2:75320/159‑1 (MQ=255)
                                                                                                                                                                                                   ggggAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATACAAACGGTTTATCCCCGCTGGCGCgggg                                                 >  1:439613/1‑250 (MQ=255)
                                                                                                                                                                                                     ggAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGCCGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAACTTCCAGGCCCGATCCAAACGGTTTATCCCCGCTGTCGCGGGGAAc                                              >  2:310860/1‑251 (MQ=255)
                                                                                                                                                                                                                             cttgcGGGGTGAGTCCGGTTTATCCTCGCTTGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCATGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCTCTCAGAAATTCCAGCCCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACa                          <  1:26815/243‑1 (MQ=255)
                                                                                                                                                                                                                               ttCGGGGTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGGTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGCTCCAAACGGTTTATCCCCGCTGGCGCGGGGACCTCTCAACATTATCAATTACCACCGAc                    >  2:353558/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cGGGGTGATTCCGGGTTAACTCCGCGGGCGCGGGGAAGTCTATTTTCCGTATCGCCGGATGTAGAAAGCTGACGGTTGATCCCCGCTGGTGCGGGGAACTCTCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGTAACTCGCGACCGCTCAGAAATGCCAGTCCCGTTCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCCCAACATTTTCAATTTCAACCGACa                   <  2:229421/250‑1 (MQ=255)
                                                                                                                                                                                                                                    ggTGAGTCCGGTTTATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGa               <  1:198506/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 tATCCTCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATACCCGCTGGCGCGGGGAAATCGCAGGCGGCGACGCGAAGGGTATGAGCGGTTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTAAGAAATTCCAGACCCGATACAAACGGTTTATCCCCGCTGGCGCGGGGAAGTCTCAACAATATCAATTACAAACGACAGGGAGACAGGTTTATcc  >  1:54779/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CCCCGCTAACGCGGGGAACTCGTAGTCCATCATTCCACCTATGTCTGAACTCCCGGTTTATCCCCGCTGGCGCGGGGAACTCCCGGGGGATAATGTTTACGGTCATGCGCCCCCCGGTTTATCCCCGCTGGCGCGGGGAACTCTGGGCGGCTTGCCTTGCAGCCAGCTCCAGCAGCGGTTTATCCCCGCTGGCGCGGGGAACTCAAGCTGGCTGGCAATCTCTTTCGGGGTGAGTCCGGTTTATCCCCGCTGGCGCGGGGAACTCTAGTTTCCGTATCTCCGGATTTATAAAGCTGACGGTTTATCCCCGCTGGCGCGGGGAACTCGCAGGCGGCGACGCGCAGGGTATGCGCGATTCGCGGTTTATCCCCGCTGGCGCGGGGAACTCGCGACCGCTCAGAAATTCCAGACCCGATCCAAACGGTTTATCCCCGCTGGCGCGGGGAACTCTCAACATTATCAATTACAACCGACAGGGAGCCCGGTTTATCC  >  NC_000913/2877831‑2878322

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