Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,711,864 T→C 10.9% 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,8640TC10.9% 165.4 / 10.3 64F235L (TTC→CTC) ubiEubiquinone/menaquinone biosynthesis C‑methyltransferase UbiE
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base C (2/5);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

GACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCCGCG  >  NC_002947/5711726‑5712011
                                                                                                                                          |                                                                                                                                                   
gacgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatca                                                                                                                                           <  3:130525/149‑5 (MQ=255)
 acgCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatcat                                                                                                                                          <  8:169306/149‑6 (MQ=255)
   gCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatcatgc                                                                                                                                        <  6:240808/149‑8 (MQ=255)
       aCTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTAc                                                                                                                                    >  3:173538/1‑149 (MQ=255)
         tCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACca                                                                                                                                  >  6:299564/1‑149 (MQ=255)
           gTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACcac                                                                                                                                 >  7:263151/1‑148 (MQ=255)
           gTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACcac                                                                                                                                 >  2:308077/1‑148 (MQ=255)
           gTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatcatg                                                                                                                                         >  1:97957/1‑134 (MQ=255)
           gTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatcatg                                                                                                                                         <  2:102005/140‑7 (MQ=255)
           gTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatcatg                                                                                                                                         <  2:97957/140‑7 (MQ=255)
           gTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAAGCCGGCCTCGAccatcatg                                                                                                                                         >  1:102005/1‑134 (MQ=255)
            ttCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCacaa                                                                                                                               <  3:320113/149‑1 (MQ=255)
             tCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCt                                                                                                                                      >  6:88448/1‑141 (MQ=255)
             tCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCt                                                                                                                                      <  5:88448/141‑1 (MQ=255)
               gCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCacaaca                                                                                                                             <  2:144581/148‑1 (MQ=255)
               gCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACAt                                                                                                                            >  4:124906/1‑149 (MQ=255)
                 cTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCacaac                                                                                                                              >  7:92331/1‑145 (MQ=255)
                 cTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCacaac                                                                                                                              <  8:92331/145‑1 (MQ=255)
                      tGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGc                                                                                                                     >  4:19560/1‑149 (MQ=255)
                       gCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGa                                                                                                                                                 >  4:233348/1‑120 (MQ=255)
                       gCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGa                                                                                                                                                 <  3:233348/120‑1 (MQ=255)
                        ccGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCg                                                                                                                    <  4:173538/148‑1 (MQ=255)
                               ctggcAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGtt                                                                                                            <  5:299564/149‑1 (MQ=255)
                                                aaTGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGAt                                                                                           <  7:56412/149‑1 (MQ=255)
                                                  tGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATc                                                                                          >  6:162234/1‑148 (MQ=255)
                                                      gCCGAGAGTTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCg                                                                                                        <  6:224989/130‑1 (MQ=255)
                                                      gCCGAGAGATACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCg                                                                                                        >  5:224989/1‑130 (MQ=255)
                                                           gagCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGAt                                                                                <  5:162234/149‑1 (MQ=255)
                                                             gCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGt                                                                              <  4:263502/149‑1 (MQ=255)
                                                                aCCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCAGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGc                                                                           >  1:92242/1‑149 (MQ=255)
                                                                 ccGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGc                                                                           <  8:263151/148‑1 (MQ=255)
                                                                 ccGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCt                                                                          >  2:94446/1‑149 (MQ=255)
                                                                  cGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTg                                                                         <  1:94446/149‑1 (MQ=255)
                                                                         tCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGgc                                                                  >  1:280093/1‑149 (MQ=255)
                                                                           gCCGAGTCGATCCACATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGa                                                                                 <  3:180759/132‑1 (MQ=255)
                                                                           gCCGAGTCGATCCACATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGa                                                                                 >  4:180759/1‑132 (MQ=255)
                                                                              gAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTc                                                             >  3:16090/1‑149 (MQ=255)
                                                                                       cGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTc                                                    <  5:51256/149‑1 (MQ=255)
                                                                                       cGCATGCACCCGGACCAGGAAACCCTGAAGACCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGt                                                     <  3:19560/148‑1 (MQ=255)
                                                                                        gCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTc                                                    >  6:178874/1‑148 (MQ=255)
                                                                                         cATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCg                                                   >  8:137479/1‑148 (MQ=255)
                                                                                         cATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGa                                                  <  5:178874/149‑1 (MQ=255)
                                                                                            gCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACa                                               >  5:263189/1‑149 (MQ=255)
                                                                                             ttcccGGACCAGTAAACCCTGAATTCCATTTTGGGCGTGGCCGGCTTCCCCGGGGGATCCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAAc                                                <  2:258683/145‑1 (MQ=255)
                                                                                             cACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACAt                                              <  7:137479/149‑1 (MQ=255)
                                                                                              aCCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATg                                             >  6:172263/1‑149 (MQ=255)
                                                                                               cccGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATgg                                            <  4:16090/149‑1 (MQ=255)
                                                                                                ccGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAAc                                                >  1:258683/1‑144 (MQ=255)
                                                                                                ccGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGc                                           >  4:210683/1‑149 (MQ=255)
                                                                                                ccGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGc                                           >  1:55511/1‑149 (MQ=255)
                                                                                                 cGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGcc                                          <  3:124906/149‑1 (MQ=255)
                                                                                                       aGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAAcc                                    <  1:308077/149‑1 (MQ=255)
                                                                                                           aaCCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACcgcg                                 >  1:319703/1‑148 (MQ=255)
                                                                                                              ccTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCt                             >  1:165092/1‑149 (MQ=255)
                                                                                                               cTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCGAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTg                            >  5:92368/1‑149 (MQ=255)
                                                                                                                 gAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAAcc                                    >  4:172182/1‑139 (MQ=255)
                                                                                                                 gAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAAcc                                    <  3:172182/139‑1 (MQ=255)
                                                                                                                      gCATGATGGTCGAGGCCGGCTTCTACCGCGTCACCTTCCACAAGATGACCAGCGGGATCGTTGCCGGGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTgg                     <  6:263189/149‑1 (MQ=255)
                                                                                                                      gCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTg                      <  2:319703/148‑1 (MQ=255)
                                                                                                                       cttgatgGGCGTGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGa                    <  2:165092/147‑1 (MQ=255)
                                                                                                                              gTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCGAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACg             <  6:92368/149‑1 (MQ=255)
                                                                                                                                      cagCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGcc     <  2:92242/147‑1 (MQ=255)
                                                                                                                                        gCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCcgc   <  2:55511/149‑1 (MQ=255)
                                                                                                                                         cTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCcgcg  >  7:306510/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GACGCCTACTCGTTCGCCTTCATGCCGCTGGCTGGCAAGCTGATCACCAATGACGCCGAGAGCTACCGTTACCTCGCCGAGTCGATCCGCATGCACCCGGACCAGGAAACCCTGAAGAGCATGATGGTCGAGGCCGGCTTCGACCGCGTCACCTACCACAACATGACCAGCGGCATCGTTGCCGTGCACCGGGGGATCAAGCCCTGATGTTGCTGGCCGGGCTGCTCGCCAGCGTCGAACATGGCCTGAACCGCGTCCTGCGCCTGGACAGCACGGCCCTGCCGCG  >  NC_002947/5711726‑5712011

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: