Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,746,541 G→C 10.0% intergenic (+32/+67) PP_5042 → / ← PP_5043 hypothetical protein/PhoPQ‑activated pathogenicity‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,746,5410GC10.0% 290.7 / 18.5 98intergenic (+32/+67)PP_5042/PP_5043hypothetical protein/PhoPQ‑activated pathogenicity‑like protein
Reads supporting (aligned to +/- strand):  ref base G (44/43);  new base C (5/5);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

TAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCT  >  NC_002947/5746395‑5746685
                                                                                                                                                  |                                                                                                                                                 
tAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgg                                                                                                                                                >  1:875861/1‑150 (MQ=255)
     gCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGggcggc                                                                                                                                           <  1:1043856/150‑1 (MQ=255)
       cGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                        <  1:804173/150‑1 (MQ=255)
          cTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCCAGGGTGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGGGCCGTGCGGCCGGCGCCCg                                                                                                                                     >  2:1021554/1‑150 (MQ=255)
           tgtgCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                    >  2:680542/1‑150 (MQ=255)
               cACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                    <  2:152987/146‑1 (MQ=255)
               cACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                    >  1:152987/1‑146 (MQ=255)
                 ccATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGAGGGCGGCGCTCGATTCtgt                                                                                                                              >  1:539501/1‑150 (MQ=255)
                  cATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCtgtg                                                                                                                             >  1:863712/1‑150 (MQ=255)
                   aTTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTgg                                                                                                                            <  2:1573695/150‑1 (MQ=255)
                    tttATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGc                                                                                                                           >  2:1175870/1‑150 (MQ=255)
                          aCAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGcc                                                                                                                     >  2:1571968/1‑150 (MQ=255)
                           cAAGCCCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCg                                                                                                                    <  1:980647/150‑1 (MQ=255)
                              gcgcCCGGATACCTGCCGCAACCACCACAAGGTCGGGCCACGGCCGGGGCATTGCGCCTACAAGCCCAAAGTGGTTGGGCGTTGAGTCTTGCGCCGGCGCCTAAATCTAGCGCAGCGCGGGAGGCGCTCGCTTCTGTGAAGAGGCAGAGc                                                                                                                 >  2:722664/1‑150 (MQ=255)
                                acacGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCTCCGGGGCATAAATCGAGCGACGCTCGGGAGGCGCTCGCTTCGATGGCGATGCCGAGCTg                                                                                                               >  2:269328/4‑150 (MQ=255)
                                 cccGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTgg                                                                                                              <  1:680542/150‑1 (MQ=255)
                                 cccGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTgg                                                                                                              <  2:875861/150‑1 (MQ=255)
                                    ggATACCAGACGCAACCCCCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGACATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCa                                                                                                           <  2:1086247/150‑1 (MQ=255)
                                          cTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGTTGCCGAGCTGGGCATGAAGc                                                                                                     >  2:371249/1‑150 (MQ=255)
                                          cTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGc                                                                                                     >  1:1825425/1‑150 (MQ=255)
                                            gCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCAc                                                                                                   <  1:1021554/150‑1 (MQ=255)
                                             ccGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCGCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACt                                                                                                  <  1:1893495/150‑1 (MQ=255)
                                               gCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGa                                                                                                >  1:194443/1‑150 (MQ=255)
                                                cacccaccCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGaa                                                                                               <  2:2149310/147‑1 (MQ=255)
                                                   ccaccCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGGGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAAccc                                                                                            >  2:1991181/1‑150 (MQ=255)
                                                    caccCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCa                                                                                           >  2:1634742/1‑150 (MQ=255)
                                                      ccccAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                        <  1:1563670/103‑1 (MQ=255)
                                                      ccccAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCAcc                                                                                         >  1:540147/1‑150 (MQ=255)
                                                            ggTCGGGCCACGGCCGGGGTATTGCGCTTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTccc                                                                                   <  2:1304972/150‑1 (MQ=255)
                                                             gTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCa                                                                                  >  2:806383/1‑150 (MQ=255)
                                                                gggCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCAttt                                                                               >  1:69093/1‑150 (MQ=255)
                                                                  gCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTca                                                                             >  2:618059/1‑150 (MQ=255)
                                                                  gCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTca                                                                             >  1:1821044/1‑150 (MQ=255)
                                                                   ccACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTc                                                                                     <  2:1603104/141‑1 (MQ=255)
                                                                   ccACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTc                                                                                     >  1:1603104/1‑141 (MQ=255)
                                                                    cACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTcagc                                                                           <  2:194443/150‑1 (MQ=255)
                                                                    cACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTcagc                                                                           <  2:304164/150‑1 (MQ=255)
                                                                     aCGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCa                                                                                                    <  2:113421/124‑1 (MQ=255)
                                                                     aCGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCa                                                                                                    >  1:113421/1‑124 (MQ=255)
                                                                     aCGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTcagcc                                                                          >  1:356944/1‑150 (MQ=255)
                                                                        gCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTcagccagc                                                                       <  1:618059/150‑1 (MQ=255)
                                                                         ccGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGa                                                                                                        >  1:873689/1‑116 (MQ=255)
                                                                         ccGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGa                                                                                                        <  2:873689/116‑1 (MQ=255)
                                                                          cGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATTAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCgg                                                                     >  1:1655141/1‑150 (MQ=255)
                                                                          cGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCgg                                                                     <  1:1175870/150‑1 (MQ=255)
                                                                           ggggTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGAt                                                                                                                        >  1:1463730/1‑98 (MQ=255)
                                                                           ggggTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGAt                                                                                                                        <  2:1463730/98‑1 (MQ=255)
                                                                             ggTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTcagccagc                                                                       >  2:493424/1‑145 (MQ=255)
                                                                             ggTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTcagccagc                                                                       <  1:493424/145‑1 (MQ=255)
                                                                               gaTTGCGCCTACAAGACCCAGGTGGTTGGGCGGTGTGACTTGCGCCGGGGCATAACTAGAGCGGCCGTTGGTACAGACTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAg                                                                <  1:1571968/149‑1 (MQ=255)
                                                                                   gcgcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCCTGAGGCGTGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTCAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTcc                                                            <  1:1894680/150‑1 (MQ=255)
                                                                                   gcgcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTcc                                                            <  1:192297/150‑1 (MQ=255)
                                                                                    cgcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTg                                                                                                               >  2:1425126/1‑98 (MQ=255)
                                                                                    cgcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTg                                                                                                               <  1:1425126/98‑1 (MQ=255)
                                                                                     gcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCgg                                                          >  2:495765/1‑150 (MQ=255)
                                                                                       cTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCa                                                        >  1:1439592/1‑150 (MQ=255)
                                                                                       cTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCa                                                        >  1:1641043/1‑150 (MQ=255)
                                                                                        tACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCaa                                                       <  2:1821044/150‑1 (MQ=255)
                                                                                        tACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCaa                                                       <  1:343595/150‑1 (MQ=255)
                                                                                           aaGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCa                                                    <  2:1173456/150‑1 (MQ=255)
                                                                                             gCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgggcg                                                                                                                                             >  2:1979982/1‑60 (MQ=255)
                                                                                             gCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgggcg                                                                                                                                             <  1:1979982/60‑1 (MQ=255)
                                                                                               ccAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTAAGCCAGCGGGGAAGCTCCGGCAAGCAcccc                                                >  1:1528362/1‑150 (MQ=255)
                                                                                                 aaGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGGGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCAcccccc                                              <  1:722664/150‑1 (MQ=255)
                                                                                                   gggggtTGGGCGTTTAGTCTTGCGCCGGTGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTg                                            <  2:923638/147‑1 (MQ=255)
                                                                                                     gggtTGGGCGTGGAGTCTTGCGCCGGGGCATAAATCGAGCGCGGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGAt                                          <  1:1556912/149‑1 (MQ=255)
                                                                                                       gtTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGa                                                                                                        <  2:584342/86‑1 (MQ=255)
                                                                                                       gtTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGa                                                                                                        >  1:584342/1‑86 (MQ=255)
                                                                                                        ttGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCg                                       >  1:1131245/1‑150 (MQ=255)
                                                                                                          gggCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGAtt                                                                                                                                  <  2:1341121/57‑1 (MQ=255)
                                                                                                          gggCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGAtt                                                                                                                                  >  1:1341121/1‑57 (MQ=255)
                                                                                                            gCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtgg                                   >  2:1798366/1‑150 (MQ=255)
                                                                                                              gTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCAcccccc                                              <  1:1038402/137‑1 (MQ=255)
                                                                                                              gTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCAcccccc                                              >  2:1038402/1‑137 (MQ=255)
                                                                                                               ttGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGCTGCCGACCTGGGCATGAAGCACTGAACCCACGGCCCCCCTTTCAGCCAGCGGCGAACCTCCGGCACGCCCCCCCCTCATCCGCtgctgg                                >  2:2233039/1‑150 (MQ=255)
                                                                                                                 gAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGAt                              >  2:920569/1‑150 (MQ=255)
                                                                                                                 gAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGAt                              <  2:37056/150‑1 (MQ=255)
                                                                                                                  aGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTg                                                                                                               >  1:1221408/1‑68 (MQ=255)
                                                                                                                  aGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTg                                                                                                               <  2:1221408/68‑1 (MQ=255)
                                                                                                                      ttGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCa                         <  1:806383/150‑1 (MQ=255)
                                                                                                                       tGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAg                        <  2:540147/150‑1 (MQ=255)
                                                                                                                            gggggCATAAGGCGAGCGCCGCGCGGGCGGCGCTCGCTTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCGCCCCCCTGATCCGCTGGTGGATGCTCAGGATaa                   <  2:1742201/149‑1 (MQ=255)
                                                                                                                            cGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATaa                   >  1:1169959/1‑150 (MQ=255)
                                                                                                                            cGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATaa                   >  1:12556/1‑150 (MQ=255)
                                                                                                                               ggCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTc                >  1:71103/1‑150 (MQ=255)
                                                                                                                                cacagaaTCGAGCGCCGCCCGCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtggt                                  <  2:1181077/126‑1 (MQ=255)
                                                                                                                                cacagaaTCGAGCGCCGCCCGCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtggt                                  >  1:1181077/6‑131 (MQ=255)
                                                                                                                                 cATAAAGCGAGCGACGCCCGCGCGCCGCAACATCCTGCGGCGATCCCGAGCTCGGCACGAAGCACTCAACCAACCGCTCCCATTTCAGCCAGCGGGGAAGCTACGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTg              <  2:714400/150‑1 (MQ=255)
                                                                                                                                 acagaaTCGAGCGCCGCCCGCGCGGCGCTCGATTTTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGc             <  1:1764983/146‑1 (MQ=255)
                                                                                                                                 acagaaTCGAGCGCCGCCCGCGCGGCGCTCGATTTTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGc             >  2:1764983/5‑150 (MQ=255)
                                                                                                                                  aTAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGc             <  1:920569/150‑1 (MQ=255)
                                                                                                                                   tAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCg            >  1:183474/1‑150 (MQ=255)
                                                                                                                                       tCGAGCGCCGCCCGCGCGGCGCTCGATTTTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCaa        >  2:1466955/1‑150 (MQ=255)
                                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGaa                                                                 <  1:1611960/91‑1 (MQ=255)
                                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGaa                                                                 >  2:1611960/1‑91 (MQ=255)
                                                                                                                                             gccgcGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTccc                                                                                   >  1:656204/1‑69 (MQ=255)
                                                                                                                                             gccgcGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTccc                                                                                   <  2:656204/69‑1 (MQ=255)
                                                                                                                                             gccgcGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCt  <  2:12556/150‑1 (MQ=255)
                                                                                                                                                  gcgggcggCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCg                                                           >  1:626764/1‑88 (MQ=255)
                                                                                                                                                  gcgggcggCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCg                                                           <  2:626764/88‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCT  >  NC_002947/5746395‑5746685

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