Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,941,126 Δ1 bp 5.4% coding (118/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,1260T. 5.4% 278.0 / 9.9 93coding (118/531 nt)PP_3475type II secretion system pseudopilin PulG‑like signal peptide protein
Reads supporting (aligned to +/- strand):  ref base T (42/46);  new base . (2/3);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

GTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCG  >  NC_002947/3940983‑3941275
                                                                                                                                                |                                                                                                                                                     
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)
            cTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGt                                                                                                                                     >  2:939431/1‑150 (MQ=255)
             tGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGCCGCCCAgc                                                                                                                                               <  2:505139/139‑2 (MQ=255)
             tGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGCCGCCCAgc                                                                                                                                               >  1:505139/1‑138 (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                                                                                                                                               >  2:1570916/1‑106 (MQ=255)
                                             ggCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGAGCGCCGCCCAgc                                                                                                                                               <  1:1570916/107‑2 (MQ=255)
                                              gCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagc                                                                                                   <  2:1816814/150‑1 (MQ=255)
                                               ccAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                                  <  2:282822/150‑1 (MQ=255)
                                                cAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCCACCCACAGCAACTGGCTTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATTGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcag                                                                                                 <  2:122576/150‑1 (MQ=255)
                                                      tGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                                  <  1:875833/143‑1 (MQ=255)
                                                      tGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                                  >  2:875833/1‑143 (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                                                                                                  >  1:1548722/1‑125 (MQ=255)
                                                                        gCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagca                                                                                                  <  2:1548722/125‑1 (MQ=255)
                                                                             gcAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGt                                                                    >  2:513663/1‑150 (MQ=255)
                                                                              cAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcag                                                                                                 <  1:1525750/120‑1 (MQ=255)
                                                                              cAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATcagcag                                                                                                 >  2:1525750/1‑120 (MQ=255)
                                                                               aaGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGc                                                                                                                                         >  2:2140898/1‑79 (MQ=255)
                                                                               aaGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGc                                                                                                                                         <  1:2140898/79‑1 (MQ=255)
                                                                               aaGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTg                                                                  <  1:513663/150‑1 (MQ=255)
                                                                                 ggCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGt                                                                                                                                           >  1:1805564/1‑75 (MQ=255)
                                                                                 ggCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGt                                                                                                                                           <  2:1805564/75‑1 (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                                                           >  1:2169467/1‑138 (MQ=255)
                                                                                                  tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTAc                                                           <  2:2169467/138‑1 (MQ=255)
                                                                                                  tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAag                                               >  2:730800/1‑150 (MQ=255)
                                                                                                  tGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAag                                               <  2:2405017/150‑1 (MQ=255)
                                                                                                   gCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAaga                                              <  2:175918/150‑1 (MQ=255)
                                                                                                   gCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAaga                                              <  1:672251/150‑1 (MQ=255)
                                                                                                    ccGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAagag                                             <  2:1215967/150‑1 (MQ=255)
                                                                                                         ccACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCt                                          <  1:2256637/148‑1 (MQ=255)
                                                                                                         ccACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCt                                          >  2:2256637/1‑148 (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                                  <  1:2133054/150‑1 (MQ=255)
                                                                                                               cAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCAc                                  <  2:2086439/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)
                                                                                                                                              gcTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGt    >  1:1806159/1‑150 (MQ=255)
                                                                                                                                               cTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTc   <  1:2334013/150‑1 (MQ=255)
                                                                                                                                                tCGCCCAGATGGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCg  >  1:258431/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GTAGCGGTTGTCCTGTAGCAGGTCGGCCAGATCCTGGGGGTACTGGGCCAGCCCTGGCGAGGCGCGGTAGTAGCTGCGCAAGGCTTGGGCATACTGGCTGCCGACCCACAGCAACTGGCGTTCGCGGTCGCGCTGGGCGGCGCTCGCCCAGATCGTGCCGGTGGCCGCCAGGGCCACGCTGCTGACCGCAATCAGCAGCAGCACCCCCAGGTATGTAAAGCCGCCGTTGGCCTTACCATTCAGCGAAGAGGCTGCCATCACGCGCCCTCCCTGTGGCACCGCTCTTGATGTCG  >  NC_002947/3940983‑3941275

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: