Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,479 A→T 8.3% intergenic (+36/‑40) PP_4824 → / → PP_4825 two‑component system sensor histidine kinase/response regulator/MarC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,489,4790AT 8.3% 153.8 / 4.4 60intergenic (+36/‑40)PP_4824/PP_4825two‑component system sensor histidine kinase/response regulator/MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base A (32/23);  new base T (2/3);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

GCATGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAGGA  >  NC_002947/5489334‑5489625
                                                                                                                                                 |                                                                                                                                                   
gcaTGGACGGACACATGGCCAAACCCGAGGAGTAGTCGCAGTTGCGTGACCTGATTCAATACTGGCCCAATAACCGGGAAGCCAAGCTGGACCCGGCACAAACGTCCTGACCGATATCCAATCCTCTTGAAGGAGCGGCCACGAGacac                                                                                                                                                  >  6:53056/1‑149 (MQ=255)
  aTGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGTGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACaa                                                                                                                                                <  6:363064/149‑1 (MQ=255)
     gacggacACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAgg                                                                                                                                              <  5:240871/148‑1 (MQ=255)
        ggacACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCcgc                                                                                                                                          >  2:309571/1‑149 (MQ=255)
          acacATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGCGCGGCCTTGTGTCGCGAGGCcgctc                                                                                                                                        <  8:54064/149‑3 (MQ=255)
              aTGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGc                                                                                                                                    <  1:360380/149‑1 (MQ=255)
               tGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCa                                                                                                                                   >  2:117711/1‑149 (MQ=255)
                  ccAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATa                                                                                                                                 >  3:319716/1‑148 (MQ=255)
                   cAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATaaa                                                                                                                               >  3:300529/1‑149 (MQ=255)
                        cGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAAccc                                                                                                                            >  1:378996/1‑147 (MQ=255)
                           tGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTTACCGATATCCAATGCTATTGCAGGAGCGGCCTCGCGACACAAGTACGCGCCGGCATAAAACAAAgt                                                                                                                        >  5:111986/1‑147 (MQ=255)
                                 tGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTg                                                                                                                 >  6:200467/1‑149 (MQ=255)
                                   tCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGc                                                                                                                <  5:304406/148‑1 (MQ=255)
                                           gCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATCCTc                                                                                                        >  4:255541/1‑148 (MQ=255)
                                             gTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCAt                                                                                                     >  2:382582/1‑149 (MQ=255)
                                                   tGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAgttg                                                                                               >  5:381256/1‑149 (MQ=255)
                                                      ttCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCCAGGAACACTGCTCATGCTCCATGAgttgtt                                                                                             >  2:297408/1‑148 (MQ=255)
                                                      ttCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAgttgtt                                                                                             <  6:196958/148‑1 (MQ=255)
                                                       tCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTTGTGTCGCGAGGCcgc                                                                                                                                          <  7:174256/102‑1 (MQ=255)
                                                       tCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTTGTGTCGCGAGGCcgc                                                                                                                                          >  8:174256/1‑102 (MQ=255)
                                                              tGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTa                                                                                    >  6:387710/1‑149 (MQ=255)
                                                               gggCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTAc                                                                                   >  4:217387/1‑149 (MQ=255)
                                                                  ccAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCt                                                                                 >  8:205570/1‑148 (MQ=255)
                                                                    aaTCACCGGGAAGCCAAGGTGGACCCGGCACTAACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACCAGGCCGCGCCGGCATAAACCCAAGGAACACAGCTGATGCTCCATGAGTAGTTCAGCGTCTACCTGa                                                                               >  7:18221/1‑148 (MQ=255)
                                                                          cGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCt                                                                        <  3:217387/149‑1 (MQ=255)
                                                                          cGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCt                                                                        >  1:60805/1‑149 (MQ=255)
                                                                           gggAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTc                                                                       <  1:382582/149‑1 (MQ=255)
                                                                                gCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGc                                                                         >  3:301819/1‑142 (MQ=255)
                                                                                gCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGc                                                                         <  4:301819/142‑1 (MQ=255)
                                                                                 ccAAGGTGGACCCGGCACATAAGTCCTGACCGATATACAATGCTCTTGCAGGAGCGGCCTATCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGctc                                                                <  5:53056/149‑1 (MQ=255)
                                                                                   aaGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTa                                                              >  8:41846/1‑149 (MQ=255)
                                                                                         gACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGccc                                                        >  8:87496/1‑149 (MQ=255)
                                                                                          aCCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCg                                                       <  5:200467/149‑1 (MQ=255)
                                                                                              ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTg                                                                   <  4:212176/133‑1 (MQ=255)
                                                                                              ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTg                                                                   >  3:212176/1‑133 (MQ=255)
                                                                                              ggCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGttct                                                   >  1:263517/1‑149 (MQ=255)
                                                                                                cacaTACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGttcttc                                                 <  5:387710/149‑1 (MQ=255)
                                                                                                    tACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGc                                             >  6:220679/1‑149 (MQ=255)
                                                                                                     aCGTCCTGACCGATATCCAATGATCTTGCAGGAGCGGACTCGCGACACAAGGCCGAGCCGGGATAACCCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCATCGTCTACCTGTAAATGCTCGTGCTCTACAGAACGTTCTTCGTGCt                                            >  6:259521/1‑149 (MQ=255)
                                                                                                        tCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTc                                         >  3:20875/1‑149 (MQ=255)
                                                                                                          cTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCgt                                       >  7:157742/1‑149 (MQ=255)
                                                                                                            gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCgtgt                                     >  8:34072/1‑149 (MQ=255)
                                                                                                            gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCgt                                       >  6:308096/1‑147 (MQ=255)
                                                                                                            gACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCgt                                       <  5:308096/147‑1 (MQ=255)
                                                                                                                     ccAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGc                            >  3:333578/1‑149 (MQ=255)
                                                                                                                      cAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGcc                           >  6:66635/1‑149 (MQ=255)
                                                                                                                              ttGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGt                   >  8:332071/1‑149 (MQ=255)
                                                                                                                                gCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTgg                 >  5:19658/1‑149 (MQ=255)
                                                                                                                                 cAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGc                <  7:34072/149‑1 (MQ=255)
                                                                                                                                 cAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGc                <  2:60805/149‑1 (MQ=255)
                                                                                                                                  aGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGcc               >  6:314019/1‑149 (MQ=255)
                                                                                                                                     cgCGGCCTTGTGTCGCGAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCg                                                >  4:324719/2‑113 (MQ=255)
                                                                                                                                     cgCGGCCTTGTGTCGCGAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCg                                                <  3:324719/112‑1 (MQ=255)
                                                                                                                                         gCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAg        <  1:138334/149‑1 (MQ=255)
                                                                                                                                          ccTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTCGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAg        <  1:309571/149‑1 (MQ=255)
                                                                                                                                          ccTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGc       <  5:220679/149‑1 (MQ=255)
                                                                                                                                             cgcgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAg    <  4:333578/149‑1 (MQ=255)
                                                                                                                                              gcgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCaa     <  7:205570/147‑1 (MQ=255)
                                                                                                                                              gcgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAg    <  2:378996/148‑1 (MQ=255)
                                                                                                                                               cgACACAAGGCCGCGCCGGCATAAACCCAAGGAACTCTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAGGa  <  3:255541/149‑1 (MQ=255)
                                                                                                                                               cgACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAGGa  <  4:319716/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCATGGACGGACACATGGCCAAACCGGTGGAGTTGTCGCAGTTGCGTGAGCTGATTCAATACTGGGCCAATCACCGGGAAGCCAAGGTGGACCCGGCACATACGTCCTGACCGATATCCAATGCTCTTGCAGGAGCGGCCTCGCGACACAAGGCCGCGCCGGCATAAACCCAAGGAACACTGCTCATGCTCCATGAGTTGTTCAGCGTCTACCTGAAAATGCTCGTGCTCTACAGCCCGTTCTTCGTGCTGTCGTGTTTCATCAGCCTGACCCGTGGCCACTCCAGCAAGGA  >  NC_002947/5489334‑5489625

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: