Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,843,622 C→T 6.6% R353H (CGC→CAC)  PP_3393 ← CAIB/BAIF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,843,6220CT 6.6% 178.9 / 2.4 61R353H (CGC→CAC) PP_3393CAIB/BAIF family protein
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

CATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGGGCGGCGGATCGGCGCGAAAGGAATGCCGTT  >  NC_002947/3843474‑3843770
                                                                                                                                                    |                                                                                                                                                     
cATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTACAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGgc                                                                                                                                                       >  1:49189/1‑149 (MQ=255)
      aaGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCAGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGc                                                                                                                                                 >  1:118256/1‑149 (MQ=255)
            tcctcAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCAGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTAACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTc                                                                                                                                             >  2:224639/1‑147 (MQ=255)
             cctcAACCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGACGTGGCTGCAGGCGCTGGCCGTCCAg                                                                                                                                          >  3:42374/1‑149 (MQ=255)
              ctcAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGt                                                                                                                                         <  8:295064/149‑1 (MQ=255)
              ctcAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGt                                                                                                                                         <  3:835/149‑1 (MQ=255)
              ctcAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGt                                                                                                                                         <  2:49189/149‑1 (MQ=255)
                     cgcgTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCAGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCg                                                                                                                                  <  2:118256/149‑1 (MQ=255)
                      gcgTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCgg                                                                                                                                 <  3:43137/149‑1 (MQ=255)
                      gcgTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCgg                                                                                                                                 <  2:107784/149‑1 (MQ=255)
                      gcgTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCgg                                                                                                                                 <  3:162911/149‑1 (MQ=255)
                        gtcgtcTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGTGCTGGCCGTCCAGTGATAGCGGc                                                                                                                                >  6:279010/1‑148 (MQ=255)
                         tcgtcTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCaa                                                                                                                              <  1:109102/149‑1 (MQ=255)
                                ttttCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACgg                                                                                                                        >  3:192919/1‑148 (MQ=255)
                                  ttCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCAt                                                                                                                     >  6:130749/1‑149 (MQ=255)
                                       cACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGgtg                                                                                                                >  8:22838/1‑149 (MQ=255)
                                        aCGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGgtgt                                                                                                               >  3:275505/1‑149 (MQ=255)
                                          gcccgccGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGgtgtgt                                                                                                             <  5:130749/149‑1 (MQ=255)
                                               ccgcACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTg                                                                                                        >  8:257535/1‑149 (MQ=255)
                                                   acacaGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGAGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGt                                                                                                    >  7:98166/1‑149 (MQ=255)
                                                       aGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTTTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGAGCTGGCCGTCCAGTGATAGCGGCAATaa                                                                                                                           <  5:150486/122‑1 (MQ=255)
                                                       aGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTTATAGCGGAAATaa                                                                                                                           >  6:150486/1‑122 (MQ=255)
                                                       aGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTc                                                                                                <  4:42374/149‑1 (MQ=255)
                                                         cTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGTGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCa                                                                                              <  5:279010/149‑1 (MQ=255)
                                                         cTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTc                                                                                                >  2:18559/1‑147 (MQ=255)
                                                            ggCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGtt                                                                                           >  2:246153/1‑149 (MQ=255)
                                                             gCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGtcg                                                                                          <  8:155160/149‑3 (MQ=255)
                                                             gCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTg                                                                                          >  6:9247/1‑149 (MQ=255)
                                                              cAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTg                                                                                          >  7:211203/1‑148 (MQ=255)
                                                               aaTGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCa                                                                                                                      >  8:124040/1‑119 (MQ=255)
                                                                        gTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCa                                                                               >  2:247631/1‑149 (MQ=255)
                                                                           gCTGTATCCTAGTTCACGCATCTCTTGGCGTGTGTGCTCACCATTCCTAGCGATGGCACGCCGTGGCTGCAGGCGCTGTCCGTCGAGTGATAGCGGCAATAACGGGATCGGTGTGTGGCTGCGGTCTTCCAGTTGCAGCTCGGCCATg                                                                             <  4:275505/148‑1 (MQ=255)
                                                                           gCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGTGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGc                                                                            >  3:6181/1‑149 (MQ=255)
                                                                                   ccAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGctc                                                                   >  6:234725/1‑149 (MQ=255)
                                                                                    cAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGctct                                                                  <  8:98166/149‑1 (MQ=255)
                                                                                      gTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGt                                                                <  8:211203/149‑1 (MQ=255)
                                                                                        tCACGCATCACTTGGCGGGAGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAg                                                              >  2:116908/1‑149 (MQ=255)
                                                                                                 aCTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTgggg                                                     >  8:43033/1‑149 (MQ=255)
                                                                                                   ttGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGct                                                                    <  5:190752/132‑1 (MQ=255)
                                                                                                   ttGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGct                                                                    >  6:190752/1‑132 (MQ=255)
                                                                                                   ttGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGt                                                   >  6:86053/1‑149 (MQ=255)
                                                                                                      gCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGTGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCt                                                <  4:6181/149‑1 (MQ=255)
                                                                                                       cGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTc                                               <  5:264111/149‑1 (MQ=255)
                                                                                                             tgCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGaaaa                                          <  7:22838/148‑1 (MQ=255)
                                                                                                               cTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGct                                       >  5:196221/1‑149 (MQ=255)
                                                                                                                  aCCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGctct                                     <  1:18559/148‑1 (MQ=255)
                                                                                                                  aCCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCtt                                    <  6:196221/149‑1 (MQ=255)
                                                                                                                  aCCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCtt                                    >  3:105144/1‑149 (MQ=255)
                                                                                                                    cAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCtt                                    >  1:52259/1‑147 (MQ=255)
                                                                                                                    cAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCtt                                    >  2:254344/1‑147 (MQ=255)
                                                                                                                         cTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGGGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGggc                              <  4:192919/148‑1 (MQ=255)
                                                                                                                          tCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGc                                        >  5:218370/1‑137 (MQ=255)
                                                                                                                          tCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGc                                        <  6:218370/137‑1 (MQ=255)
                                                                                                                            gcgATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCTGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGggcggc                           <  1:254344/148‑1 (MQ=255)
                                                                                                                            gcgATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGggcggc                           <  7:257535/148‑1 (MQ=255)
                                                                                                                              gATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGGGCGGCGGa                        <  5:234725/149‑1 (MQ=255)
                                                                                                                              gATCGCCCGCCGTGGCTGCAGGCGCTGACCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGGGCGGCGGa                        >  1:8737/1‑149 (MQ=255)
                                                                                                                                   gccgccgTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGGGCGGCGGATCGgc                   <  4:105144/148‑1 (MQ=255)
                                                                                                                                   cccgccGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCGCGGCCATGCCTGCCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGGGCGGCGGATCgg                    <  1:246153/149‑1 (MQ=255)
                                                                                                                                    ccgccgTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTGGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGGGCGGCGGATCGgcg                  >  7:200418/1‑149 (MQ=255)
                                                                                                                                                 aGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTc                                                                                                >  3:183546/1‑59 (MQ=255)
                                                                                                                                                 aGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTc                                                                                                <  4:183546/59‑1 (MQ=255)
                                                                                                                                                    cgcTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTGGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGGGCGGCGGATCGGCGCGAAAGGAATGCCGtt  <  8:200418/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATCTCAAGGCCTCCTCAAGCCGCGTCGTCTGTTTTCAGCACGCCCGCCGCACACAGCTGGGCAATGTGCGCGTCGCTGTAACCCAGTTCACGCATCACTTGGCGGGTGTGCTCACCAATCCTCGCGATCGCCCGCCGTGGCTGCAGGCGCTGGCCGTCCAGTGATAGCGGCAATAACGGCATCGGTGTGTGGCTGCCGTCTTCCAGTTGCAGCTCGGCCATGCC‑ACCGCTCTGTAGCAGGTGGGGGTCCTCGAAAAGCTCTTCAGGGCGGCGGATCGGCGCGAAAGGAATGCCGTT  >  NC_002947/3843474‑3843770

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑