Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,287,469 C→G 8.4% A89A (GCG→GCC)  PP_3758 ← two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,287,4690CG 8.4% 201.8 / 7.3 69A89A (GCG→GCC) PP_3758two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (2/4);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAGGCAGCCT  >  NC_002947/4287326‑4287617
                                                                                                                                               |                                                                                                                                                    
tGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTAGCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCcgccgc                                                                                                                                                 >  2:125692/1‑149 (MQ=255)
  ccGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCgcgc                                                                                                                                               >  5:189253/1‑149 (MQ=255)
  ccGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCgcgc                                                                                                                                               >  3:346015/1‑149 (MQ=255)
    ggCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGGGCcgcgaga                                                                                                                                             <  7:145644/149‑4 (MQ=255)
        aaGCACCTTGCGCTGGCGGAACAGGTCGACGAACCCCGAGAACTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCCTCTCCCGTCCCGCcacgaggccga                                                                                                                                         >  6:251561/1‑142 (MQ=255)
        aaGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACg                                                                                                                                          >  7:97158/1‑148 (MQ=255)
        aaGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACg                                                                                                                                          >  5:96086/1‑148 (MQ=255)
         aGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGaa                                                                                                                                        <  8:97158/149‑1 (MQ=255)
           cACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGa                                                                                                                                      <  6:96086/149‑1 (MQ=255)
              cTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGa                                                                                                                                            <  1:303333/140‑1 (MQ=255)
              cTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGa                                                                                                                                            >  2:303333/1‑140 (MQ=255)
                 gcgcTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGACATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATc                                                                                                                                 <  4:181590/148‑1 (MQ=255)
                  cgcTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCg                                                                                                                                <  2:344133/148‑1 (MQ=255)
                       gCGGAACAGGTCGACTAACACCGAGACCTTGCTTTTTACCGCCAGTGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGAt                                                                                                                            <  5:251561/147‑1 (MQ=255)
                       gCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGAtgt                                                                                                                          >  6:246133/1‑149 (MQ=255)
                         ggAACAGGTCGACGAACACCGAGACCTTGCTTTTTAGCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGc                                                                                                                        <  1:125692/149‑1 (MQ=255)
                         ggAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGAtgtgg                                                                                                                        >  3:257114/1‑148 (MQ=255)
                         ggAACAGGTCGACCAACACCGAGACCTTGCTTTTTACCGCCACGTTGTCGAGCGGCTTGTGAAGAAAGTCCACGCCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTAGCGGCCCGCCGCGCTCACGAAGACGATCGGGATGTGc                                                                                                                        >  7:139926/1‑149 (MQ=255)
                              aGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGACGTGCTTGGt                                                                                                                   >  4:213726/1‑149 (MQ=255)
                                  cgacgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGGGCcgcgaga                                                                                                                                             >  7:120543/1‑116 (MQ=255)
                                  cgacgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGGGCcgcgaga                                                                                                                                             <  8:120543/119‑4 (MQ=255)
                                   gacgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTc                                                                                                                  >  3:139253/1‑145 (MQ=255)
                                   gacgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTc                                                                                                                  <  4:139253/145‑1 (MQ=255)
                                   gacgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGACCTTGAAGGAATAATTCATCTCGCGGCCCGACGCGCTGACGAAGACGATCGGGATGTGCTTGGTATTGt                                                                                                              >  4:270172/1‑149 (MQ=255)
                                    acgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTGTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTc                                                                                                             <  1:112579/149‑1 (MQ=255)
                                    acgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCGCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTc                                                                                                             <  1:235643/149‑1 (MQ=255)
                                    acgaACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTc                                                                                                             <  4:74950/149‑1 (MQ=255)
                                                aCCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCAt                                                                                                 <  6:189253/149‑1 (MQ=255)
                                                  cTTGCTTTTTACCGCCAGGGTGTCGAGCGGATTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTATTGTAGGTACCGAGCAt                                                                                                 >  6:340525/1‑147 (MQ=255)
                                                       tttttACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGt                                                                                                              >  6:260484/1‑129 (MQ=255)
                                                       tttttACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGt                                                                                                              <  5:260484/129‑1 (MQ=255)
                                                       tttttACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCg                                                                                          >  2:289136/1‑149 (MQ=255)
                                                         tttACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCggg                                                                                                                              <  6:181048/111‑1 (MQ=255)
                                                         tttACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCggg                                                                                                                              >  5:181048/1‑111 (MQ=255)
                                                          ttACGGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCTTAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTTCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGc                                                                                        <  8:139926/148‑1 (MQ=255)
                                                            aCCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACCGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCg                                                                                          >  3:87347/1‑144 (MQ=255)
                                                            aCCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACAGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCg                                                                                          <  4:87347/144‑1 (MQ=255)
                                                               gCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTc                                                                                  <  6:326438/149‑1 (MQ=255)
                                                                 cAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCg                                                                                 >  3:39979/1‑148 (MQ=255)
                                                                   gggTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTc                                                                                                             >  2:358841/1‑118 (MQ=255)
                                                                                cTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCg                                                                 <  1:231772/149‑1 (MQ=255)
                                                                                    tgCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCAt                                                             <  1:289136/149‑1 (MQ=255)
                                                                                      cccccactCCCCGGCCCCCCTCCCGCACCCCATCACCCCATAATTCATCTCGCGGCCCGCCGCGCAGACGAAGACGATCCGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCt                                                           <  4:257114/142‑1 (MQ=255)
                                                                                      cAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCt                                                           >  2:15579/1‑149 (MQ=255)
                                                                                        gAAAGTCCACGGCCCCGCTCTCGTAGCCATTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGAAGAAGACGATCGGGATGTGATTGGTCTTGTGGGTACCGCGCAGCAGGTAGGCAAGGTGGAATCAGGTCATGCCCGTCATcgag                                                         >  8:288588/1‑146 (MQ=255)
                                                                                         aaaGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGca                                                        >  6:264086/1‑149 (MQ=255)
                                                                                           aGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGcaca                                                      >  8:25309/1‑149 (MQ=255)
                                                                                            gTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACAt                                                     >  1:186625/1‑149 (MQ=255)
                                                                                               cACGGCCCCGCTCTCGTAGCCCTTGAAGTCATAATTCATCTCGCGTCCCGCCGCGCTGACGAAGACGATCGTGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCga                                                  <  3:270172/149‑1 (MQ=255)
                                                                                               cACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTc                                                                                  <  3:294769/117‑1 (MQ=255)
                                                                                               cACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTc                                                                                  >  4:294769/1‑117 (MQ=255)
                                                                                                aCGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCgag                                                 <  2:186625/149‑1 (MQ=255)
                                                                                                  ggCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGaa                                               <  5:246133/149‑1 (MQ=255)
                                                                                                        gCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGc                                           >  5:132873/1‑147 (MQ=255)
                                                                                                             cGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGc                                     >  1:85644/1‑148 (MQ=255)
                                                                                                               tAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGaa                                  >  5:308361/1‑149 (MQ=255)
                                                                                                               tAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGa                                   >  3:184094/1‑148 (MQ=255)
                                                                                                                   ccTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGccc                                                                  <  6:9086/113‑1 (MQ=255)
                                                                                                                   ccTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGccc                                                                  >  5:9086/1‑113 (MQ=255)
                                                                                                                       gAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGACGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGc                           <  3:213726/148‑1 (MQ=255)
                                                                                                                        aaGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTc                         <  6:132873/149‑1 (MQ=255)
                                                                                                                          ggCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCca                       <  7:25309/149‑1 (MQ=255)
                                                                                                                             aTAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCcagca                    <  5:264086/149‑1 (MQ=255)
                                                                                                                             aTAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCcagca                    <  2:85644/149‑1 (MQ=255)
                                                                                                                               aaTTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCg                                                                           <  7:26080/92‑1 (MQ=255)
                                                                                                                               aaTTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCg                                                                           >  8:26080/1‑92 (MQ=255)
                                                                                                                                   cATCTGGCGGGAGGCCGCGCTGTCGAAGACGATGGGGTTGTGCTTGGTGTTGTCGGTACCGCGCATCAGCTGGGCCAGTTCGTAACCGTTCATGCCCGGCATCTGCACATCGTGAATGGGCAGTGCGAACTCGTGCTCCAGCAGCAAc               <  7:273621/148‑1 (MQ=255)
                                                                                                                                    tcagcgcgGCGGGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAAc               >  7:268231/5‑147 (MQ=255)
                                                                                                                                    tcagcgcgGCGGGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAAc               <  8:268231/143‑1 (MQ=255)
                                                                                                                                         tcgGCCCGCCGCGCTTACGAAGACGATCTTGATGTGCTTGTTCTTGTCGGTACCGCTCATCAGCTCGGACAGTTCGAAACAGTTCATGCCCGGCATCTGCACATCGAGAATTGCCAGTGCGAACTCGTGCTCCAGCAGCAACGAc            <  5:340525/144‑1 (MQ=255)
                                                                                                                                               cgccgcGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCg                              >  4:48310/1‑121 (MQ=255)
                                                                                                                                               cgccgcGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCg                              <  3:48310/121‑1 (MQ=255)
                                                                                                                                               cgccgcGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAGGCAGCCt  >  5:264864/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGACGAAGACGATCGGGATGTGCTTGGTCTTGTCGGTACCGCGCATCAGCTCGGCCAGTTCGAAACCGTTCATGCCCGGCATCTGCACATCGAGAATGGCCAGTGCGAACTCGTGCTCCAGCAGCAACGACAAGGCAGCCT  >  NC_002947/4287326‑4287617

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