Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,919,686 C→G 6.3% P91R (CCC→CGC) ‡ 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,6860CG 6.3% 240.4 / 6.7 79P91R (CCC→CGC) ‡PP_5190type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (35/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.43e-01

ACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTC  >  NC_002947/5919542‑5919833
                                                                                                                                                |                                                                                                                                                   
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                                                                                                      <  4:362921/127‑1 (MQ=255)
                                                                 ttGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCa                                                                                                      >  3:362921/1‑127 (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                                                                                                                                            <  3:19881/81‑3 (MQ=255)
                                                                         ccTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcga                                                                                                                                            >  4:19881/1‑79 (MQ=255)
                                                                                cAGTGGCTGGGAGAGATTGTCGTCCAACCCTACCTGCGCATGGACCCCTTGCAAATCTACCTGCCCCAGGTCACCAGGCTGATCTCCCCATCCTTCGCCCAGCGCCACGGCATCCTCTTCGTGGCGGCCGATGCCTCCAGCGTCACCGt                                                                 <  2:272592/149‑1 (MQ=255)
                                                                                cAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGGGCATCGACCCCATGCAACTCGACCTGCCCCAGGTAACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGACACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGt                                                                 >  2:110751/1‑149 (MQ=255)
                                                                                 aGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAg                                                                          <  5:38220/139‑1 (MQ=255)
                                                                                 aGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAg                                                                          >  6:38220/1‑139 (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                                                                <  5:33116/111‑1 (MQ=255)
                                                                                                                       aTCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTc                                                                >  6:33116/1‑111 (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           >  5:295962/1‑149 (MQ=255)
                                                                                                                                      cTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCt           <  1:298571/149‑1 (MQ=255)
                                                                                                                                         gACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCg         <  3:372088/148‑1 (MQ=255)
                                                                                                                                            cTGCCCGAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCg     <  5:71118/149‑1 (MQ=255)
                                                                                                                                               ccccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCACCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTc  <  1:50585/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTC  >  NC_002947/5919542‑5919833

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: