Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,119,196 G→C 6.1% E324Q (GAG→CAG)  PP_4505 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,119,1960GC 6.1% 214.2 / 2.5 65E324Q (GAG→CAG) PP_4505sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base C (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.98e-01

AAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCG  >  NC_002947/5119056‑5119337
                                                                                                                                              |                                                                                                                                             
aaGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAg                                                                                                                                       <  5:250841/149‑1 (MQ=255)
    ccTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAg                                                                                                                                   <  8:219077/149‑1 (MQ=255)
       gATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGtt                                                                                                                                >  6:208711/1‑149 (MQ=255)
         tCAGCGAAAGCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGttt                                                                                                                               <  1:112307/148‑1 (MQ=255)
           aGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                             <  3:144587/148‑1 (MQ=255)
           aGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAg                                                                                                                            <  5:179824/149‑1 (MQ=255)
             cGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAT‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCt                                                                                                                          >  4:304891/1‑149 (MQ=255)
             cGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCt                                                                                                                          >  2:256757/1‑149 (MQ=255)
              gAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCtt                                                                                                                         >  6:3388/1‑149 (MQ=255)
                    cTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAgg                                                                                                                    <  5:208711/148‑1 (MQ=255)
                     tCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGt                                                                                                                  >  3:227485/1‑149 (MQ=255)
                       gCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTg                                                                                                                <  1:256757/149‑1 (MQ=255)
                         ccTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGc                                                                                                              >  1:302199/1‑149 (MQ=255)
                          cTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGcc                                                                                                             <  4:227485/149‑1 (MQ=255)
                             gCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAAc                                                                                                          >  4:31320/1‑149 (MQ=255)
                             gCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTCCAGGCGCacgtcc                                                                                                                                   <  3:189734/124‑6 (MQ=255)
                             gCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTCCAGGCGCacgtcc                                                                                                                                   >  4:189734/1‑119 (MQ=255)
                              cccGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACg                                                                                                         <  8:156514/149‑1 (MQ=255)
                              cccGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACg                                                                                                         <  3:31320/149‑1 (MQ=255)
                                        tGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAg                                                                                                <  6:301637/148‑1 (MQ=255)
                                        tGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCCCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAg                                                                                                >  5:301637/1‑148 (MQ=255)
                                         gCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGc                                                                                               <  6:89260/148‑1 (MQ=255)
                                         gCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGc                                                                                               >  5:89260/1‑148 (MQ=255)
                                         gCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGAAGAGGACGTGTGAGGCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAg                                                                                                >  1:265377/1‑149 (MQ=255)
                                           cGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGCC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTg                                                                                            >  6:70608/1‑149 (MQ=255)
                                               cAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGt                                                                                        <  5:3388/149‑1 (MQ=255)
                                                 cgCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTTCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTg                                                                                      <  7:76115/148‑1 (MQ=255)
                                                  gCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGa                                                                                     <  4:314145/149‑1 (MQ=255)
                                                  gCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGa                                                                                     >  7:126270/1‑149 (MQ=255)
                                                          cGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgc                                                                             >  6:87060/1‑149 (MQ=255)
                                                          cGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgc                                                                             >  6:60207/1‑149 (MQ=255)
                                                              ccGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTgcg                                                                         >  7:104170/1‑149 (MQ=255)
                                                                 gTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCaa                                                                      >  7:221663/1‑149 (MQ=255)
                                                                        tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTg                                                               >  4:241390/1‑149 (MQ=255)
                                                                              cccTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCtt                                                          <  5:60207/148‑1 (MQ=255)
                                                                                cTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCaa                                                       >  7:59764/1‑149 (MQ=255)
                                                                                 tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAc                                                      >  7:240283/1‑149 (MQ=255)
                                                                                  aCTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGAc                                                                                    >  4:170759/1‑118 (MQ=255)
                                                                                  aCTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGAc                                                                                    <  3:170759/118‑1 (MQ=255)
                                                                                   cTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAc                                                      <  4:205164/147‑1 (MQ=255)
                                                                                   cTCCGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCCCCGGATCAGGAC‑‑GTGCGCCTGGACGCGCAGCCACGTTTGAGCTTTCACGCTAGGCCAACGCTGATCGACCGTGCCCTAGACAACATCCAGCGCAATGCCCTGCGATACAAccc                                                    >  4:75968/1‑149 (MQ=255)
                                                                                    tCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcg                                                   >  8:288315/1‑149 (MQ=255)
                                                                                      ggCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcggc                                                 <  2:265377/149‑1 (MQ=255)
                                                                                      ggCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcgg                                                  <  8:126270/148‑1 (MQ=255)
                                                                                      ggCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcg                                                   <  8:240283/147‑1 (MQ=255)
                                                                                                gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCg                                       >  8:204428/1‑149 (MQ=255)
                                                                                                  aaGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGt                                     >  3:21172/1‑149 (MQ=255)
                                                                                                  aaGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGt                                     <  8:221663/149‑1 (MQ=255)
                                                                                                            tgctgAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAgc                           >  2:144846/1‑149 (MQ=255)
                                                                                                              ctgAGCGCGCCGGATCAGGAT‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAgcgc                         <  3:304891/149‑1 (MQ=255)
                                                                                                                gAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGcgccg                       >  6:206063/1‑149 (MQ=255)
                                                                                                                     cgcCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCcgcgcgc                  <  4:21172/149‑1 (MQ=255)
                                                                                                                        cGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGc                                             <  8:229005/119‑1 (MQ=255)
                                                                                                                        cGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGc                                             >  7:229005/1‑119 (MQ=255)
                                                                                                                         ggATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCACGCGCGAGc              <  8:104170/149‑1 (MQ=255)
                                                                                                                               ggAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGAcc        <  5:70608/149‑1 (MQ=255)
                                                                                                                               ggAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGAcc        >  4:49112/1‑149 (MQ=255)
                                                                                                                               ggAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGAcc        <  5:206063/149‑1 (MQ=255)
                                                                                                                                   ggctGCGCCTCCTGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGcc                                        <  3:256392/112‑1 (MQ=255)
                                                                                                                                   ggctGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGcc                                        >  4:256392/4‑115 (MQ=255)
                                                                                                                                       gcgcCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAgg           <  7:63177/140‑1 (MQ=255)
                                                                                                                                       gcgcCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAgg           >  8:63177/1‑140 (MQ=255)
                                                                                                                                       gcgcCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCg  <  5:230910/149‑1 (MQ=255)
                                                                                                                                          ccTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCa                                           <  8:150567/105‑1 (MQ=255)
                                                                                                                                          ccTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCa                                           >  7:150567/1‑105 (MQ=255)
                                                                                                                                              |                                                                                                                                             
AAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGAC‑‑GTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCG  >  NC_002947/5119056‑5119337

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