Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,650,718 T→A 6.1% Q35H (CAA→CAT)  PP_3216 ← cysteine dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,650,7180TA 6.1% 193.7 / 2.7 65Q35H (CAA→CAT) PP_3216cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base A (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CCCCGGGCCCCAGACGAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAACCCTAGCAAACCT  >  NC_002947/3650574‑3650866
                                                                                                                                                |                                                                                                                                                    
ccccGGGCCCCAGACGAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTGATGGGCGACCAGATCGGCCAGCAGGCCTTg                                                                                                                                                    <  1:39846/147‑1 (MQ=255)
ccccGGGCCCCAGACGAAGCTCACCACCGAGAAACGTTGGCCCGAGTCGCAGTGCAGCAGGTATTGCTCATAGCGATGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGc                                                                                                                                                  >  8:340882/1‑149 (MQ=255)
  ccGGGCCCCAGACGAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCt                                                                                                                                                <  8:149668/149‑1 (MQ=255)
    gggCCCCAGACGAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGc                                                                                                                                              <  5:212428/149‑1 (MQ=255)
      gCCCCAGACGAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCAAGGCCTGctg                                                                                                                                            <  5:268898/149‑3 (MQ=255)
       caccAGACGAAGCTCACCACCGAGAATCGTTGGCTCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCAAGGCCTGctgg                                                                                                                                           <  4:79600/147‑4 (MQ=255)
          cAGACGAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTc                                                                                                                                           >  6:163175/1‑146 (MQ=255)
             aCGAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAg                                                                                                                                     >  8:213369/1‑149 (MQ=255)
               gAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAggg                                                                                                                                   >  2:111687/1‑149 (MQ=255)
                aaGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGt                                                                                                                                  <  1:111687/149‑1 (MQ=255)
                 aGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTg                                                                                                                                 <  7:354790/149‑1 (MQ=255)
                     caccacCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGct                                                                                                                                             <  7:196432/133‑2 (MQ=255)
                     caccacCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGct                                                                                                                                             >  8:196432/1‑132 (MQ=255)
                       ccaccGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCg                                                                                                                           >  4:130213/1‑149 (MQ=255)
                        caccGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGt                                                                                                                          >  6:380203/1‑149 (MQ=255)
                          ccGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCa                                                                                                                        <  5:163175/149‑1 (MQ=255)
                                 aCGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGtt                                                                                                                 >  6:118372/1‑149 (MQ=255)
                                  cGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTg                                                                                                                <  7:340882/149‑1 (MQ=255)
                                    ttGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGTGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGt                                                                                                              >  4:15402/1‑149 (MQ=255)
                                       gcgcGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCa                                                                                                            <  6:431355/148‑1 (MQ=255)
                                        cgcgAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGTGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAag                                                                                                          <  3:15402/149‑1 (MQ=255)
                                        cgcgAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAag                                                                                                          <  3:130213/149‑1 (MQ=255)
                                         gcgAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCCGGCTGTTGGTCAag                                                                                                          >  7:329731/1‑148 (MQ=255)
                                                gCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTgg                                                                                                  >  3:66386/1‑149 (MQ=255)
                                                    tGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAg                                                                                              <  2:67448/149‑1 (MQ=255)
                                                     gcagcaGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAg                                                                                              <  3:231705/148‑1 (MQ=255)
                                                      cagcagGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAgg                                                                                             <  8:40280/148‑1 (MQ=255)
                                                      cagcagGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGCTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAggt                                                                                            >  1:85105/1‑149 (MQ=255)
                                                       agcagGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAggtg                                                                                           >  2:9366/1‑149 (MQ=255)
                                                         cagGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGACAAGCAGCGTGGCGAggtggt                                                                                         <  1:134707/149‑1 (MQ=255)
                                                              aTTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGAt                                                                                     >  2:417955/1‑148 (MQ=255)
                                                                tGCTGATAGCGCTGCGGGTCGGGCTGGGCGTGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATg                                                                                    >  1:35755/1‑147 (MQ=255)
                                                                    gATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGATTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTgg                                                                              >  4:241728/1‑149 (MQ=255)
                                                                     aTAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGgc                                                                             >  8:221931/1‑149 (MQ=255)
                                                                      tAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGAt                                                                                     >  8:110258/1‑140 (MQ=255)
                                                                      tAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGAt                                                                                     <  7:110258/140‑1 (MQ=255)
                                                                         cgcTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCaa                                                                         <  3:20340/149‑1 (MQ=255)
                                                                                  tCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCgg                                                                <  4:66386/149‑1 (MQ=255)
                                                                                          ggCGAGGTCGTTCGGCAGCCAGTCGTCATGTGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAggg                                                        <  5:219932/149‑1 (MQ=255)
                                                                                                gTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCgg                                                  <  5:380203/149‑1 (MQ=255)
                                                                                                 tCGTTCGGCAGCCAGTCGTCCTGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGt                                                 >  6:161861/1‑149 (MQ=255)
                                                                                                    ttCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCACAAGCGCTGCGGGTGCAGGGGTTCGGCCag                                              >  2:187398/1‑148 (MQ=255)
                                                                                                      cGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATg                                             <  1:9366/148‑1 (MQ=255)
                                                                                                        gCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGtt                                          <  7:221931/149‑1 (MQ=255)
                                                                                                        gCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGtt                                          >  6:3371/1‑149 (MQ=255)
                                                                                                            ccAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAgg                                      <  7:213369/149‑1 (MQ=255)
                                                                                                               gtcgtcATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTTGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGCTTCAGgctc                                   >  3:71238/1‑149 (MQ=255)
                                                                                                               gtcgtcATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGgctc                                   >  3:53697/1‑149 (MQ=255)
                                                                                                               gtcgccATGGCCCCCCAGCTCGGCCAGCAGGCCTTGGCCTGCGTCGAGCCGGGGGCCTTCGTCAGCCTGGCGGTCACGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGgctc                                   <  8:196686/149‑1 (MQ=255)
                                                                                                                  gtcATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGgctctg                                 <  5:118372/148‑1 (MQ=255)
                                                                                                                   tcATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCTAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGgctctgc                                >  6:302100/1‑148 (MQ=255)
                                                                                                                     aTGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGgctctgctc                              >  2:378696/1‑148 (MQ=255)
                                                                                                                        ggCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGGAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAg                          <  8:50772/149‑1 (MQ=255)
                                                                                                                         gCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAgg                         >  1:41997/1‑149 (MQ=255)
                                                                                                                              cAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGttt                      >  4:309350/1‑147 (MQ=255)
                                                                                                                              cAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCa                    <  5:3371/149‑1 (MQ=255)
                                                                                                                                   cGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTaa               >  2:96988/1‑149 (MQ=255)
                                                                                                                                    ggCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAAc              >  3:97156/1‑149 (MQ=255)
                                                                                                                                     gCCAGCAGGCCTTGGCCTGCGTCGCGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAAcc             >  5:118195/1‑149 (MQ=255)
                                                                                                                                     gCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAAcc             <  4:71838/149‑1 (MQ=255)
                                                                                                                                      gacgcagGCCAAGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAAccc            >  5:16648/4‑149 (MQ=255)
                                                                                                                                      gacgcagGCCAAGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAAccc            >  3:73251/4‑149 (MQ=255)
                                                                                                                                         gcagGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAACCCTAg         <  5:302100/149‑1 (MQ=255)
                                                                                                                                         gcagGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAACCCTAg         <  1:417955/149‑1 (MQ=255)
                                                                                                                                              ccTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAACCCTAGCAAAc    <  1:294954/149‑1 (MQ=255)
                                                                                                                                                ttGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAACCCTAGCAAACCt  >  3:416870/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCCCGGGCCCCAGACGAAGCTCACCACCGAGAAACGTTGGCGCGAGTCGCAGTGCAGCAGGTATTGCTGATAGCGCTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGCGCAAGCGCTGCGGGTGCAGGGGTTCGGTCATGGTTCAGGCTCTGCTCTCAGGTTTCAGTTAACCCTAGCAAACCT  >  NC_002947/3650574‑3650866

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