Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,027,746 T→A 7.7% E228D (GAA→GAT PP_3552 ← two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,027,7460TA 7.7% 181.1 / 2.8 61E228D (GAA→GATPP_3552two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base A (3/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGTTCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGT  >  NC_002947/4027600‑4027888
                                                                                                                                                  |                                                                                                                                               
tCACGTTCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTc                                                                                                                                               >  1:283704/1‑149 (MQ=255)
tCACGTTCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCATc                                                                                                                                               >  2:344372/1‑149 (MQ=255)
  aCGTTCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGa                                                                                                                                            >  3:40820/1‑149 (MQ=255)
     ttCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATgc                                                                                                                                         <  4:40820/149‑1 (MQ=255)
             ccTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCATCGAAGCGCaggaaa                                                                                                                                 <  1:344372/149‑6 (MQ=255)
                 cAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATg                                                                                                                              <  3:113228/148‑1 (MQ=255)
                       aCGCGGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCggg                                                                                                                       <  3:31933/149‑1 (MQ=255)
                                 gCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGac                                                                                                             >  8:153680/1‑149 (MQ=255)
                                       gCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTc                                                                                                       >  6:57291/1‑149 (MQ=255)
                                         aTAGTTGACACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTccac                                                                                                    <  5:221048/149‑1 (MQ=255)
                                          tAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTc                                                                                                       <  7:277638/146‑1 (MQ=255)
                                          tAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTc                                                                                                       >  8:277638/1‑146 (MQ=255)
                                            gTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTccacca                                                                                                  >  7:66043/1‑149 (MQ=255)
                                             ttGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCcaccac                                                                                                 >  7:106347/1‑149 (MQ=255)
                                                    ccGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTccacc                                                                                                   <  8:173146/140‑1 (MQ=255)
                                                    ccGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTccacc                                                                                                   >  7:173146/1‑140 (MQ=255)
                                                            aGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGt                                                                                  >  2:265926/1‑149 (MQ=255)
                                                                gCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCTAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCa                                                                              >  5:118694/1‑149 (MQ=255)
                                                                  tgattgatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACg                                                                            >  2:275248/1‑149 (MQ=255)
                                                                    attgatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGa                                                                          <  1:169005/149‑1 (MQ=255)
                                                                    attgatCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGa                                                                          >  2:169005/1‑149 (MQ=255)
                                                                           tcGTGGGCGATGCCCGAGCCCAGCTCGCCCAGCGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATACGCAACGTGATGTCGAGGGACACGCTGACAACCTCCACCACCTCACCGCTGTAGATCTAACTGATCGACCg                                                                   >  8:31869/1‑149 (MQ=255)
                                                                             gTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGcc                                                                     <  1:179391/145‑1 (MQ=255)
                                                                             gTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGc                                                                 >  3:134341/1‑149 (MQ=255)
                                                                             gTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGc                                                                 >  8:315722/1‑149 (MQ=255)
                                                                               gggCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGAGCTTCGATGCGCAAGGTGATGTCGCGGGACACGATTACCACCGCACCCACCGCACCGGGGTAGGTCTCACGGCTAGCCCGGCTg                                                               >  3:295381/1‑149 (MQ=255)
                                                                                ggCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTg                                                               <  1:265926/148‑1 (MQ=255)
                                                                                  cGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGc                                                             >  2:155228/1‑148 (MQ=255)
                                                                                       cccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCAGCCGGGCGAtccg                                                       <  7:153680/149‑4 (MQ=255)
                                                                                        ccGAGGCCAGCTCGCCTAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCAGCCGGGCGAtccgt                                                      <  6:118694/149‑5 (MQ=255)
                                                                                        ccGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCg                                                      <  1:275248/149‑1 (MQ=255)
                                                                                           aGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAAc                                                   <  5:57291/149‑1 (MQ=255)
                                                                                             gCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAAcc                                                  >  6:168712/1‑148 (MQ=255)
                                                                                               cAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTa                                                                                     >  5:45855/1‑111 (MQ=255)
                                                                                               cAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTa                                                                                     <  6:45855/111‑1 (MQ=255)
                                                                                                  cTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGt                                            <  5:168712/149‑1 (MQ=255)
                                                                                                       ccAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACTGATCGCCCGCAATGAGATATCTAACCAACGTTAGcct                                       >  6:50494/1‑149 (MQ=255)
                                                                                                         aGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctcc                                     >  1:12355/1‑149 (MQ=255)
                                                                                                         aGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctc                                      <  4:72963/148‑1 (MQ=255)
                                                                                                         aGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctc                                      >  3:72963/1‑148 (MQ=255)
                                                                                                         aGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctc                                      <  2:59195/148‑1 (MQ=255)
                                                                                                         aGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCACCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGcctc                                      >  1:59195/1‑148 (MQ=255)
                                                                                                               gCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCg                                                                   <  4:224244/113‑1 (MQ=255)
                                                                                                               gCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCg                                                                   >  3:224244/1‑113 (MQ=255)
                                                                                                                cGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTccacc                                                                                                   <  4:245905/80‑1 (MQ=255)
                                                                                                                cGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTccacc                                                                                                   >  3:245905/1‑80 (MQ=255)
                                                                                                                 gAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCt                             <  8:66043/149‑1 (MQ=255)
                                                                                                                  aGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGtct                                                                                >  7:305060/1‑97 (MQ=255)
                                                                                                                  aGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGtct                                                                                <  8:305060/97‑1 (MQ=255)
                                                                                                                      gggCACTGTGCGCCAGGCTTTCCTGCGCTTCGGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTg                            >  1:219281/1‑146 (MQ=255)
                                                                                                                      gggCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTg                            <  2:219281/145‑1 (MQ=255)
                                                                                                                          agtgtgCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTAccc                    <  7:31869/147‑1 (MQ=255)
                                                                                                                            tgtgCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCgg                  >  5:279610/1‑149 (MQ=255)
                                                                                                                                cgcCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCCCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAgg              >  4:26028/1‑149 (MQ=255)
                                                                                                                                      gCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAgg              <  4:122411/143‑1 (MQ=255)
                                                                                                                                      gCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAgg              >  3:122411/1‑143 (MQ=255)
                                                                                                                                          ccttGCGCATCGAAGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGtgtg    >  2:54222/4‑149 (MQ=255)
                                                                                                                                           cttGCGCATCGAAGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCACACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGtgtg    <  2:309495/146‑1 (MQ=255)
                                                                                                                                           cttGCGCATCGAAGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCACACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGtgtg    >  1:309495/3‑148 (MQ=255)
                                                                                                                                           ccTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGtgtg    <  2:283704/148‑1 (MQ=255)
                                                                                                                                            cTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGt  <  4:134341/149‑1 (MQ=255)
                                                                                                                                                  ttCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc         <  3:338288/136‑1 (MQ=255)
                                                                                                                                                  ttCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGATCTCACTGATCGCCCTAAAATAAAAAAAAAAACAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATc         >  4:338288/1‑136 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCACGTTCCAGGCCCTGCAGGTAACGCTGACTGGCGTTGGCATAGTTGACCACCGCGGCCAGTGGCTGATTGATCTCGTGGGCGATGCCCGAGGCCAGCTCGCCCAGGGTAGCGAGGCGGGCACTGTGCGCCAGGCTTTCCTGCGCTTCGATGCGCAAGGTGATGTCGCGGGACACGCTGACCACCTCCACCACCGCACCGGTGTAGGTCTCACGGATCGCCCGGCTGGCGATTTCGAACCAACGGTAGCCTCCTGTCGCCTGGCGTACCCGGCAGGTCATCGTGTGGT  >  NC_002947/4027600‑4027888

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: