Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,919,687 C→G 6.3% P91P (CCC→CCG) ‡ PP_5190 → type II secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,919,6870CG 6.3% 249.5 / 6.8 79P91P (CCC→CCG) ‡PP_5190type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base G (2/3);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTCG  >  NC_002947/5919542‑5919834
                                                                                                                                                 |                                                                                                                                                   
aCCCGCTGGAGCACATCGCCGCCTGCGGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAg                                                                                                                                                  <  3:3572/149‑1 (MQ=255)
 cccGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAgg                                                                                                                                                 <  1:17560/149‑1 (MQ=255)
 cccGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAgg                                                                                                                                                 <  2:75622/149‑1 (MQ=255)
        gAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcgag                                                                                                                                            <  6:2505/147‑4 (MQ=255)
        gAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCAcc                                                                                                                                            <  5:23130/147‑1 (MQ=255)
         aGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGc                                                                                                                                         <  4:45475/149‑1 (MQ=255)
          gCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGcc                                                                                                                                        <  7:341718/149‑1 (MQ=255)
             caTCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTg                                                                                                                                      >  4:7201/1‑148 (MQ=255)
                   cgccTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTccc                                                                                                                               <  1:234204/149‑1 (MQ=255)
                    gccTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTcccc                                                                                                                              <  2:38417/149‑1 (MQ=255)
                           gCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTc                                                                                                                       <  3:7201/149‑1 (MQ=255)
                            ccTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCg                                                                                                                      <  2:332054/149‑1 (MQ=255)
                              tGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCCTCGACCCCAGGCAACTCGACCTGCCCCAGGGCACCAGCCTGATCTCCACAGCAGTCGcc                                                                                                                    >  8:266730/1‑149 (MQ=255)
                               ggAAACCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGGCCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGccc                                                                                                                   <  4:288483/149‑1 (MQ=255)
                                  aaGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcgag                                                                                                                                            >  8:17576/1‑118 (MQ=255)
                                  aaGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcgag                                                                                                                                            <  7:17576/121‑4 (MQ=255)
                                      cgccAGCACCCCGGCCAGCCACTGGACTTGGAGGACCTTTGCAAGTGGCTGGCAGACAAGGTCGGCCAACCCTACCTGAGCATCGACCCCATGCAACTCGACATGCCCAAGGTCACCAGCCTGATCTCCCCCGCCTTCGCCCAGCGCCa                                                                                                            >  8:267481/1‑149 (MQ=255)
                                         cAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACgg                                                                                                         >  6:341283/1‑149 (MQ=255)
                                               cccGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGCCCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCt                                                                                                   >  3:277182/1‑149 (MQ=255)
                                                ccGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTc                                                                                                  <  6:67077/149‑1 (MQ=255)
                                                     cAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCAtcgt                                                                                             >  1:264649/1‑149 (MQ=255)
                                                      aGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCAtcgtc                                                                                            >  3:335758/1‑149 (MQ=255)
                                                         cccTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTcgcg                                                                                         >  5:369412/1‑149 (MQ=255)
                                                               aCTTGGAGCACCTTTCCCAGTCGCTGTCCGACAAGGTCCGCCAACCCTACCTGCGCATCGACCCCAAGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGAt                                                                                   <  8:315981/149‑1 (MQ=255)
                                                                 ttGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCa                                                                                                       >  3:362921/1‑127 (MQ=255)
                                                                 ttGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCa                                                                                                       <  4:362921/127‑1 (MQ=255)
                                                                   ggAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCCGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCt                                                                               >  6:171465/1‑149 (MQ=255)
                                                                   ggAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCAAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGcc                                                                                <  2:246219/148‑1 (MQ=255)
                                                                     aGCACCATTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCAACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCAGATCTCCACACCCTACGCCAATCGCCACGGCATCCTCATCGTCGCGGCCGATcccac                                                                              >  6:103597/1‑146 (MQ=255)
                                                                       cACCTTTGCCAGTGGCTTGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCGCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCa                                                                            <  6:82226/148‑1 (MQ=255)
                                                                       cACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATcc                                                                                                    >  2:325526/1‑124 (MQ=255)
                                                                       cACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATcc                                                                                                    <  1:325526/124‑1 (MQ=255)
                                                                        aCCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGc                                                                          <  5:313033/149‑1 (MQ=255)
                                                                         ccTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcga                                                                                                                                             >  4:19881/1‑79 (MQ=255)
                                                                         ccTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcga                                                                                                                                             <  3:19881/81‑3 (MQ=255)
                                                                                cAGTGGCTGGGAGAGATTGTCGTCCAACCCTACCTGCGCATGGACCCCTTGCAAATCTACCTGCCCCAGGTCACCAGGCTGATCTCCCCATCCTTCGCCCAGCGCCACGGCATCCTCTTCGTGGCGGCCGATGCCTCCAGCGTCACCGt                                                                  <  2:272592/149‑1 (MQ=255)
                                                                                cAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGGGCATCGACCCCATGCAACTCGACCTGCCCCAGGTAACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGACACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGt                                                                  >  2:110751/1‑149 (MQ=255)
                                                                                 aGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAg                                                                           >  6:38220/1‑139 (MQ=255)
                                                                                 aGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAg                                                                           <  5:38220/139‑1 (MQ=255)
                                                                                 aGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTc                                                                 <  5:341283/149‑1 (MQ=255)
                                                                                  gTGGCTGGCCCACAAGGTCGGCCAACCCTACCTGCGAATAGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGAACTCCCCCCACTTCGCCCAGCGCCACGGCATCCTCAACGTCGCCGCCGATGCCGCCCGCGACACCggac                                                                >  6:253061/1‑146 (MQ=255)
                                                                                   tggctggcCGACAAGGTCGTCCAACCCTACCTGCGCATCGACCCCATTCAACTCGACCAGCCCCAGGTCACCAGCCTGATCTCCCCAGCCCTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGc                                                               <  6:215175/149‑1 (MQ=255)
                                                                                      ttggcCGACCTTGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCCCCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAg                                                            >  2:19265/2‑149 (MQ=255)
                                                                                       tggcCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAgc                                                           >  4:372088/1‑149 (MQ=255)
                                                                                                  gTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCACCAGGTCACCAGCCTGATCTCCCCAGTCTTCGCCCAGCGCCCCGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCttt                                                >  7:156566/1‑149 (MQ=255)
                                                                                                   tCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTc                                               >  5:17441/1‑149 (MQ=255)
                                                                                                      gCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAg                                            <  1:219892/149‑1 (MQ=255)
                                                                                                       ccaaggcTACCTGCTGATCGACCCCATGCATCTCGACCTGCCCCAGGTCACCCGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAgg                                           <  5:103597/143‑1 (MQ=255)
                                                                                                       ccAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAgg                                           <  2:264649/149‑1 (MQ=255)
                                                                                                       ccAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAg                                            >  3:22908/1‑148 (MQ=255)
                                                                                                       ccAACCCATCCTGCGCATCGACCCCAAGCAACACGTCCTGCCCCCGGTCACCAGCCTGATCTCCCCAGGCTTCGCCCTGCTCCACGCCATCCTCATCGTCTCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCa                                             <  8:156566/147‑1 (MQ=255)
                                                                                                          aCCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATg                                        <  7:271471/149‑1 (MQ=255)
                                                                                                           cccTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGa                                       >  3:164056/1‑149 (MQ=255)
                                                                                                             cTACCTGCGCATCGCCCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGCTGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTACCAGGATCACt                                     >  8:98292/1‑149 (MQ=255)
                                                                                                             cTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACCGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGACACCGTCGCCACCGCCCAGCCTTTCCAGGATCACt                                     >  6:381635/1‑149 (MQ=255)
                                                                                                              tACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTg                                    >  4:83025/1‑149 (MQ=255)
                                                                                                                ccTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTcc                                                                                                                                <  6:120748/55‑1 (MQ=255)
                                                                                                                ccTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTcc                                                                                                                                >  5:120748/1‑55 (MQ=255)
                                                                                                                ccTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTgg                                   <  5:381635/148‑1 (MQ=255)
                                                                                                                 cTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTgg                                   >  2:298571/1‑147 (MQ=255)
                                                                                                                   gcgcATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAgg                               >  2:247259/1‑149 (MQ=255)
                                                                                                                     gcATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCg                                                                                                                      <  6:137494/60‑1 (MQ=255)
                                                                                                                     gcATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCg                                                                                                                      >  5:137494/1‑60 (MQ=255)
                                                                                                                      cATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCg                            <  1:110751/149‑1 (MQ=255)
                                                                                                                       aTCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTc                                                                 >  6:33116/1‑111 (MQ=255)
                                                                                                                       aTCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTc                                                                 <  5:33116/111‑1 (MQ=255)
                                                                                                                            ccccATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGa                           >  6:91743/1‑144 (MQ=255)
                                                                                                                            ccccATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGa                           <  5:91743/144‑1 (MQ=255)
                                                                                                                               cATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGccc                   >  1:251217/1‑149 (MQ=255)
                                                                                                                                aTGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCg                  >  6:278224/1‑149 (MQ=255)
                                                                                                                                 ggCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGc                 <  6:369412/148‑1 (MQ=255)
                                                                                                                                   cAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCa                <  4:22908/148‑1 (MQ=255)
                                                                                                                                    aaCTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCCCCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGc              >  3:116120/1‑149 (MQ=255)
                                                                                                                                    aaCTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGc              >  8:327418/1‑149 (MQ=255)
                                                                                                                                     aCTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCa                                 >  8:178582/1‑129 (MQ=255)
                                                                                                                                      cTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCt            <  1:298571/149‑1 (MQ=255)
                                                                                                                                      cTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCt            >  5:295962/1‑149 (MQ=255)
                                                                                                                                         gACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCg          <  3:372088/148‑1 (MQ=255)
                                                                                                                                            cTGCCCGAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCg      <  5:71118/149‑1 (MQ=255)
                                                                                                                                               ccccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCACCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTc   <  1:50585/149‑1 (MQ=255)
                                                                                                                                                 ccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGt    <  2:375351/146‑1 (MQ=255)
                                                                                                                                                 ccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGt    >  1:375351/1‑146 (MQ=255)
                                                                                                                                                 ccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTCg  >  2:328472/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTCG  >  NC_002947/5919542‑5919834

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