Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,701,515 C→G 5.9% A649P (GCG→CCG)  kdpB ← K+ transporting ATPase subunit KdpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,701,5150CG 5.9% 275.8 / 5.9 85A649P (GCG→CCG) kdpBK+ transporting ATPase subunit KdpB
Reads supporting (aligned to +/- strand):  ref base C (37/43);  new base G (3/2);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

AATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACACCCAGTTGCGGGTAGAT  >  NC_002947/4701369‑4701662
                                                                                                                                                  |                                                                                                                                                    
aaTGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCTCGTAGATCAGCAGGTTGCGCCGCAACAGGTGcgcc                                                                                                                                                     <  1:61035/148‑1 (MQ=255)
             gATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTg                                                                                                                                      <  6:145204/149‑1 (MQ=255)
              aTAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGc                                                                                                                                     <  2:348149/149‑1 (MQ=255)
              aTAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATAAGCAGGTTGCGCCGAAACAGGTGCGCCGCACTGGCCGCCTg                                                                                                                                      >  7:75051/1‑148 (MQ=255)
               tAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGc                                                                                                                                     <  4:397309/148‑1 (MQ=255)
                aaGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCAc                                                                                                                                   >  8:288777/1‑149 (MQ=255)
                aaGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGGAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCAc                                                                                                                                   >  6:414344/1‑149 (MQ=255)
                  gCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCTGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCAcgc                                                                                                                                 <  4:214694/149‑1 (MQ=255)
                     gTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCa                                                                                                                              >  4:167252/1‑149 (MQ=255)
                        aTGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCAcc                                                                                                                            <  2:211397/148‑1 (MQ=255)
                           gcggcCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCAc                                                                                                                        >  6:69113/1‑149 (MQ=255)
                           gcggcCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCAc                                                                                                                        <  5:69113/149‑1 (MQ=255)
                                   ctcAGACCAGGTTCAGGGCGTTGAGCAGTAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAgcgc                                                                                                                 >  2:184020/1‑148 (MQ=255)
                                   ctcAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAgcgc                                                                                                                 >  6:91118/1‑148 (MQ=255)
                                   ctcAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAgcgc                                                                                                                 >  2:376751/1‑148 (MQ=255)
                                      aGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGCAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGc                                                                                                             <  5:414344/149‑1 (MQ=255)
                                       gACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCg                                                                                                            >  2:268612/1‑149 (MQ=255)
                                       gACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCg                                                                                                            <  5:91118/149‑1 (MQ=255)
                                        aCCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCGTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCgg                                                                                                           <  7:298752/149‑1 (MQ=255)
                                         ccAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCggg                                                                                                          <  1:376751/149‑1 (MQ=255)
                                         ccAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCggg                                                                                                          <  3:207974/149‑1 (MQ=255)
                                           aGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGAt                                                                                                        >  3:211901/1‑149 (MQ=255)
                                            ggTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATc                                                                                                       <  4:211901/149‑1 (MQ=255)
                                             gTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCAcccc                                                                                                                          >  6:374377/1‑129 (MQ=255)
                                             gTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCAcccc                                                                                                                          <  5:374377/129‑1 (MQ=255)
                                             gTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCa                                                                                                      <  7:276308/149‑1 (MQ=255)
                                              ttCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAg                                                                                                     >  7:11092/1‑149 (MQ=255)
                                                 aGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCa                                                                                                   <  5:376875/148‑1 (MQ=255)
                                                    gCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTg                                                                                                                                              <  1:204916/102‑1 (MQ=255)
                                                    gCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTg                                                                                                                                              >  2:204916/1‑102 (MQ=255)
                                                      gTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAg                                                                                                                    >  2:363780/1‑126 (MQ=255)
                                                      gTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAg                                                                                                                    <  1:363780/126‑1 (MQ=255)
                                                         gAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGa                                                                                          >  2:171554/1‑149 (MQ=255)
                                                             agcagGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAg                                                                                      >  5:181122/1‑149 (MQ=255)
                                                              gcagGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGgcc                                                                                                                                           >  1:196609/1‑95 (MQ=255)
                                                              gcagGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGgcc                                                                                                                                           <  2:196609/95‑1 (MQ=255)
                                                              gcagGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAgc                                                                                     <  6:181122/149‑1 (MQ=255)
                                                               cagGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCggg                                                                                                          <  3:167930/127‑1 (MQ=255)
                                                               cagGTCGATCAGCTTGATCCCGGCAAACGGCACAACAATGCAACCCAGCCAGTAGATCAGCAGGTTGCGCAGCAACAGGTGCGCCGCAATGGCCGCCTGAACGCGCACCCAACGCAGCGCCAGCggg                                                                                                          >  4:167930/1‑127 (MQ=255)
                                                                  gTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGt                                                                                 >  3:397938/1‑149 (MQ=255)
                                                                  gTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGt                                                                                 >  2:121060/1‑149 (MQ=255)
                                                                     gATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGa                                                                              <  8:11092/149‑1 (MQ=255)
                                                                      aTCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGAt                                                                                         <  1:57762/137‑1 (MQ=255)
                                                                      aTCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGAt                                                                                         >  2:57762/1‑137 (MQ=255)
                                                                          gCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAacac                                                                          >  5:127431/1‑148 (MQ=255)
                                                                              gATCCCTGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGAAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATc                                                                      <  2:6849/148‑1 (MQ=255)
                                                                                  ccGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGa                                                                 >  3:64654/1‑149 (MQ=255)
                                                                                  ccGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGa                                                                 >  7:374900/1‑149 (MQ=255)
                                                                                   cGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGa                                                                 >  5:188448/1‑148 (MQ=255)
                                                                                   cGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGAc                                                                <  6:188448/149‑1 (MQ=255)
                                                                                   cGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGAc                                                                <  1:121060/149‑1 (MQ=255)
                                                                                   cGGCAAACGGCACCACAATGCCACACAGCACGTAGATCATCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCTCCAGCGGGATCAGCACGATAATCATCACCGCGATGAACACGATCGCCGa                                                                 >  4:119906/1‑148 (MQ=255)
                                                                                    ggCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACa                                                               <  1:184020/149‑1 (MQ=255)
                                                                                     gCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAg                                                              <  4:64654/149‑1 (MQ=255)
                                                                                       aaaCGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGa                                                            <  2:26822/149‑1 (MQ=255)
                                                                                         aCGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCACACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGAACAGCGCGATGAACACGATCGCCGACAGGATc                                                          >  5:156935/1‑149 (MQ=255)
                                                                                          cGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATcg                                                         <  2:222533/149‑1 (MQ=255)
                                                                                           ggCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCACCAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATcgc                                                        >  3:54883/1‑149 (MQ=255)
                                                                                            gcaagacaATGCCACCCAGCCCGTAGATGAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCACACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATcgcg                                                       <  6:156935/144‑1 (MQ=255)
                                                                                              accacAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGct                                                     >  7:216255/1‑149 (MQ=255)
                                                                                                cacaATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGctct                                                   >  6:383799/1‑149 (MQ=255)
                                                                                                   aaTGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCg                                                >  6:305870/1‑149 (MQ=255)
                                                                                                        cACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAg                                                              >  1:291114/1‑131 (MQ=255)
                                                                                                        cACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAg                                                              <  2:291114/130‑1 (MQ=255)
                                                                                                        cACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTg                                           >  8:42447/1‑149 (MQ=255)
                                                                                                        cACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTg                                           <  8:374900/149‑1 (MQ=255)
                                                                                                        cACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTg                                           <  4:397938/149‑1 (MQ=255)
                                                                                                          ccgaGCCCGTAGATCAGCTGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCTCGCGCACCCCACGCAGGGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGTCAGGATCGCGCTCTGCGGGCTGGc                                         <  4:54883/146‑1 (MQ=255)
                                                                                                              gCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATcgcg                                                       >  7:274219/1‑131 (MQ=255)
                                                                                                              gCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATcgcg                                                       <  8:274219/131‑1 (MQ=255)
                                                                                                               cccGTAGATCAGCATGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGTATCAGCACGATAATGATCAGCGCGTTGAACACGATCTCCGACAGGATCGCGCTCTGCGGGCTGGCAAGtt                                    <  6:271023/149‑1 (MQ=255)
                                                                                                               cccGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGtt                                    <  6:127431/149‑1 (MQ=255)
                                                                                                                 cGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGc                                  <  8:216255/149‑1 (MQ=255)
                                                                                                                  gTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTgg                                          <  5:253750/140‑1 (MQ=255)
                                                                                                                  gTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTgg                                          >  6:253750/1‑140 (MQ=255)
                                                                                                                     gATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCa                              <  7:288777/149‑1 (MQ=255)
                                                                                                                           cagGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGa                        <  5:305870/149‑1 (MQ=255)
                                                                                                                             cagtGCGGCGCACCTGTTGCGGCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACa                                                               >  6:231909/4‑108 (MQ=255)
                                                                                                                             cagtGCGGCGCACCTGTTGCGGCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACa                                                               <  5:231909/105‑1 (MQ=255)
                                                                                                                                tGCGGCGCACCTGTTGCGGCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGcaca                   >  7:235057/1‑149 (MQ=255)
                                                                                                                                tGCGGCGCACCTGTTGCGGCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGcaca                   >  8:102353/1‑149 (MQ=255)
                                                                                                                                   gccgcAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACAcc                 >  4:169929/1‑148 (MQ=255)
                                                                                                                                      gcAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACACCCAGt             <  3:167252/149‑1 (MQ=255)
                                                                                                                                        aaCAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACACCCAg              <  2:402928/146‑1 (MQ=255)
                                                                                                                                        aaCAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACACCCAg              >  1:402928/1‑146 (MQ=255)
                                                                                                                                             ttGCGGCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACACCCAGTTGCgg        <  8:235057/146‑1 (MQ=255)
                                                                                                                                                cgccgcACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACACCCAGTTGCGGGTAg    <  3:119906/148‑1 (MQ=255)
                                                                                                                                                 gccgcACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACACCCAGTTGCGGGTAGAt  >  2:30047/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AATGCCGGGCGTAGATAAGCAGTCATGGCGGCCTCCTCAGACCAGGTTCAGGGCGTTGAGCAGCAGGTCGATCAGCTTGATCCCGGCAAACGGCACCACAATGCCACCCAGCCCGTAGATCAGCAGGTTGCGCCGCAACAGGTGCGCCGCACTGGCCGCCTGCACGCGCACCCCACGCAGCGCCAGCGGGATCAGCACGATAATGATCAGCGCGTTGAACACGATCGCCGACAGGATCGCGCTCTGCGGGCTGGCCAGTTGCATCAGGTTGAGCACACCCAGTTGCGGGTAGAT  >  NC_002947/4701369‑4701662

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: