Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,919,685 C→G 6.3% P91A (CCC→GCC) ‡ 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,6850CG 6.3% 260.0 / 6.4 79P91A (CCC→GCC) ‡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 = 9.93e-01

GCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTC  >  NC_002947/5919539‑5919833
                                                                                                                                                  |                                                                                                                                                    
gCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGc                                                                                                                                                      <  7:246357/147‑1 (MQ=255)
   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                                                                                                                                           <  7:17576/121‑4 (MQ=255)
                                     aaGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcgag                                                                                                                                           >  8:17576/1‑118 (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                                                                                                                                            >  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                                                                          <  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                                                                >  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           >  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)
                                                                                                                                                  |                                                                                                                                                    
GCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTC  >  NC_002947/5919539‑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: