Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,980,171 T→C 9.0% A272A (GCT→GCC 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,980,1710TC 9.0% 184.0 / 6.6 66A272A (GCT→GCCPP_5242GAF domain‑containing sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base C (2/4);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAGGCCTT  >  NC_002947/5980023‑5980311
                                                                                                                                                    |                                                                                                                                             
gACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCt                                                                                                                                               >  2:120984/1‑149 (MQ=255)
     cTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTg                                                                                                                                          >  1:305305/1‑149 (MQ=255)
                   gCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCCCACTGCCGCTGGCAGTgag                                                                                                                            <  4:170445/149‑3 (MQ=255)
                     aGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCCCGGCCGCGGGCAGGGGgcc                                                                                                                          >  1:202783/1‑149 (MQ=255)
                       ccGCGGGCGGCCAGCGCAGCGGACAGCAAGGGGCCGCACTCCCGCAGCCTGGGCCTGGCAACGGCCGTGGCCAACGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCCCACTGCCGCTGGCAGTGAgccaa                                                                                                                        <  4:189903/149‑3 (MQ=255)
                       ccGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGgccgc                                                                                                                        <  8:127741/149‑1 (MQ=255)
                       ccGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGgccgc                                                                                                                        <  8:123435/149‑1 (MQ=255)
                        cgcgGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGgccgcc                                                                                                                       >  6:9191/1‑149 (MQ=255)
                           gggcggCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCca                                                                                                                     >  1:116717/1‑148 (MQ=255)
                            ggcggcCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCcac                                                                                                                    >  2:305928/1‑148 (MQ=255)
                               ggcCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACaaaa                                                                                                                >  6:165636/1‑149 (MQ=255)
                               ggcCAGCGAAGCGGACAGCAAGGCGCCGCACTCCCTCAGCCTGGGCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACaaaa                                                                                                                <  5:233661/149‑1 (MQ=255)
                                gcCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAAc                                                                                                               <  5:9191/149‑1 (MQ=255)
                                 ccAGCGGAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAAc                                                                                                               >  8:290894/1‑148 (MQ=255)
                                         gCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCAcc                                                                                                      >  1:216500/1‑149 (MQ=255)
                                               aGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCAc                                                                                                <  2:205132/149‑1 (MQ=255)
                                                 cAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCAccc                                                                                              <  8:3817/149‑1 (MQ=255)
                                                    ggTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAg                                                                                           <  6:233535/149‑1 (MQ=255)
                                                     gTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAgg                                                                                          <  2:216500/149‑1 (MQ=255)
                                                       gCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATg                                                                                                           <  4:142786/130‑1 (MQ=255)
                                                       gCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATg                                                                                                           >  3:142786/1‑130 (MQ=255)
                                                              tCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCCCACTGCCGCTGGCAGTGAgccaa                                                                                                                        <  4:145138/110‑3 (MQ=255)
                                                              tCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCCCACTGCCGCTGGCAGTGAgccaa                                                                                                                        >  3:145138/1‑108 (MQ=255)
                                                                ccTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGAcc                                                                               >  5:259351/1‑149 (MQ=255)
                                                                ccTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGAcc                                                                               >  6:40832/1‑149 (MQ=255)
                                                                ccTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGAcc                                                                               <  2:305305/149‑1 (MQ=255)
                                                                ccTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGAcc                                                                               >  4:301197/1‑149 (MQ=255)
                                                                 cTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGAccc                                                                              <  2:116717/149‑1 (MQ=255)
                                                                    aGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCgg                                                                            <  5:165636/148‑1 (MQ=255)
                                                                    aGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCgg                                                                            <  7:290894/148‑1 (MQ=255)
                                                                          gTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATg                                                                                                           >  8:40686/1‑111 (MQ=255)
                                                                          gTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATg                                                                                                           <  7:40686/111‑1 (MQ=255)
                                                                            ccTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGCGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCCGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCg                                                                   >  5:261235/1‑149 (MQ=255)
                                                                             cTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCgg                                                                  >  2:134737/1‑149 (MQ=255)
                                                                                gCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCa                                                               >  2:110368/1‑149 (MQ=255)
                                                                                 cAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCAc                                                              >  6:48825/1‑149 (MQ=255)
                                                                                 cAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCAc                                                              <  1:134737/149‑1 (MQ=255)
                                                                                  aaaGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCAcg                                                             >  4:222883/1‑149 (MQ=255)
                                                                                     ggCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCAcgcgc                                                          <  5:48825/149‑1 (MQ=255)
                                                                                            gCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCCCCCCAGGCTGCGCCGCCCCGGTAGCGCCCGGCCACGCGCACCTGtt                                                  >  4:82059/1‑149 (MQ=255)
                                                                                             ccAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCg                                                                                                                                    <  3:129732/67‑1 (MQ=255)
                                                                                             ccAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCg                                                                                                                                    >  4:129732/1‑67 (MQ=255)
                                                                                              cAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGTCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTg                                                >  4:224001/1‑149 (MQ=255)
                                                                                              cAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTg                                                >  3:306036/1‑149 (MQ=255)
                                                                                              cAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTg                                                >  5:96763/1‑149 (MQ=255)
                                                                                               aaaGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTgg                                               <  3:222883/149‑1 (MQ=255)
                                                                                                     ccGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCgcgg                                         <  1:305928/149‑1 (MQ=255)
                                                                                                         aCGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGcc                                     >  5:152345/1‑149 (MQ=255)
                                                                                                             ccGCTACCAGA‑‑GCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCa                                >  1:45370/1‑148 (MQ=255)
                                                                                                              cGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGc                                 >  1:98326/1‑148 (MQ=255)
                                                                                                                   ccAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCaaaa                             <  5:290673/147‑1 (MQ=255)
                                                                                                                            gcgcgcCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCAccc                  >  5:90293/1‑149 (MQ=255)
                                                                                                                                gcCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCctgct              >  7:2384/1‑149 (MQ=255)
                                                                                                                                    ttGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTgctgct           >  5:120812/1‑148 (MQ=255)
                                                                                                                                      gCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCctgct              <  4:299533/143‑1 (MQ=255)
                                                                                                                                      gCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCctgct              >  3:299533/1‑143 (MQ=255)
                                                                                                                                       cAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAg       <  4:306036/149‑1 (MQ=255)
                                                                                                                                        aGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAgg      >  1:292994/1‑149 (MQ=255)
                                                                                                                                            gACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAGGCCtt  <  6:259351/149‑1 (MQ=255)
                                                                                                                                            gACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAGGCCt   <  8:2384/148‑1 (MQ=255)
                                                                                                                                              ggcggCCCACTGCCGCTGGCAGTGAGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTAt                       <  6:298890/123‑1 (MQ=255)
                                                                                                                                              ggcggCCCACTGCCGCTGGCAGTGAGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTAt                       >  5:298890/4‑126 (MQ=255)
                                                                                                                                                 ggCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTCGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAg       <  1:71762/140‑1 (MQ=255)
                                                                                                                                                 ggCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCgcgg                                         <  6:149347/105‑1 (MQ=255)
                                                                                                                                                 ggCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCgcgg                                         >  5:149347/1‑105 (MQ=255)
                                                                                                                                                 ggCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAg       >  2:71762/1‑139 (MQ=255)
                                                                                                                                                  gCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAGGCCt   <  2:273106/142‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GACAGCTGCATGCTGGTGTGCAGCCGCGGGCGGCCAGCGAAGCGGACAGCAAGGTGCCGCACTCCCTCAGCCTGGTCCTGGCAAAGGCCTTGGCCAAAGAGCCGGACGCCCGCTACCAGAGCGCGCGCGCCTTTGCAGTGGACTTGGCTCACTGCCAGCGGCAGTGGGCCGCCCACAAAACCATGGCACCCTTCACCCCAGGCTGCGCCGACCCGGTAGCGCCCGGCCACGCGCGCCTGTTTGGCCGCGGCGCCGAGCAAAAGGCTATCACCCTGCTGCTCAAGGCCTT  >  NC_002947/5980023‑5980311

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: