Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,781,385 G→C 9.1% intergenic (+75/‑68) PP_3338 → / → PP_3339 cytochrome bd‑type quinol oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,781,3850GC 9.1% 367.2 / 3.1 146intergenic (+75/‑68)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (64/67);  new base C (3/12);  total (67/79)
Fisher's exact test for biased strand distribution p-value = 5.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

GCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGACC  >  NC_002947/3781236‑3781534
                                                                                                                                                     |                                                                                                                                                      
gCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGTCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcg                                                                                                                                                        <  2:5434944/150‑1 (MQ=255)
 ccGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcgc                                                                                                                                                       <  2:2382866/150‑1 (MQ=255)
   gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcac                                                                                                                                                     >  2:3260515/1‑150 (MQ=255)
   gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAGGCGcac                                                                                                                                                     >  2:3260516/1‑150 (MQ=255)
           gccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGACAAAAATGCGCACCACCGAAc                                                                                                                                             >  1:1765047/1‑150 (MQ=255)
                  agCCACGGTTTCTTGCGGGTAGGCACCCGTTTCATCCTGCCCTTTTGGCTTGGCTGATCCGACTTTTCTTGAGGGCCGCCGGCATGCGCAGCAATTGCCCCCTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACa                                                                                                                                      <  2:2471358/150‑1 (MQ=255)
                    ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACCAAAAGGCGCACCACGGAACATGTACaac                                                                                                                                    >  2:2883964/1‑148 (MQ=255)
                    ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCAAAAACGCCCCCCACGGAACGTTTAAAGc                                                                                                                                   >  2:491908/1‑150 (MQ=255)
                    ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAAGCGCACCACCGCACATGTaaatc                                                                                                                                    >  1:2187213/1‑146 (MQ=255)
                      aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcgc                                                                                                                                                       <  2:4256533/129‑1 (MQ=255)
                      aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcgc                                                                                                                                                       >  1:4256533/1‑129 (MQ=255)
                       aGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGACGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAAGCGCACCACCGAAAATGTACAccccc                                                                                                                                >  1:4184922/1‑149 (MQ=255)
                           ttaaTGCTGGTAGCCACAATGTACATCCTGCCCAGTTCGCTTGGCTCATCCGACTGTTCATGATGCCCGCCGCCATGCGCACCAATTGACCCCTCGCCTTGTTGCGCCTGCGACAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAAc                                                                                                                            <  2:2274822/147‑1 (MQ=255)
                                gCTGTTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGAATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGAACCACGGAACTTGTACAGCCCGCAACGGATg                                                                                                                       <  2:901322/150‑1 (MQ=255)
                                 cTGGTCGGCACACTGTCCAGCCTGCCCCGTTCGCTTGGCTGACCCGACTGTTCATGCGGGCCGCGGGCTGGCGCAGCCATTGAGCCCGCGCGGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                                      <  2:3014815/150‑1 (MQ=255)
                                 cTGGTAGGCACACTGTTCACCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCCCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                                      <  2:819285/150‑1 (MQ=255)
                                  tGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGAACATGTACAGCCCCCAACCGATGCg                                                                                                                     >  2:5603248/1‑150 (MQ=255)
                                   ggTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCgg                                                                                                                    <  2:2799692/150‑1 (MQ=255)
                                   ggTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACGTGTACAGCCCGCAACCGATGCgg                                                                                                                    >  2:3278412/1‑150 (MQ=255)
                                         caccCTGGTCATCCTGCCCTGTGCGCTTGGCTGATCCGACTTTTCAGGATGGCCGCCGGAATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCCACAACACGGCGCACCACGGAACTTGTACAGCCCGCATCCGTTGCgggctgta                                                                                                              <  2:4088891/150‑7 (MQ=255)
                                         cacaGTGTTCATCCTGCCCATTTCGGTTGGGTGATCCGTCTGTTCATTTTGGCCGCCGGCTTGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCCGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAAt                                                                                                              <  2:1659732/150‑1 (MQ=255)
                                          acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGGCAAAAAGGCGCACCACGGAACATGTACAGCCCGCAACGGCTGCGGCTTaaag                                                                                                             >  2:5421280/1‑148 (MQ=255)
                                          acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGGCAAAAAGGCGCACCACGGAACGTTTAAAGCCCGCAACCGAGGCGGCGTaaag                                                                                                             >  2:1611648/1‑148 (MQ=255)
                                            acTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGCCCAAAACGCGCACCACGGAACATGTACAGCCCCCCACCGGAGCGGCGTAATgcg                                                                                                           >  2:3955871/1‑150 (MQ=255)
                                             cTGTTCATCGTGCCCTGTTCGGTTGGCTGATCCGCCTTTTCATGTTGGCCGCGGGCTTGCGCTGCCTTTGCGCCCTGGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgc                                                                                                          <  2:2468419/150‑1 (MQ=255)
                                             cTGTTCATCCTGCCCCTTTCGCTGGGCTGAGCCGACGGTGCTTGAGGGCCGCCGGCCGGCGCAGCAATGGCGCCGTGGCCGTGGGGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgc                                                                                                          <  2:4745410/150‑1 (MQ=255)
                                             cTGTTCATCCTGCCCAGTTCGCTTGGCTTATCCGACGGTTCTTGTTGGCCGCCGGCATGCGCAGCAATTGAGCCCTCGCGGAGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgc                                                                                                          <  2:4850023/150‑1 (MQ=255)
                                              tGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGGACGTGTACAGCCCGCAACGGATGCGGCGTAATgcgcg                                                                                                         >  2:4219374/1‑150 (MQ=255)
                                               gTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgcgc                                                                                                        >  2:5520325/1‑150 (MQ=255)
                                                    tCCTGCTCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCAAAAACGCGCACCCCCGGACATGTACAACCCGCCACCGAGGCGGCGTAATGCGCGcccaat                                                                                                   >  2:2086869/1‑146 (MQ=255)
                                                    tCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAACATGTACAGCCCGCAACCGATGCGGCTTAATGCGCGCCaaaa                                                                                                   >  2:3568214/1‑147 (MQ=255)
                                                    tCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCAGCGACAAAAAAGCCCACCACCGCACAGGTACAGCCCCCAACCCAAGCGGCGTAAAGCGCGcccaaa                                                                                                   >  2:1962480/1‑146 (MQ=255)
                                                      cTGCCCTTTTCGCTTGGCTTACCCTCCTGTTCATGTTGGCCGCCGGCATGCGCAGCCATTGAGCCCCCGCCGTGGTGCGCCTGCCAAAAAACGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTc                                                                                                 <  2:1404624/150‑1 (MQ=255)
                                                      cTGCCCCGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCAATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTc                                                                                                 <  2:1832005/150‑1 (MQ=255)
                                                       tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCGGCGCCAAAAAGGCGCCCCACGGCACGTGTACCGCCCCCCCCCGCAGCGGCGTAAGGCGCGCCATAttgg                                                                                                >  2:537762/1‑148 (MQ=255)
                                                       tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCAAAAAGGCGCCCCCCCGCACCTGTACAGCCCGCCACCGATGCGGCGGAATGCGCGcccgattag                                                                                                >  2:3351973/1‑143 (MQ=255)
                                                          cctgTTCGGTTGGCTGACCCTCCTGTTATTTTTGCCCGCCGGCATGCGCCGCAATTGCGCCCTGGCGGTGGTGCGCCTCCTAAAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGTGGCTTAATGCGCGCCATATTCGCTc                                                                                             <  1:1225710/147‑1 (MQ=255)
                                                            tgTTGGCTTGGCTGAGCCGCGTGTTCTTGTGGGCCGCGGGACTGGGCCCCCTTTGCGCCCTGGCGGTGGGGCGCCTGCCAAAAAATGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCgg                                                                                           <  1:2684419/149‑1 (MQ=255)
                                                            tgTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCgg                                                                                           <  2:3943409/149‑1 (MQ=255)
                                                                 gCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAAATTGTTCAACCCCCAACGGGTGCGGCTTAATTCGCGCCACGTTCGGGCGGGTgtc                                                                                      >  2:4390704/1‑148 (MQ=255)
                                                                    tGGCTGATCCGACTGTTCATGCGGGCCGCCGGAATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCa                                                                                   <  2:4321939/150‑1 (MQ=255)
                                                                      gCTGATCCGACTGTTCTTGAGGGCCGCCGGCATGCGCAGCCATTGAGCCCCCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATg                                                                                 <  2:5370523/150‑1 (MQ=255)
                                                                      gCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATg                                                                                 <  2:5601251/150‑1 (MQ=255)
                                                                          aTCCTAGTGTTCATGTGGGCCGCCGGCATGCGCAGCCATTGCGCCCCCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGTCTCATGCttt                                                                             <  2:187681/150‑1 (MQ=255)
                                                                          aTCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCttt                                                                             <  2:814872/150‑1 (MQ=255)
                                                                             cGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                                          <  2:5007313/150‑1 (MQ=255)
                                                                             cGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                                          <  2:5007409/150‑1 (MQ=255)
                                                                             cGACTGTTATTTTTGGCGGCCGGATGGGGCAGCAATTGCGCCCTCGCGGTGGTGCGCCTGCGAAAAATTGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                                          <  1:893608/150‑1 (MQ=255)
                                                                              gACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGGCAAAAAAGCGCCCCACCGGACATTTAAAACCCGCAACCGATTAGGCGTTATTCGCGCCAAAATAGCTCGGGTGACTCAAGCTTTgacg                                                                         >  2:564522/1‑150 (MQ=255)
                                                                                  gTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGCCAAAAAGGCGCACCACCGCACGTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCCAATTCGCGCGGGGGACTCATGCTTTcacgtgca                                                                     >  2:3657338/1‑150 (MQ=255)
                                                                                      atgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGCCAAAAAGGCGCCCCACCGGACGTGTTCAGCCCGCAACCGCTGCGGCGGAATTCGCGCCCAAATAGCGCGGGTGACACAATCTTTcacgtgcacctg                                                                 >  2:3812795/1‑150 (MQ=255)
                                                                                      atgaggGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                                                 <  2:4403852/150‑1 (MQ=255)
                                                                                       tgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCAAAAAGGCCCACCACCGGACTTGTACAGCCCCCAACCGATGCGGCTTAATGCGCGCCCAAATCGCTCGGGTGACTCATGCTTTCacgtgcacgtgc                                                                >  2:3816853/1‑150 (MQ=255)
                                                                                              cgccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCAAAAAGGCGCACCACCGCCCGAGTAAACCCCCCAACCGCTTCGGGGTAATGCGCGCCCGAGTCGCGCGGGGGGCGCATGCTTTCACGGGTAAGTGcacagagt                                                         >  2:1075483/1‑146 (MQ=255)
                                                                                              cgccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGGCAAAAACGCGCACCACCGACCGTGTAAACCCCCCCCACGGTGCGGCGTAATTCGCGCCCTAGTAGCGCGGGGGGCACATGCTTTAACCGGTAAGTGcacaaagt                                                         >  2:1075442/1‑146 (MQ=255)
                                                                                               gccgGCTTGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                        <  1:5076618/150‑1 (MQ=255)
                                                                                               gccgGCATGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                        <  1:2731651/150‑1 (MQ=255)
                                                                                               gccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                        <  2:1260689/150‑1 (MQ=255)
                                                                                                 cgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                                      <  2:4673357/150‑1 (MQ=255)
                                                                                                 cgGCAGGCGCAGCCATTGCGCCCTCGCCGTGGTGGGCCTTCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                                      <  2:1815857/150‑1 (MQ=255)
                                                                                                   gcttGCGCCGCCTTTGCGCCCTCGCCGTGGGGCGCCTGCGACAAAAAGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                                                    <  2:1928016/147‑1 (MQ=255)
                                                                                                   gCATGCGCTGCAATTTTGCCTTGGCGGTGGGGCGCCTGCCAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                                                    <  1:4812986/150‑1 (MQ=255)
                                                                                                   gCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                                                    <  2:4526579/150‑1 (MQ=255)
                                                                                                    cATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAATGCGCACCACGGAACTTTTACAGCCCGCAACCGATGCGGCGTAATGCGCGCCAAATTAGCGCGGGTGTCTCATGCTTTAACGTGCACGTGCACAACATTGTCTg                                                   >  1:5628309/1‑150 (MQ=255)
                                                                                                             ccATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgc                                          <  2:1765047/150‑1 (MQ=255)
                                                                                                             ccATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgc                                          <  2:3686270/150‑1 (MQ=255)
                                                                                                               aTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACCGATGCGGCGTAAAGCGCGCCAAAATCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCGTTGTCTGTTCgcctgcca                                        >  2:5627458/1‑149 (MQ=255)
                                                                                                               aTTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAACGCGCACCACGGCACTTGTACAGCCCCCAAACGCTTCGGCGGAATTCGCGCCAAATTCGCGCGGGTGACGCATGCGGTAACCTGCACGTGCACAACATTGTCGGGTCgcctgccg                                        >  2:4025838/1‑149 (MQ=255)
                                                                                                                   aGCCCTCGCGGTGGGGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc                                    <  2:2849846/150‑1 (MQ=255)
                                                                                                                   aGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc                                    <  1:5144774/150‑1 (MQ=255)
                                                                                                                        tCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt                               <  1:5044508/150‑1 (MQ=255)
                                                                                                                            ggtggtgCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTc                           <  1:3098642/149‑1 (MQ=255)
                                                                                                                              tggtgCGCCTGCGACAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCaa                         <  1:3196240/150‑1 (MQ=255)
                                                                                                                               ggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc                        <  1:5603248/150‑1 (MQ=255)
                                                                                                                                 tgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                                <  1:895852/140‑1 (MQ=255)
                                                                                                                                   cgcCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgc                    >  2:3540350/1‑150 (MQ=255)
                                                                                                                                     ccTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgctg                  <  2:1874888/150‑1 (MQ=255)
                                                                                                                                     ccTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgctg                  <  2:4709040/150‑1 (MQ=255)
                                                                                                                                      cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg                 >  1:3374797/1‑150 (MQ=255)
                                                                                                                                        gCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCg               >  1:1727819/1‑150 (MQ=255)
                                                                                                                                         cGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGa              <  1:1294519/150‑1 (MQ=255)
                                                                                                                                         cGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGa              >  2:1742805/1‑150 (MQ=255)
                                                                                                                                          gACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGAt             >  1:1541823/1‑150 (MQ=255)
                                                                                                                                          gACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGAt             >  1:369696/1‑150 (MQ=255)
                                                                                                                                           aCAAAATGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATg            <  1:2421136/150‑1 (MQ=255)
                                                                                                                                           aCAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATg            >  2:1561191/1‑150 (MQ=255)
                                                                                                                                            cAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGc           >  1:260910/1‑150 (MQ=255)
                                                                                                                                             caacaGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc          <  1:5392352/146‑1 (MQ=255)
                                                                                                                                             aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc          >  1:1834060/1‑150 (MQ=255)
                                                                                                                                             aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc          >  2:4631415/1‑150 (MQ=255)
                                                                                                                                             aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc          >  1:4297575/1‑150 (MQ=255)
                                                                                                                                             aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc          >  2:2913771/1‑150 (MQ=255)
                                                                                                                                               aacggCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGc        <  1:4942819/147‑1 (MQ=255)
                                                                                                                                               aaaGCCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGc        <  1:1944054/150‑1 (MQ=255)
                                                                                                                                                 aGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCtt      <  2:2563286/150‑1 (MQ=255)
                                                                                                                                                 aGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCtt      <  2:2563208/150‑1 (MQ=255)
                                                                                                                                                 aGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCtt      <  1:150067/150‑1 (MQ=255)
                                                                                                                                                 aGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCtt      <  1:5392106/150‑1 (MQ=255)
                                                                                                                                                     gcACCACGGCCCTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATAGTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCGTTGTCTGTTCGCCTGCCTCGCCGTGCTGTGCAACCTGCTGGCGATGCCGCTGGAcc  >  2:3928214/1‑150 (MQ=255)
                                                                                                                                                     gcACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGAcc  >  2:1322702/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGACC  >  NC_002947/3781236‑3781534

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: