Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,979,975 T→A 7.2% L207* (TTA→TAA)  PP_5242 → GAF domain‑containing sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,979,9750TA 7.2% 206.4 / 5.4 69L207* (TTA→TAA) PP_5242GAF domain‑containing sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base A (3/2);  total (36/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 = 8.65e-01

CGTGCCGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAA  >  NC_002947/5979830‑5980119
                                                                                                                                                 |                                                                                                                                                
cGTGCCGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCTTGGGGGCAATCCTCTACCAGTTAc                                                                                                                                                <  7:51576/148‑1 (MQ=255)
 gTGCCGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTg                                                                                                                                              >  8:30053/1‑149 (MQ=255)
   gCCGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTgg                                                                                                                                             >  1:230498/1‑148 (MQ=255)
       aCCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGc                                                                                                                                        <  7:179225/149‑1 (MQ=255)
       aCCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTAACTGGTAgag                                                                                                                                        >  8:8926/1‑147 (MQ=255)
        cccAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGc                                                                                                                                        >  6:59527/1‑148 (MQ=255)
           agcagAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTAACTGGTAgaggatt                                                                                                                                    <  7:8926/149‑7 (MQ=255)
             cagAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTg                                                                                                                                  >  8:260453/1‑149 (MQ=255)
                        aGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAg                                                                                                                                          >  8:179858/1‑130 (MQ=255)
                        aGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAg                                                                                                                                          >  8:175033/1‑130 (MQ=255)
                        aGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAg                                                                                                                                          <  7:175033/130‑1 (MQ=255)
                        aGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAg                                                                                                                                          <  7:179858/130‑1 (MQ=255)
                        aGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCgg                                                                                                                       >  4:177944/1‑149 (MQ=255)
                            aaCAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGgcgc                                                                                                                   <  2:24211/149‑1 (MQ=255)
                               aGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGAGACCGACGAAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGa                                                                                                                 >  2:238879/1‑148 (MQ=255)
                                 ggCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGacac                                                                                                              <  6:340453/149‑1 (MQ=255)
                                     ccccGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTg                                                                                                          <  5:59527/149‑1 (MQ=255)
                                      cccGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGc                                                                                                         >  8:180536/1‑149 (MQ=255)
                                             ggCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGc                                                                                                   <  2:230498/148‑1 (MQ=255)
                                             ggCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGCCGGGCGCGACACCCTGCACTGGc                                                                                                   <  2:166300/147‑1 (MQ=255)
                                                     ggCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTg                                                                                           >  8:173844/1‑148 (MQ=255)
                                                     ggCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTg                                                                                           <  3:177944/148‑1 (MQ=255)
                                                      gCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTgg                                                                                                    >  6:30238/1‑138 (MQ=255)
                                                      gCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTgg                                                                                                    <  5:30238/138‑1 (MQ=255)
                                                      gCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCa                                                                                         <  7:173844/149‑1 (MQ=255)
                                                       ccTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCt                                                                                            <  3:255783/145‑1 (MQ=255)
                                                       ccTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCt                                                                                            >  4:255783/1‑145 (MQ=255)
                                                           ccTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTg                                                                                    <  3:224815/149‑1 (MQ=255)
                                                           ccTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTg                                                                                    <  3:229313/149‑1 (MQ=255)
                                                             tGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGgt                                                                                  >  5:219606/1‑149 (MQ=255)
                                                                  ccGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCa                                                                             <  7:260453/149‑1 (MQ=255)
                                                                    ggAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGc                                                                           >  5:119164/1‑149 (MQ=255)
                                                                     gAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGcc                                                                          >  8:172894/1‑149 (MQ=255)
                                                                      aGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCcg                                                                         >  3:37300/1‑149 (MQ=255)
                                                                           gTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCATTCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCgcgggc                                                                    >  5:82136/1‑149 (MQ=255)
                                                                           gTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCgcgggc                                                                    <  7:172894/149‑1 (MQ=255)
                                                                           gTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCgcggg                                                                     >  6:237566/1‑148 (MQ=255)
                                                                           gTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCgcggg                                                                     <  5:237566/148‑1 (MQ=255)
                                                                           gTTTGCCAGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCgcgggc                                                                    <  8:229709/149‑1 (MQ=255)
                                                                                  cGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGc                                                             >  3:198734/1‑149 (MQ=255)
                                                                                        gTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCg                                                       <  6:119164/149‑1 (MQ=255)
                                                                                           gcagcagcGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGa                                                     <  7:30053/148‑1 (MQ=255)
                                                                                           gcagcagcGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGAc                                                    <  1:238879/149‑1 (MQ=255)
                                                                                             agcagcGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAg                                                  >  8:269638/1‑149 (MQ=255)
                                                                                                 gcGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACGGCAAg                                              <  6:219606/149‑1 (MQ=255)
                                                                                                            aGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTcc                                  >  5:17897/1‑149 (MQ=255)
                                                                                                              tGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAg                                              <  8:343674/136‑1 (MQ=255)
                                                                                                              tGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAg                                              >  7:343674/1‑136 (MQ=255)
                                                                                                              tGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTccc                                 <  7:180536/149‑1 (MQ=255)
                                                                                                               gACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCt                                >  4:197153/1‑149 (MQ=255)
                                                                                                                 cATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCt                                >  6:111364/1‑147 (MQ=255)
                                                                                                                 cATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCt                                >  4:80043/1‑147 (MQ=255)
                                                                                                                  aTCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAg                             >  3:21613/1‑149 (MQ=255)
                                                                                                                   tCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCt                                <  3:152260/145‑1 (MQ=255)
                                                                                                                   tCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCt                                >  4:152260/1‑145 (MQ=255)
                                                                                                                    cTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGcc                           >  7:65513/1‑149 (MQ=255)
                                                                                                                     tATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCt                          >  2:234494/1‑149 (MQ=255)
                                                                                                                       tGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTg                         <  4:198734/148‑1 (MQ=255)
                                                                                                                                ggCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCaa               <  6:82136/149‑1 (MQ=255)
                                                                                                                                    aTCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACTCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGc           <  2:1509/149‑1 (MQ=255)
                                                                                                                                      ccTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCt         >  2:84411/1‑149 (MQ=255)
                                                                                                                                        tctACCAGTTACTGGTAGGCGAACCGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTg       >  5:335553/1‑149 (MQ=255)
                                                                                                                                        cctACCAGTAACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAg             >  6:95227/2‑143 (MQ=255)
                                                                                                                                        cctACCAGTAACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAg             <  5:95227/142‑1 (MQ=255)
                                                                                                                                        cctACCAGTAACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTg       >  3:153710/2‑149 (MQ=255)
                                                                                                                                         ctACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTc                               <  5:3247/124‑1 (MQ=255)
                                                                                                                                         ctACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTc                               >  6:3247/1‑124 (MQ=255)
                                                                                                                                          tACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGc     >  4:11001/1‑149 (MQ=255)
                                                                                                                                             cAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCaa  <  4:37300/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGTGCCGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAA  >  NC_002947/5979830‑5980119

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: