Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,711,864 T→C 6.3% F235L (TTC→CTC)  ubiE → ubiquinone/menaquinone biosynthesis C‑methyltransferase UbiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,711,8640TC 6.3% 177.8 / 2.7 63F235L (TTC→CTC) ubiEubiquinone/menaquinone biosynthesis C‑methyltransferase UbiE
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base C (2/2);  total (33/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 = 4.98e-01

AAGGCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCCGCGG  >  NC_002947/5711717‑5712012
                                                                                                                                                   |                                                                                                                                                     
aaGGCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCtt                                                                                                                                                      <  6:137702/149‑1 (MQ=255)
   gCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGCCCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGa                                                                                                                                                   >  5:85799/1‑149 (MQ=255)
    ccTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGAc                                                                                                                                                  <  4:79464/149‑1 (MQ=255)
       acgacgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACcgc                                                                                                                                               >  1:129874/1‑149 (MQ=255)
        cgacgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACcgc                                                                                                                                               <  6:122833/148‑1 (MQ=255)
         gacgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatca                                                                                                                                             <  4:25346/149‑5 (MQ=255)
          acgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTc                                                                                                                                            <  4:277462/149‑1 (MQ=255)
          acgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatcat                                                                                                                                            <  4:194811/149‑6 (MQ=255)
           cgCCTACTCGTTCGCCTTCATGCCGGTGGCTGGCAAGCTGATCACCTATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAACCCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTct                                                                                                                                           <  6:309057/149‑2 (MQ=255)
            gCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCAc                                                                                                                                          <  3:128138/149‑1 (MQ=255)
              cTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTc                                                                                                                                                     >  8:270844/1‑136 (MQ=255)
              cTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTc                                                                                                                                                     <  7:270844/136‑1 (MQ=255)
                  tCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACca                                                                                                                                    <  3:125160/149‑1 (MQ=255)
                    gTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAACCCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGCCCGCGTCACCTAcccc                                                                                                                                   >  6:55002/1‑146 (MQ=255)
                        gcactCATGCCGCTGGCTGGCAAGCTGATCACCACTCACGCCGAGACCTACCGTTACCTCCCCGAGTCGATCCGCATGCACCCGGACCAGGAACCCCTGATGAGCAGGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACAt                                                                                                                              <  6:215201/145‑1 (MQ=255)
                               tGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGc                                                                                                                      >  7:67748/1‑149 (MQ=255)
                                 ccGCTGGATGGCAAGCTGATCACCAATGACGCAGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATACACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCg                                                                                                                     <  1:49788/148‑1 (MQ=255)
                                         tggcAAGCTGATCACCAATGACGCAGAGAGCTACCGTTACCTCGCCGAGTCGCTCCGCATGCACCCGGACCAGGAACCCCTGAAGAGCATCATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGAGGCATCCTTg                                                                                                            >  6:22096/1‑149 (MQ=255)
                                            cAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCg                                                                                                         >  6:316378/1‑149 (MQ=255)
                                              aGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTg                                                                                                       <  2:56696/149‑1 (MQ=255)
                                               gCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGc                                                                                                      >  7:133993/1‑149 (MQ=255)
                                                      aCCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCgggg                                                                                               >  5:78101/1‑149 (MQ=255)
                                                      aCCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCggg                                                                                                <  6:297909/148‑1 (MQ=255)
                                                       ccAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCTGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCgggg                                                                                               >  2:287939/1‑148 (MQ=255)
                                                       ccAAAGACTCCAATAGCTACCCATACCTCGCCGACTCCAACCGCAAGCACCCGGACAAGGAAACCCTCAATACCAACATGCTCAAGGCCAGCTTCCACAGCGTCACCTACCACAACATGACCAGCGGCATAGATTCCGAACACCgcgc                                                                                               >  4:3092/1‑145 (MQ=255)
                                                              cGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTTAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGc                                                                                       <  8:211053/149‑1 (MQ=255)
                                                               gCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGcc                                                                                      >  2:93683/1‑149 (MQ=255)
                                                                 cGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCt                                                                                    >  8:34536/1‑149 (MQ=255)
                                                                          ccGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCt                                                                           <  6:78101/149‑1 (MQ=255)
                                                                                  tCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGTAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCgg                                                                                                 <  1:104405/119‑1 (MQ=255)
                                                                                  tCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCgg                                                                                                 >  2:104405/1‑119 (MQ=255)
                                                                                     ccGAGTCGATCCGCATGCACCCGGACCAGTAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACAACAACATCTCCCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGgctg                                                                 <  1:93683/149‑1 (MQ=255)
                                                                                            gATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAg                                                         >  1:191389/1‑149 (MQ=255)
                                                                                              tCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCg                                                       >  2:174449/1‑149 (MQ=255)
                                                                                                   aTGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGaa                                                  <  7:34536/149‑1 (MQ=255)
                                                                                                    tGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGaa                                                  >  8:263054/1‑148 (MQ=255)
                                                                                                     gCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACa                                                >  7:260320/1‑149 (MQ=255)
                                                                                                       aCCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATg                                              >  5:317246/1‑149 (MQ=255)
                                                                                                       aCCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATg                                              >  2:158285/1‑149 (MQ=255)
                                                                                                       aCCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATg                                              >  2:179121/1‑149 (MQ=255)
                                                                                                                     aCCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTc                                <  1:174449/149‑1 (MQ=255)
                                                                                                                       ccTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCt                              >  2:268090/1‑149 (MQ=255)
                                                                                                                           aaGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTgc                            <  6:85799/147‑1 (MQ=255)
                                                                                                                              accATGAGGGGCGAGGCCGGCTGCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTg                       <  8:99510/147‑1 (MQ=255)
                                                                                                                                cttgatgGGCGAGTGCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGAGCAAGCCCTGATGTGGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGa                     <  5:22096/147‑1 (MQ=255)
                                                                                                                                cATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGTGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGa                     >  4:261075/1‑149 (MQ=255)
                                                                                                                                cATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGa                     >  7:203984/1‑149 (MQ=255)
                                                                                                                                     tgGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTgc                            >  5:287709/1‑137 (MQ=255)
                                                                                                                                     tgGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTgc                            <  6:287709/137‑1 (MQ=255)
                                                                                                                                     cggTCGAAGCCGGCCTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCa                >  7:279601/2‑149 (MQ=255)
                                                                                                                                      ggTCGAAGCCGGCCTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCAc               >  6:90261/1‑149 (MQ=255)
                                                                                                                                        tCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACgg             >  6:47818/1‑149 (MQ=255)
                                                                                                                                        tCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACgg             >  4:38846/1‑149 (MQ=255)
                                                                                                                                        tCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACgg             <  1:268090/149‑1 (MQ=255)
                                                                                                                                        tCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACgg             <  1:179121/149‑1 (MQ=255)
                                                                                                                                         cGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTGGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGc            <  5:55002/149‑1 (MQ=255)
                                                                                                                                             gCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTg        >  3:70765/1‑149 (MQ=255)
                                                                                                                                               cGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGc                                                        <  4:141003/99‑1 (MQ=255)
                                                                                                                                               cGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGc                                                        >  3:141003/1‑99 (MQ=255)
                                                                                                                                               cGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTc                                >  3:285569/1‑123 (MQ=255)
                                                                                                                                                ggCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCcg     <  5:47818/149‑1 (MQ=255)
                                                                                                                                                  cTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCcgcg   <  6:317246/149‑1 (MQ=255)
                                                                                                                                                  cTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCcgc    >  6:112623/1‑148 (MQ=255)
                                                                                                                                                  cTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCcgc    <  5:112623/148‑1 (MQ=255)
                                                                                                                                                   ttCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCCGCgg  >  8:175979/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGGCCTACGACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCCGCGG  >  NC_002947/5711717‑5712012

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: