Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,941,121:1 +A 5.6% coding (123/531 nt) PP_3475 ← type II secretion system pseudopilin PulG‑like signal peptide protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,941,1211.A 5.6% 362.9 / 13.4 89coding (123/531 nt)PP_3475type II secretion system pseudopilin PulG‑like signal peptide protein
Reads supporting (aligned to +/- strand):  ref base . (40/44);  new base A (2/3);  total (42/47)
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.64e-01

GCCTGCGGGTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCT  >  NC_002947/3940975‑3941267
                                                                                                                                                   |                                                                                                                                                  
gCCTGCGGGTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGgcggcg                                                                                                                                                 >  1:53942/1‑150 (MQ=255)
 ccTGCGGGTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGgcgc                                                                                                                                                >  1:1111499/1‑150 (MQ=255)
  cTGCGGGTAGCGGTTGTCCTGGAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGgcgcc                                                                                                                                               >  2:760967/1‑149 (MQ=255)
        gTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCa                                                                                                                                         <  2:53942/150‑1 (MQ=255)
         tAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGCCGCCCAg                                                                                                                                        <  1:689561/150‑1 (MQ=255)
          aGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGa                                                                                                                                       <  1:786932/150‑1 (MQ=255)
                  tCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCg                                                                                                                               >  1:1816814/1‑150 (MQ=255)
                  tCCTGTAGCAGGTCGGCCAGATCCTGGGGCTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTCGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCTAGATCGTGccc                                                                                                                               >  2:656338/1‑149 (MQ=255)
                    cTGTAGCAGGTCGGCCAGATCCTGGGGGTAGTGGTCCAGCCCTGGCGCGGCGCGGTAGTAGCTGCGCAACGCTTTGGCATACTCGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCTCCCCGATCGTGCCGGt                                                                                                                             >  2:715988/1‑150 (MQ=255)
                    cTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGgcgc                                                                                                                                                <  2:716861/131‑1 (MQ=255)
                    cTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGgcgc                                                                                                                                                >  1:716861/1‑131 (MQ=255)
                    cTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGt                                                                                                                             >  2:939431/1‑150 (MQ=255)
                     tGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGCCGCCCAgc                                                                                                                                       >  1:505139/1‑138 (MQ=255)
                     tGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGCCGCCCAgc                                                                                                                                       <  2:505139/139‑2 (MQ=255)
                     tGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTg                                                                                                                            <  1:746515/150‑1 (MQ=255)
                      gTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTgg                                                                                                                           >  1:121334/1‑150 (MQ=255)
                        aGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGgcc                                                                                                                         <  1:49355/150‑1 (MQ=255)
                            ggTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCa                                                                                                                     <  2:670475/150‑1 (MQ=255)
                              tCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAgg                                                                                                                   <  2:2128109/150‑1 (MQ=255)
                              tCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGATAGGCCATCCTGGATGCACGCCCACCGCATCTGGCAAGACAGTTCGGCCGGGGCGGTGCTGACGCAGATCGTGCCGGTGGCCGCCAgg                                                                                                                   <  2:1480421/150‑1 (MQ=255)
                               cGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAggg                                                                                                                  <  1:939431/150‑1 (MQ=255)
                               cGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAggg                                                                                                                  <  2:1688970/150‑1 (MQ=255)
                                ggCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGc                                                                                                                 >  2:806810/1‑150 (MQ=255)
                                 gCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCg                                                                                                                               >  1:521409/1‑135 (MQ=255)
                                 gCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCg                                                                                                                               <  2:521409/135‑1 (MQ=255)
                                 gCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGGCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGcc                                                                                                                >  2:672251/1‑150 (MQ=255)
                                  ccAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGTATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCa                                                                                                               >  2:2039810/1‑150 (MQ=255)
                                   cAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCAc                                                                                                              >  1:2405017/1‑150 (MQ=255)
                                       tCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACgctg                                                                                                          >  2:1643149/1‑150 (MQ=255)
                                              ggTACTGGGCCCGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGCCGTTCGCGGTCGCGCTGGGCAGCGCTGGCCCACGTCGTGCCGGTGGCCGCCAGGGCCACGCTGCGGACCg                                                                                                   >  2:690696/1‑150 (MQ=255)
                                                  cTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGTATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAAt                                                                                               <  1:2039810/150‑1 (MQ=255)
                                                  cTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAAt                                                                                               <  2:959266/150‑1 (MQ=255)
                                                     ggCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGCCGCCCAgc                                                                                                                                       <  1:1570916/107‑2 (MQ=255)
                                                     ggCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGCCGCCCAgc                                                                                                                                       >  2:1570916/1‑106 (MQ=255)
                                                      gCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagc                                                                                           <  2:1816814/150‑1 (MQ=255)
                                                       ccAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                          <  2:282822/150‑1 (MQ=255)
                                                        cAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCCACCCACAGCAACTGGCTTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATTGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcag                                                                                         <  2:122576/150‑1 (MQ=255)
                                                              tGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                          >  2:875833/1‑143 (MQ=255)
                                                              tGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                          <  1:875833/143‑1 (MQ=255)
                                                                gCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCAccc                                                                                 >  2:2371532/1‑150 (MQ=255)
                                                                  gAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGt                                                                                                                                   >  2:2130591/1‑98 (MQ=255)
                                                                  gAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGt                                                                                                                                   <  1:2130591/98‑1 (MQ=255)
                                                                  gAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGgccgc                                                                                                                       <  2:2299312/110‑1 (MQ=255)
                                                                  gAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGgccgc                                                                                                                       >  1:2299312/1‑110 (MQ=255)
                                                                  gAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCAccccc                                                                               >  1:2403689/1‑150 (MQ=255)
                                                                   aGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCAcccccc                                                                              >  2:969459/1‑149 (MQ=255)
                                                                        cgGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTAt                                                                         <  1:1643149/150‑1 (MQ=255)
                                                                          gtagtaGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCATGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGt                                                                       >  1:1467868/1‑150 (MQ=255)
                                                                             gtagCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTaaa                                                                    <  1:2371532/150‑1 (MQ=255)
                                                                                gCTGCGCCAGGCTTGGGGCTGCTGGCGGGCGCCCCTCATCAACGTTCTTTCTCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAgcc                                                                 <  2:1467868/150‑1 (MQ=255)
                                                                                gCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                          <  2:1548722/125‑1 (MQ=255)
                                                                                gCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                          >  1:1548722/1‑125 (MQ=255)
                                                                                     gcAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGt                                                            >  2:513663/1‑150 (MQ=255)
                                                                                      cAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcag                                                                                         <  1:1525750/120‑1 (MQ=255)
                                                                                      cAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcag                                                                                         >  2:1525750/1‑120 (MQ=255)
                                                                                       aaGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGc                                                                                                                                 <  1:2140898/79‑1 (MQ=255)
                                                                                       aaGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGc                                                                                                                                 >  2:2140898/1‑79 (MQ=255)
                                                                                       aaGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTg                                                          <  1:513663/150‑1 (MQ=255)
                                                                                         ggCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGt                                                                                                                                   <  2:1805564/75‑1 (MQ=255)
                                                                                         ggCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGt                                                                                                                                   >  1:1805564/1‑75 (MQ=255)
                                                                                         ggCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGc                                                        <  1:969459/150‑1 (MQ=255)
                                                                                          gCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGcc                                                       >  1:1215967/1‑150 (MQ=255)
                                                                                             tGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTa                                                    <  2:1915102/150‑1 (MQ=255)
                                                                                              gggCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTAc                                                   >  1:2006320/1‑150 (MQ=255)
                                                                                                      tGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGc                                           >  1:2248692/1‑150 (MQ=255)
                                                                                                      tGGCTGCCGACCCACAGCAACTGGCGTGCGCGGTCGCGCGGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGc                                           <  1:715988/150‑1 (MQ=255)
                                                                                                          tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTAc                                                   <  2:2169467/138‑1 (MQ=255)
                                                                                                          tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTAc                                                   >  1:2169467/1‑138 (MQ=255)
                                                                                                          tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAag                                       <  2:2405017/150‑1 (MQ=255)
                                                                                                          tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAag                                       >  2:730800/1‑150 (MQ=255)
                                                                                                           gCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAaga                                      <  2:175918/150‑1 (MQ=255)
                                                                                                           gCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAaga                                      <  1:672251/150‑1 (MQ=255)
                                                                                                            ccGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAagag                                     <  2:1215967/150‑1 (MQ=255)
                                                                                                                 ccACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCt                                  >  2:2256637/1‑148 (MQ=255)
                                                                                                                 ccACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCt                                  <  1:2256637/148‑1 (MQ=255)
                                                                                                                 ccACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGc                                >  2:1536291/1‑150 (MQ=255)
                                                                                                                 ccACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCCCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGc                                >  2:2133054/1‑150 (MQ=255)
                                                                                                                  cacaGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCGAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCAGTCAGCCAAGATGCTGcc                               >  2:1721589/1‑150 (MQ=255)
                                                                                                                   acaGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCa                              >  1:98088/1‑150 (MQ=255)
                                                                                                                     aGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTAc                                                   >  1:1902961/1‑127 (MQ=255)
                                                                                                                     aGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTAc                                                   <  2:1902961/127‑1 (MQ=255)
                                                                                                                       cAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCt                                  >  1:467344/1‑142 (MQ=255)
                                                                                                                       cAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCAc                          <  2:2086439/150‑1 (MQ=255)
                                                                                                                       cAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCAc                          <  1:2133054/150‑1 (MQ=255)
                                                                                                                       cAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCAAGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCt                                  <  2:467344/142‑1 (MQ=255)
                                                                                                                       cAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGAGCCTGTGTCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCAc                          >  2:340357/1‑150 (MQ=255)
                                                                                                                           tGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGATTGCGGTcagcagc                                                                                        >  2:2160175/1‑84 (MQ=38)
                                                                                                                           tGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACgcgc                      >  1:968462/1‑150 (MQ=255)
                                                                                                                                 tCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGccctccc                <  2:1111499/150‑1 (MQ=255)
                                                                                                                                    cgGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCtgt             <  2:968462/150‑1 (MQ=255)
                                                                                                                                    cgGTCGCGCTGGGCGGCCCTCGCCCAGATTGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCtgt             <  2:2361331/150‑1 (MQ=255)
                                                                                                                                     ggTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCtgtg            <  1:1536291/150‑1 (MQ=255)
                                                                                                                                       tCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGc          >  2:2334013/1‑150 (MQ=255)
                                                                                                                                          cgcTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCAcc       <  2:121334/150‑1 (MQ=255)
                                                                                                                                               ggcggcGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGctct  <  1:806810/150‑1 (MQ=255)
                                                                                                                                               ggcggcGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGctct  >  1:278830/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCTGCGGGTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCT  >  NC_002947/3940975‑3941267

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