Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,978,158 T→A 8.5% D9V (GAC→GTC)  PP_5240 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,978,1580TA 8.5% 273.7 / 14.5 96D9V (GAC→GTC) PP_5240two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base T (43/43);  new base A (3/5);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

Only 100 of 102 total aligned reads displayed.

CCCGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGCTGG  >  NC_002947/5978011‑5978307
                                                                                                                                                   |                                                                                                                                                     
cccGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCAtcgtcgacg                                                                                                                                                     >  2:2351962/1‑150 (MQ=255)
 ccGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTcga                                                                                                                                                    >  2:1456743/1‑150 (MQ=255)
  cGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTcgac                                                                                                                                                   >  2:2260514/1‑150 (MQ=255)
  cGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTcgac                                                                                                                                                   >  1:681315/1‑150 (MQ=255)
      aTTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCCGGTTGGCCCGGCTTTTGCGCTCCGAGGCATCGTCGACGACGATg                                                                                                                                               <  1:1261201/150‑1 (MQ=255)
         ttCATGGCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGcac                                                                                                                                            >  2:603982/1‑150 (MQ=255)
         ttCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGcac                                                                                                                                            <  2:681315/150‑1 (MQ=255)
         ttCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGACGACGATGcct                                                                                                                                            <  1:2351962/150‑3 (MQ=255)
           cAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCacac                                                                                                                                          <  1:2118209/150‑1 (MQ=255)
               tCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCct                                                                                                                                      >  1:142233/1‑150 (MQ=255)
                 cAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCctct                                                                                                                                    >  2:377749/1‑150 (MQ=255)
                 cAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGAACGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCctct                                                                                                                                    >  1:1262053/1‑150 (MQ=255)
                       cacacaGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCAt                                                                                                                              >  1:2129033/1‑150 (MQ=255)
                        acacaGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGAt                                                                                                                                                >  2:1603417/1‑131 (MQ=255)
                        acacaGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGAt                                                                                                                                                <  1:1603417/131‑1 (MQ=255)
                        acacaGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATg                                                                                                                             >  1:1960973/1‑150 (MQ=255)
                         cacaGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATgc                                                                                                                            >  2:1255926/1‑150 (MQ=255)
                            aGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACTACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgcc                                                                                                                         >  1:1751018/1‑150 (MQ=255)
                            aGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgcc                                                                                                                         >  1:176054/1‑150 (MQ=255)
                            aGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgcc                                                                                                                         >  2:2108007/1‑150 (MQ=255)
                            aGCCCGCCTCGCCCGCCAGATGCGAAGCCACCAAGGCGTGCCCGGCGGCAAACGAAAAGGTCTCGACCCCCGCCGAGATCAGCAGGTTGGCCAGGCGTTTGCGCTCCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgcc                                                                                                                         <  1:1329687/150‑1 (MQ=255)
                             ccccGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGcgccg                                                                                                                        <  1:953637/149‑1 (MQ=255)
                               ccGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCgccgcc                                                                                                                      >  2:1777732/1‑150 (MQ=255)
                                cGCATCGCCCGCCAGTGGCGAAGCCAGGAATTCGTGCCCGGCGGAAAAGGACAAGGTCTCGAACCCCGCCGAGCGCAGCATGTTGGCCCGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCGGCTCATGCGCCGccc                                                                                                                     >  2:1565639/1‑150 (MQ=255)
                                 gCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCa                                                                                                                    >  2:171420/1‑150 (MQ=255)
                                  ccTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAg                                                                                                                   >  2:616919/1‑150 (MQ=255)
                                      gcccgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGCCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTc                                                                                                               <  1:1103291/150‑1 (MQ=255)
                                        ccgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGc                                                                                                             <  1:2393291/150‑1 (MQ=255)
                                         cgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCa                                                                                                            <  2:1470341/150‑1 (MQ=255)
                                         cgccAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCa                                                                                                            >  2:1311745/1‑150 (MQ=255)
                                            cAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGAc                                                                                                         >  2:221981/1‑150 (MQ=255)
                                            cAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGAc                                                                                                         >  2:1837898/1‑150 (MQ=255)
                                            cAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGAc                                                                                                         >  1:1174515/1‑150 (MQ=255)
                                             aGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGAcc                                                                                                        <  1:316640/150‑1 (MQ=255)
                                                         aGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCgagga                                                                                            >  2:2160557/1‑150 (MQ=255)
                                                         aGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCgagga                                                                                            <  2:1234808/150‑1 (MQ=255)
                                                         aGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCgagga                                                                                            >  1:966569/1‑150 (MQ=255)
                                                             aTGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGtt                                                                                        <  2:607781/150‑1 (MQ=255)
                                                               gCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTtgt                                                                                      <  2:142233/150‑1 (MQ=255)
                                                               gCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTtgt                                                                                      <  2:1395566/150‑1 (MQ=255)
                                                                  tGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACAccc                                                                                                                                        >  2:52589/1‑97 (MQ=255)
                                                                  tGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACAccc                                                                                                                                        <  1:52589/97‑1 (MQ=255)
                                                                    cccGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGCCGACTATGCACATTCCTCTGCGCATGCGCCGCCCAGGCTCGCAGACTGCGCTCCGAGAAGGATGTGCATg                                                                                 <  2:1732921/150‑1 (MQ=255)
                                                                     ccGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCacac                                                                                                                                          <  2:1988366/92‑1 (MQ=255)
                                                                     ccGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCacac                                                                                                                                          >  1:1988366/1‑92 (MQ=255)
                                                                     ccGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATgg                                                                                <  1:1118357/150‑1 (MQ=255)
                                                                         cggAAAACGACACGGTCTCGAACCCCGCCGCGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGc                                                                            <  1:2056944/150‑1 (MQ=255)
                                                                           gAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCaa                                                                          <  1:1255926/150‑1 (MQ=255)
                                                                            aaaaCGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAg                                                                         <  1:2306209/150‑1 (MQ=255)
                                                                            aaaaCGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAg                                                                         <  1:322803/150‑1 (MQ=255)
                                                                                   cAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCcagc                                                                  <  1:603982/150‑1 (MQ=255)
                                                                                    aaGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCcagca                                                                 >  1:2172392/1‑150 (MQ=255)
                                                                                      ggTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCagcagc                                                               <  1:974680/150‑1 (MQ=255)
                                                                                        tctcGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCagcga                                                             >  2:2004403/1‑150 (MQ=255)
                                                                                          tcGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTcc                                                                     >  2:33623/1‑140 (MQ=255)
                                                                                          tcGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTcc                                                                     <  1:33623/140‑1 (MQ=255)
                                                                                          tcGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTATTGCGCACCGAGGCATCGTCGACGACGATGc                                                                                                                                              <  2:318762/67‑1 (MQ=255)
                                                                                          tcGAACCCCGCCGAGCGCAGCAGGTTGGCCAAGCTTTTGCGCACCGAGGCATCGTCGACGACGATGc                                                                                                                                              >  1:318762/1‑67 (MQ=39)
                                                                                            gAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGGCGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCTGAGGAGGTTGTGCATGACCGCAAGATCCAGCAGCGAGCgg                                                         >  1:701271/1‑150 (MQ=255)
                                                                                            gAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCgg                                                         <  2:1686902/150‑1 (MQ=255)
                                                                                            gAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCgg                                                         >  2:250354/1‑150 (MQ=255)
                                                                                             aaCCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCgagga                                                                                            <  1:316923/114‑1 (MQ=255)
                                                                                             aaCCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCgagga                                                                                            >  2:316923/1‑114 (MQ=255)
                                                                                              aCCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGa                                                       <  1:2160557/150‑1 (MQ=255)
                                                                                                cccGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGGCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACt                                                     >  2:1411573/1‑150 (MQ=255)
                                                                                                  cgccgAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGc                                                   <  1:1837898/150‑1 (MQ=255)
                                                                                                    ccgAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCAt                                                 <  2:176054/150‑1 (MQ=255)
                                                                                                      gAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGtt                                                                                        >  2:1491088/1‑109 (MQ=255)
                                                                                                      gAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGtt                                                                                        <  1:1491088/109‑1 (MQ=255)
                                                                                                          gcagcaGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGtt                                           <  2:1174515/150‑1 (MQ=255)
                                                                                                          gcagcaGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGtt                                           <  2:1027010/150‑1 (MQ=255)
                                                                                                            agcagGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCa                                         >  2:1788694/1‑150 (MQ=255)
                                                                                                                 gTTGGCTCGGCTTTTGCGCACCAAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCatgatg                                    <  2:2206244/150‑1 (MQ=255)
                                                                                                                 gTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCatgatg                                    <  2:2129033/150‑1 (MQ=255)
                                                                                                                     gCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGc                                                   <  1:523127/131‑1 (MQ=255)
                                                                                                                     gCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGc                                                   >  2:523127/1‑131 (MQ=255)
                                                                                                                      ccAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttctt                               <  1:1456743/150‑1 (MQ=255)
                                                                                                                      ccAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttctt                               <  2:2172392/150‑1 (MQ=255)
                                                                                                                       cAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttcttc                              >  1:132364/1‑150 (MQ=255)
                                                                                                                       cAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGttcttc                              <  2:1548783/150‑1 (MQ=255)
                                                                                                                        aGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCCTGATGTtcttct                             >  2:305237/1‑150 (MQ=255)
                                                                                                                        aGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTtcttct                             >  1:1106468/1‑150 (MQ=255)
                                                                                                                         ggCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTcttctt                            <  1:669886/150‑1 (MQ=255)
                                                                                                                            ttttGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCGTGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCAt                                                                                  >  1:1870892/1‑93 (MQ=255)
                                                                                                                            ttttGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCAt                                                                                  <  2:1870892/93‑1 (MQ=255)
                                                                                                                                     ccGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGgt                <  1:124269/150‑1 (MQ=255)
                                                                                                                                      gtgtgCATCGTCGACGACGATGCCCACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATc                                                <  2:2027462/113‑1 (MQ=255)
                                                                                                                                      gtgtgCATCGTCGACGACGATGCCCACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATc                                                >  1:2027462/5‑117 (MQ=255)
                                                                                                                                      cGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGgtc               <  1:1311745/150‑1 (MQ=255)
                                                                                                                                      cGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGgtc               >  1:553985/1‑150 (MQ=255)
                                                                                                                                      cGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGgtc               <  1:377749/150‑1 (MQ=255)
                                                                                                                                           cATCGTCGACGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGTTCGTCAc          <  2:90888/150‑1 (MQ=255)
                                                                                                                                           cATCGTCGACGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCAc          >  1:90888/1‑150 (MQ=255)
                                                                                                                                            aTCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCAcc         >  2:699054/1‑150 (MQ=255)
                                                                                                                                            aTCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCAcc         >  2:2119611/1‑150 (MQ=255)
                                                                                                                                                tcgtcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGc     >  1:235228/1‑150 (MQ=255)
                                                                                                                                                tcgtcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGc     >  2:1025201/1‑150 (MQ=255)
                                                                                                                                                  gtcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGCTg   >  2:2280186/1‑150 (MQ=255)
                                                                                                                                                   tcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGCTgg  >  2:713803/1‑150 (MQ=255)
                                                                                                                                                   tcgACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATCATCCTCGTTTTGTGCACCTTGGTCGTCACCTCGCTgg  >  2:1397989/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCCGGCATTTTCAGGTCCAGCAACACACAGCCCGCCTCGCCCGCCAGTGGCGAAGCCAGGAATGCGTGCCCGGCGGAAAACGACAAGGTCTCGAACCCCGCCGAGCGCAGCAGGTTGGCCAGGCTTTTGCGCACCGAGGCATCGTCGTCGACGATGCACACCCCTCTGCTCATGCGCCGCCCAGGCTCGCAGACCGCGCTCCGAGGAGGTTGTGCATGGCCGCAAGCTCCAGCAGCGAGCGGGACTGCATCTTGTTCATGATGTTCTTCTTGTGCACCTTGGTCGTCACCTCGCTGG  >  NC_002947/5978011‑5978307

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