Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,104,368 T→G 5.8% V299G (GTC→GGC)  PP_2714 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,104,3680TG 5.8% 200.2 / 2.1 64V299G (GTC→GGC) PP_2714sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base G (2/2);  total (35/29)
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.92e-01

ACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACATCGACGAAGGGGACCACCACCTGC  >  NC_002947/3104223‑3104511
                                                                                                                                                  |                                                                                                                                               
aCGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTc                                                                                                                                                >  4:214287/1‑147 (MQ=255)
         gCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACt                                                                                                                                     <  5:127620/149‑1 (MQ=255)
         gCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGAGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGGCGAGCCAAGACt                                                                                                                                     >  8:135079/1‑149 (MQ=255)
          ccGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTg                                                                                                                                    <  3:272074/149‑1 (MQ=255)
             tATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGACAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCTCGACGCGGGCCAAggtcagc                                                                                                                                  <  8:130026/148‑7 (MQ=255)
                   gcCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGa                                                                                                                           <  5:146257/149‑1 (MQ=255)
                     cTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGc                                                                                                                         >  3:355415/1‑149 (MQ=255)
                      tGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGc                                                                                                                         >  7:247595/1‑148 (MQ=255)
                        cACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCa                                                                                                                       <  4:355415/148‑1 (MQ=255)
                            cgcAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCAc                                                                                                                   >  1:105224/1‑149 (MQ=255)
                             gcAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCAccc                                                                                                                 >  5:79217/1‑149 (MQ=255)
                             gcAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCGTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCAccc                                                                                                                 <  7:227345/149‑1 (MQ=255)
                                gAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCTCGACGCGGGCCAAggt                                                                                                                                      <  5:317819/125‑3 (MQ=255)
                                gAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCTCGACGCGGGCCAAggt                                                                                                                                      >  6:317819/1‑123 (MQ=255)
                                      gcTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGAcc                                                                                                        <  8:23978/149‑1 (MQ=255)
                                             gccAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCAcc                                                                                                                  <  8:109753/132‑1 (MQ=255)
                                             gccAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCAcc                                                                                                                  >  7:109753/1‑132 (MQ=255)
                                               cAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAgg                                                                                               >  1:203805/1‑149 (MQ=255)
                                                 gCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGgtc                                                                                             <  2:105224/149‑1 (MQ=255)
                                                 gCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGgtc                                                                                             <  6:79217/149‑1 (MQ=255)
                                                       atgCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCAcc                                                                                       >  2:189815/1‑149 (MQ=255)
                                                            gAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAAc                                                                                  >  4:279294/1‑149 (MQ=255)
                                                               cgcCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCt                                                                               >  6:109811/1‑149 (MQ=255)
                                                                  aaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGc                                                                            >  7:238568/2‑149 (MQ=255)
                                                                          ccTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGa                                                                    >  3:82280/1‑149 (MQ=255)
                                                                          ccTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGa                                                                    >  7:285535/1‑149 (MQ=255)
                                                                          ccTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGa                                                                    >  7:82326/1‑149 (MQ=255)
                                                                          ccTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGa                                                                    <  3:279294/149‑1 (MQ=255)
                                                                                  tttCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATg                                                             >  1:137001/1‑149 (MQ=255)
                                                                                  tttCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATg                                                             >  1:31750/1‑149 (MQ=255)
                                                                                    tCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATgg                                                            <  2:137001/148‑1 (MQ=255)
                                                                                    tCTTGTGGGCGGGGGCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATgg                                                            <  8:247595/148‑1 (MQ=255)
                                                                                       tgtgGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATa                                                         <  4:82280/148‑1 (MQ=255)
                                                                                        gtgGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATa                                                         <  8:285535/147‑1 (MQ=255)
                                                                                           ggCGGGTTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTTCCCCCATGGAAACTGGa                                                    >  7:288664/1‑149 (MQ=255)
                                                                                                 gTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGCCCCCATGGATACTGGACCGTgg                                              >  7:277385/1‑149 (MQ=255)
                                                                                                    gATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGcc                                            <  5:109811/148‑1 (MQ=255)
                                                                                                    gATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCt                                           >  4:140905/1‑149 (MQ=255)
                                                                                                     aTCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTg                                          <  3:214287/149‑1 (MQ=255)
                                                                                                      tCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTg                                          >  7:129037/1‑148 (MQ=255)
                                                                                                       cGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCCCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTCACCCCATGTATAATGGACAGTgccaatg                                         >  6:283638/1‑142 (MQ=255)
                                                                                                        ggCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGa                                        <  8:238568/148‑1 (MQ=255)
                                                                                                           tGACTCGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGCGCAGCGAGCCAGCACCCGGATCGACCTGTACTAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAAccc                                    <  5:283638/149‑1 (MQ=255)
                                                                                                           tGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAAccc                                    <  2:203805/149‑1 (MQ=255)
                                                                                                                 gggTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGaaa                                                                                   >  4:213470/1‑96 (MQ=255)
                                                                                                                 gggTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGaaa                                                                                   <  3:213470/96‑1 (MQ=255)
                                                                                                                   gtggtCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCATCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTg                              >  6:199066/1‑147 (MQ=255)
                                                                                                                   gtggtCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCATCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTg                              <  5:199066/147‑1 (MQ=255)
                                                                                                                   gtggtCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGa                            <  1:189815/149‑1 (MQ=255)
                                                                                                                    tggtAAACCAGTTGCTGACCTGGGCCGGCGTCGAGCCAAGACTGGGGCAGCGATCCAGCACCCGGATGGACCTGGCCGAGGTCGTGACAGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTTGAc                           <  8:288664/149‑1 (MQ=255)
                                                                                                                      gtCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACAt                         >  3:284097/1‑149 (MQ=255)
                                                                                                                       tCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGCGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGACGTCACCGAAACCCTGGCTGAACTCACCCCATGGAAACTGGCCCGTCGCCTGGAACCCACGCTGGAAATc                        >  8:192174/1‑149 (MQ=255)
                                                                                                                        cAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGa                                                                    <  6:261860/104‑1 (MQ=255)
                                                                                                                        cAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGa                                                                    >  5:261860/1‑104 (MQ=255)
                                                                                                                         aaCCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACATcga                      >  3:1621/1‑149 (MQ=255)
                                                                                                                           ccAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACATcgac                     >  1:191544/1‑148 (MQ=255)
                                                                                                                                   tGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACATCGACGAAGGGGac            >  8:45423/1‑149 (MQ=255)
                                                                                                                                     cccTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACGCTCGCTGGACATCGACGAAGGGGac            <  7:192174/146‑1 (MQ=255)
                                                                                                                                      ccTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACAGTGGCCTGGAACCCTCGCTGGACATCGACGAAGGGGaccac         <  8:82326/149‑1 (MQ=255)
                                                                                                                                         tGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCATGGCTGAACTGACCCCATGGATACAGGACAGAGGCCTGGAACCCACGCTGGACATCCACCAAGTGGacaaccac      >  7:371098/1‑149 (MQ=255)
                                                                                                                                          ggCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACATCGACGAAGGGGaccaccacc     >  5:164775/1‑149 (MQ=255)
                                                                                                                                             ccGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACATCGACGAAGGGGACCACCACCTg   >  5:268338/1‑148 (MQ=255)
                                                                                                                                             ccGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACATCGACGAAGGGGACCACCACCTGc  <  4:284097/149‑1 (MQ=255)
                                                                                                                                             ccGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCGGGAACCCTCGCTGGACATCGACGAAGGGGACCACCACCTg   >  5:64978/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACCCCATGGATACTGGACCGTGGCCTGGAACCCTCGCTGGACATCGACGAAGGGGACCACCACCTGC  >  NC_002947/3104223‑3104511

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: