Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,984,442 G→A 6.6% intergenic (+38/+51) PP_5242 → / ← ycaC‑II GAF domain‑containing sensor histidine kinase/putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,984,4420GA 6.6% 170.9 / 2.0 62intergenic (+38/+51)PP_5242/ycaC‑IIGAF domain‑containing sensor histidine kinase/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base A (2/2);  total (34/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 = 4.23e-01

GATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGC  >  NC_002947/5984299‑5984587
                                                                                                                                               |                                                                                                                                                 
gATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCgg                                                                                                                                              <  1:311286/149‑1 (MQ=255)
 aTCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGc                                                                                                                                             >  7:365215/1‑149 (MQ=255)
      tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCATCGCGGGTAAACCCCCTCCGACGGGCGGCTGTAg                                                                                                                                        >  6:239581/1‑149 (MQ=255)
        gTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTACGACCGTCGGCTGTAggg                                                                                                                                      >  8:169878/1‑149 (MQ=255)
            ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                  <  7:88714/149‑1 (MQ=255)
            ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                  <  3:167152/149‑1 (MQ=255)
            ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                  <  4:13147/149‑1 (MQ=38)
                  tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTGCGACGGTCGGCGGTTGGGGCGGGGTTAc                                                                                                                            >  7:133737/1‑149 (MQ=39)
                  tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTg                                                                                                                                           >  5:336252/1‑134 (MQ=255)
                  tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTg                                                                                                                                           <  6:336252/134‑1 (MQ=255)
                  tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTa                                                                                                                                         <  2:44441/136‑1 (MQ=255)
                  tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTa                                                                                                                                         >  1:44441/1‑136 (MQ=255)
                  tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCACCTACGGACGGCTGTAGGGGCGGGTTTAc                                                                                                                            >  5:190770/1‑149 (MQ=255)
                      gCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTAGCGGGTAAAACCGCGCCTTCGGTGGGCGGTATGGGCGGGTTTAACcgc                                                                                                                        >  5:340024/1‑149 (MQ=255)
                       cccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCcgc                                                                                                                        >  3:332753/1‑148 (MQ=255)
                                ggCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGt                                                                                                                                          <  3:352428/121‑1 (MQ=255)
                                ggCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGt                                                                                                                                          >  4:352428/1‑121 (MQ=255)
                                ggCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTGGGGGCGGGTTTACCCGCGAAGAGGCCa                                                                                                              <  6:151231/149‑1 (MQ=255)
                                 gCAGCTGCCGGCGCTGGTGTCGCAGTGGGGGGGATGCCGGATTCGTGTCAGGCTGCCTTTGAACCCCGCATCGCCTGTGCCGGCCTCTTCGCGGGAAAACCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAg                                                                                                             <  5:239581/149‑1 (MQ=37)
                                 gCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGCCGACCGTAGGAGc                                                                                                                                    <  2:339307/126‑1 (MQ=38)
                                 gCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGCCGACCGTAGGAGc                                                                                                                                    >  1:339307/1‑126 (MQ=38)
                                  cAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAGGGTCGGCTGTAGGGGCGGGTTTACCcg                                                                                                                         >  1:133541/1‑136 (MQ=255)
                                  cAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCcg                                                                                                                         <  2:133541/136‑1 (MQ=255)
                                       gCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCTGGAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAgg                                                                                                       >  8:78273/1‑149 (MQ=255)
                                        cAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCGGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAgg                                                                                                       <  1:303650/148‑1 (MQ=255)
                                                  tgtCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAagag                                                                                                                  <  1:50060/127‑1 (MQ=255)
                                                  tgtCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAagag                                                                                                                  >  2:50060/1‑127 (MQ=255)
                                                          gggggTGGATGCCGTATTCTTGCCAGTCTGCCTTCCAACCCCGAATCGCGTGTGCGGCCCTCTGCTCGGTTAACCCGGCTCCTACGGTCGGTTTTAGGGGGGGGTTAACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAgg                                                                                      <  6:340024/145‑1 (MQ=255)
                                                            gggTGCATGCCGGATTCGTGTCAGTCTCCCTTTGAATAACCCATATCCTGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCa                                                                                   <  3:215644/148‑1 (MQ=37)
                                                                   tGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAGCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCg                                                                           >  4:119408/1‑149 (MQ=255)
                                                                   tGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTACTACGGTCGGCTGTAGGGGCGGGGTTAACCGCGAAGAGGCCAGGACAGGTTTACATCa                                                                                              >  1:349669/1‑130 (MQ=255)
                                                                     ccGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCg                                                                           >  1:13185/1‑147 (MQ=255)
                                                                     ccGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGt                                                                          <  4:332753/148‑1 (MQ=255)
                                                                              tgtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTAGCGGGTAAAACCGCTACTACGGTCGGCTGTAGGGGCGGGTTTACCAGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtca                                                                 >  5:123319/1‑148 (MQ=255)
                                                                                tCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcatca                                                              <  2:13185/149‑1 (MQ=255)
                                                                                      tGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTAAATAAGCGCGAGGACAGCCCTGCGGTCCTGGTCCTCCGGTTGc                                                        >  8:70088/1‑149 (MQ=25)
                                                                                        ccTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGc                                                        <  7:169878/147‑1 (MQ=255)
                                                                                        ccTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACCGTCGGGTGGAGGGGCGGGGTTACCCGCGAAGAGGGCAGTACCGGGTTTCAAAAGCCCCAGGCCAGCGCAGCGTAGCTGTTCATCAGGTTGCgg                                                      >  2:84404/1‑149 (MQ=37)
                                                                                               aTCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                               <  8:365215/149‑1 (MQ=255)
                                                                                                tCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTggg                                              >  7:220420/1‑149 (MQ=255)
                                                                                                      atcgCCTGGGCCGGCCTCTTCGCGGGTAACCCCGCCCCTACGGGCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGt                                                                    <  5:248172/118‑1 (MQ=37)
                                                                                                          ccTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCATCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGc                                    <  1:84404/149‑1 (MQ=255)
                                                                                                                 cGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcag                               >  2:316325/1‑147 (MQ=255)
                                                                                                                 cGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagca                             <  8:220420/149‑1 (MQ=255)
                                                                                                                   gCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAgg                           <  7:70088/149‑1 (MQ=255)
                                                                                                                     ctctTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGt                          >  3:354790/1‑148 (MQ=255)
                                                                                                                      tctTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtc                                                                  >  2:328785/1‑107 (MQ=255)
                                                                                                                        ttCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCGGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTAcc                      >  5:206940/1‑149 (MQ=255)
                                                                                                                             gggTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGcc                 >  6:161237/1‑149 (MQ=255)
                                                                                                                             gggTAAACCCGCTCCTACGGCCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGcc                 <  1:45735/149‑1 (MQ=255)
                                                                                                                              ggTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGccc                >  7:204089/1‑149 (MQ=255)
                                                                                                                                tAAACCCGCGCCTACGGGCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCtt              <  3:119408/149‑1 (MQ=255)
                                                                                                                                   ccccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGa           <  8:204089/148‑1 (MQ=255)
                                                                                                                                    cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGGTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt          >  4:205096/1‑149 (MQ=25)
                                                                                                                                     ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg         >  2:270999/1‑149 (MQ=255)
                                                                                                                                     ccGCTCCTACGGTAGGCTGTAGGGGCGGTTTTACCCTCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTTGGGTAGTTGGGTATACGCTGCGACAGCATGTTACCCAGCCCTTCGATg         >  4:118411/1‑149 (MQ=255)
                                                                                                                                      cGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt        >  2:242686/1‑149 (MQ=255)
                                                                                                                                       gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTc       <  6:123319/149‑1 (MQ=255)
                                                                                                                                       gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTc       >  3:172361/1‑149 (MQ=255)
                                                                                                                                         tCCTACGGTCGGCTGTTGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGa                                                                                        <  3:176961/66‑1 (MQ=255)
                                                                                                                                         tCCTACGGTCGGCTCTAGCGGCGGGTTTACCCGCGAAGAGGCCAGTACAGCTTTACATCACCGCGa                                                                                        >  4:176961/1‑66 (MQ=25)
                                                                                                                                            tACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTgc  >  5:371465/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTTGC  >  NC_002947/5984299‑5984587

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