Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,884,638 C→T 6.8% L175L (CTG→TTG)  PP_3429 → two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,884,6380CT 6.8% 172.6 / 2.8 59L175L (CTG→TTG) PP_3429two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base T (2/2);  total (29/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 = 9.61e-01

CGCAACCGGCGGCGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGTTTCACGCGCCG  >  NC_002947/3884491‑3884785
                                                                                                                                                   |                                                                                                                                                   
cgcAACCGGCGGCGACAGTACGTGACACGTCATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTCACAGGCACGCAAGATGGAGGCCATCGGCCCACt                                                                                                                                                    >  6:269105/1‑149 (MQ=255)
    aCCGGCGGCGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGcc                                                                                                                                                >  6:205398/1‑149 (MQ=255)
       ggcggcGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgat                                                                                                                                             <  5:316531/149‑3 (MQ=255)
        gcggcgACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCcggcg                                                                                                                                            <  8:41853/149‑1 (MQ=255)
             gACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTc                                                                                                                                       <  5:205398/149‑1 (MQ=255)
             gACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTc                                                                                                                                       >  4:18583/1‑149 (MQ=255)
              aCAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCg                                                                                                                                      <  4:228458/149‑1 (MQ=255)
                aGTACGTGACACGTGATCAACTCCATTGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGcc                                                                                                                                    >  8:37529/1‑149 (MQ=255)
                    cGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCGGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCAt                                                                                                                                 <  7:316821/148‑1 (MQ=255)
                    cGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCAt                                                                                                                                 <  6:325/148‑1 (MQ=255)
                     gTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGa                                                                                                                               <  3:18583/149‑1 (MQ=255)
                      tGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGAc                                                                                                                              >  2:266486/1‑149 (MQ=255)
                        acacGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACtt                                                                                                                            <  5:248190/149‑1 (MQ=255)
                         cacGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTc                                                                                                                           <  4:145215/149‑1 (MQ=255)
                          acGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgatg                                                                                                                                            <  4:20590/131‑4 (MQ=255)
                          acGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgatg                                                                                                                                            >  3:20590/1‑128 (MQ=255)
                          acGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTc                                                                                                                           >  2:322954/1‑148 (MQ=255)
                                  aaCTCCATGGCATGCAGCAGCGGCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTg                                                                                                                  <  5:269105/149‑1 (MQ=255)
                                       cATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCAc                                                                                                             <  1:266486/149‑1 (MQ=255)
                                           gCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAg                                                                                                          >  5:12423/1‑148 (MQ=255)
                                           gCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAg                                                                                                          >  6:246163/1‑148 (MQ=255)
                                            caagcagcaGCGCCTGGCCAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCa                                                                                                        <  8:283035/146‑1 (MQ=255)
                                                cagcagCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAg                                                                                                    >  4:74270/1‑149 (MQ=255)
                                                    agcGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCgg                                                                                                 >  3:262797/1‑148 (MQ=255)
                                                        ccTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCt                                                                                            >  2:255900/1‑149 (MQ=255)
                                                             aCAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGaa                                                                                       <  7:100179/149‑1 (MQ=255)
                                                               aGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACGTGCTCACCAGCATCAGCTGCCGCtgatgctg                                                                                     <  5:246163/149‑8 (MQ=255)
                                                               aGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACt                                                                                     >  2:249408/1‑149 (MQ=255)
                                                                    gCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATcc                                                                                <  2:49967/149‑1 (MQ=255)
                                                                      cTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAg                                                                              <  1:322954/149‑1 (MQ=255)
                                                                                 gaagaaGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCACCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTgcg                                                                   >  7:14837/1‑149 (MQ=255)
                                                                                     aagTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCaa                                                               <  3:74270/149‑1 (MQ=255)
                                                                                          cAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCg                                                          <  1:255900/149‑1 (MQ=255)
                                                                                            ggACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCg                                                          >  1:174676/1‑147 (MQ=255)
                                                                                             gACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAg                                                       >  6:117540/1‑149 (MQ=255)
                                                                                                  cgGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACg                                                  >  7:317351/1‑149 (MQ=255)
                                                                                                           cGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGAt                                          <  2:174676/148‑1 (MQ=255)
                                                                                                            gCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATgg                                        >  7:215393/1‑149 (MQ=255)
                                                                                                             cccGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGc                                       <  8:317351/149‑1 (MQ=255)
                                                                                                                  tGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCt                                  <  4:262797/149‑1 (MQ=255)
                                                                                                                      aGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGcc                                                                 <  2:257729/114‑1 (MQ=255)
                                                                                                                        gCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCAt                            >  8:28731/1‑149 (MQ=255)
                                                                                                                            gCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCTGCGCCTGCGCCAAGGCCGCAGCCACGGCCTGGATGGCGTTCTGTGCATGGGc                        >  2:86397/1‑149 (MQ=255)
                                                                                                                            gCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGc                        >  2:146720/1‑149 (MQ=255)
                                                                                                                            gCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGc                        >  3:172986/1‑149 (MQ=255)
                                                                                                                            gCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGc                        <  8:215393/149‑1 (MQ=255)
                                                                                                                               aGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCcgc                     >  2:21850/1‑149 (MQ=255)
                                                                                                                                gATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCcgc                     >  1:215785/1‑148 (MQ=255)
                                                                                                                                  tGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAg                  <  7:85510/149‑1 (MQ=255)
                                                                                                                                    gAGGCCATCGGCCAACTGGCCGGCGGCGTCACCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGc                <  8:14837/149‑1 (MQ=255)
                                                                                                                                     aGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCTGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGcc               >  2:81699/1‑149 (MQ=255)
                                                                                                                                      ggCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCt                                  >  1:218971/1‑129 (MQ=255)
                                                                                                                                      ggCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCt                                  <  2:218971/129‑1 (MQ=255)
                                                                                                                                        gccacGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGtt            <  5:175776/145‑1 (MQ=255)
                                                                                                                                         cgccGGCCAGTTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGttt           >  3:71238/4‑149 (MQ=255)
                                                                                                                                          aTCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGTTTc          <  2:215785/149‑1 (MQ=255)
                                                                                                                                           tCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGTTTCa         <  5:117540/149‑1 (MQ=255)
                                                                                                                                                 aaCTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGTTTCACGcgcc   >  7:105422/1‑149 (MQ=255)
                                                                                                                                                  aCTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGTTTCACGcgccg  >  3:27607/1‑149 (MQ=255)
                                                                                                                                                   cTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGc                >  4:270461/1‑134 (MQ=255)
                                                                                                                                                   cTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGc                <  3:270461/134‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCAACCGGCGGCGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGTTTCACGCGCCG  >  NC_002947/3884491‑3884785

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: ‑