Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,979,975 T→A 8.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 8.2% 174.8 / 5.5 61L207* (TTA→TAA) PP_5242GAF domain‑containing sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base A (2/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GTGCCGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGG  >  NC_002947/5979831‑5980115
                                                                                                                                                |                                                                                                                                            
gTGCCGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTg                                                                                                                                          <  5:357171/149‑1 (MQ=255)
   ccGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTAACTgg                                                                                                                                         >  2:237678/1‑147 (MQ=255)
      aCCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTAACTGGTAgag                                                                                                                                    <  7:294703/149‑3 (MQ=255)
      aCCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTAACTGGTAgag                                                                                                                                    <  8:130890/149‑3 (MQ=255)
       cccAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTAACTGGTAgagg                                                                                                                                   <  1:237678/149‑4 (MQ=255)
         cagcagAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGaa                                                                                                                                 >  7:348631/1‑149 (MQ=255)
         cagcagAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGaa                                                                                                                                 <  4:63222/149‑1 (MQ=255)
          agcagAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGGCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAAc                                                                                                                                <  3:207341/149‑1 (MQ=255)
              gagCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGcc                                                                                                                            >  5:43973/1‑149 (MQ=255)
              gagCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGcc                                                                                                                            <  2:1074/149‑1 (MQ=255)
              gagCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGcc                                                                                                                            >  3:66245/1‑149 (MQ=255)
                     tCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGcc                                                                                                                     >  7:280045/1‑149 (MQ=255)
                      cAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGc                                                                                                                      <  1:342565/147‑1 (MQ=255)
                      cAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGc                                                                                                                      >  2:342565/1‑147 (MQ=255)
                        gCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCggg                                                                                                                  >  7:160446/1‑149 (MQ=255)
                         cGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGgc                                                                                                                 <  6:43973/149‑1 (MQ=255)
                          gAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAg                                                                                                                                      <  4:117130/127‑1 (MQ=255)
                          gAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAg                                                                                                                                      >  3:117130/1‑127 (MQ=255)
                           aaCAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGCCCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGgcgc                                                                                                               >  3:316157/1‑149 (MQ=255)
                           aaCAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGgcgc                                                                                                               <  8:348631/149‑1 (MQ=255)
                           aaCAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGgcgc                                                                                                               >  6:42134/1‑149 (MQ=255)
                              aGGGCAGCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGa                                                                                                             >  1:63161/1‑148 (MQ=255)
                                ggCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGaca                                                                                                           >  8:48293/1‑148 (MQ=255)
                                ggCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGaca                                                                                                           <  7:48293/148‑1 (MQ=255)
                                     cccGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTg                                                                                                      >  4:37744/1‑148 (MQ=255)
                                        gCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACt                                                                                                  >  5:64631/1‑149 (MQ=255)
                                          tCGGCAACTGGGCCGACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTgg                                                                                                >  7:314664/1‑149 (MQ=255)
                                            ggCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGc                                                                                               >  1:232312/1‑148 (MQ=255)
                                                      ccTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCAt                                                                                    <  3:37744/149‑1 (MQ=255)
                                                      ccTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCAt                                                                                    <  6:64631/149‑1 (MQ=255)
                                                      ccTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCAt                                                                                    >  6:326909/1‑149 (MQ=255)
                                                         aCCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGc                                                                                      >  1:170756/1‑144 (MQ=255)
                                                         aCCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGc                                                                                      <  2:170756/144‑1 (MQ=255)
                                                         aCCTGGAGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCt                                                                                 >  7:94087/1‑149 (MQ=255)
                                                           cTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGTCATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGATCTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTgg                                                                               <  8:160446/149‑1 (MQ=255)
                                                            tGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGgt                                                                              >  2:190092/1‑149 (MQ=255)
                                                             ggCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTTCCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGgtg                                                                             >  4:321133/1‑149 (MQ=255)
                                                              gcgcCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTTCCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGgtg                                                                             <  3:321133/148‑1 (MQ=255)
                                                              gcgcCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGgt                                                                              <  5:42134/147‑1 (MQ=255)
                                                                         ggTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCgcggg                                                                 >  5:306614/1‑149 (MQ=255)
                                                                          gTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATg                                                                                   <  1:349961/130‑1 (MQ=255)
                                                                          gTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATg                                                                                   >  2:349961/1‑130 (MQ=255)
                                                                              gCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCGTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGggcggc                                                             <  7:27364/148‑1 (MQ=255)
                                                                               cccGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCa                                                           <  8:314664/149‑1 (MQ=255)
                                                                                ccGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAg                                                          >  8:272003/1‑149 (MQ=255)
                                                                                       gTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCg                                                   >  5:280336/1‑149 (MQ=255)
                                                                                        tAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCg                                                   >  6:169851/1‑148 (MQ=255)
                                                                                           cagcagcGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACAcc                                                                                                         >  3:138729/1‑91 (MQ=255)
                                                                                           cagcagcGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACAcc                                                                                                         <  4:138729/91‑1 (MQ=255)
                                                                                            agcagcGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAg                                              <  7:272003/149‑1 (MQ=255)
                                                                                                  gaccgacGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAg                                              >  1:25491/1‑143 (MQ=255)
                                                                                                  gaccgacGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAg                                              <  2:25491/143‑1 (MQ=255)
                                                                                                     cgacGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGcc                                     <  1:190092/149‑1 (MQ=255)
                                                                                                                   cTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGcc                       >  3:280726/1‑149 (MQ=255)
                                                                                                                      tGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTg                     <  4:280726/148‑1 (MQ=255)
                                                                                                                               ggCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTAGATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCa            >  1:300160/1‑148 (MQ=255)
                                                                                                                                  aaTCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCgcgggc                                                                <  4:309712/93‑1 (MQ=255)
                                                                                                                                  aaTCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCgcgggc                                                                >  3:309712/1‑93 (MQ=255)
                                                                                                                                       cctACCAGTAACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTg   >  3:55892/2‑149 (MQ=255)
                                                                                                                                           ccAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTgg  <  5:326909/146‑1 (MQ=255)
                                                                                                                                           ccAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTgg  <  4:316157/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GTGCCGACCCAGCAGAGCGGATCAGCGAACAGGGCACCCCGCTCGGCAACTGGGCCTACCTGGCGCCGGAGCAGGTTTGCCCGCACGGTAGCAGCAGCGACCGACGCAGTGACATCTATGCCCTGGGGGCAATCCTCTACCAGTTACTGGTAGGCGAACTGCCGCTGGCCGGGCGCGACACCCTGCACTGGCGACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGG  >  NC_002947/5979831‑5980115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑