Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,426,689 G→C 5.3% S83W (TCG→TGG)  PP_2127 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,426,6890GC 5.3% 296.2 / 6.8 95S83W (TCG→TGG) PP_2127sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (43/46);  new base C (2/3);  total (46/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.74e-01

ATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGACCGTTGCCCA  >  NC_002947/2426542‑2426834
                                                                                                                                                     |                                                                                                                                                      
aTGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGAc                                                                                                                                                      <  2:352898/150‑1 (MQ=255)
 tGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGAcc                                                                                                                                                     >  2:2131234/1‑150 (MQ=255)
   cgcCGTGGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGAAGTTGACTGTTGGCGACGTGCGGTTCCCGTCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTg                                                                                                                                                   <  1:541455/150‑1 (MQ=255)
    gcCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTgg                                                                                                                                                  >  1:2295688/1‑150 (MQ=255)
     ccGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGgc                                                                                                                                                 >  1:1005918/1‑150 (MQ=255)
     ccGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGgc                                                                                                                                                 >  1:1673023/1‑150 (MQ=255)
       gTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGgcgc                                                                                                                                               >  1:869552/1‑150 (MQ=255)
        ttGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCgcc                                                                                                                                              >  1:729578/1‑150 (MQ=255)
            aGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCgcccagg                                                                                                                                          >  1:390051/1‑150 (MQ=255)
              ccACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCgcccaggcc                                                                                                                                        <  2:2188871/150‑1 (MQ=255)
                ccTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCgcccaggccca                                                                                                                                      <  2:1337295/149‑1 (MQ=255)
                aCTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCCAGGTCGATCCACCT‑GCgcccggc                                                                                                                                          >  2:2198149/1‑143 (MQ=255)
                aCTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCGCCAGGTCGATCCACCT‑GCgcccggc                                                                                                                                          <  1:2198149/146‑4 (MQ=255)
                 cTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGGTCGAc                                                                                                                                                      <  2:1782062/133‑1 (MQ=255)
                 cTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGGTCGAc                                                                                                                                                      >  1:1782062/1‑133 (MQ=255)
                  tGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGt                                                                                                                                    >  1:1103255/1‑150 (MQ=255)
                   gCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGtt                                                                                                                                   <  2:2307958/150‑1 (MQ=255)
                   gCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGtt                                                                                                                                   <  1:1963883/150‑1 (MQ=255)
                      gTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTa                                                                                                                                <  2:1568748/150‑1 (MQ=255)
                         gCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTacaa                                                                                                                             >  2:1311316/1‑150 (MQ=255)
                         gCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCAAGGCCCAGTTGTacaa                                                                                                                             >  2:1475216/1‑150 (MQ=255)
                                 gCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCAccc                                                                                                                     >  2:1765563/1‑150 (MQ=255)
                                    gCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGc                                                                                                                  <  1:2170686/150‑1 (MQ=255)
                                        ggAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTgg                                                                                                              <  2:390051/149‑1 (MQ=255)
                                            cAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGAc                                                                                                          >  2:2076103/1‑150 (MQ=255)
                                            cAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGACGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGAc                                                                                                          <  1:692071/150‑1 (MQ=255)
                                              gTTGACCGATGTCGTCGTGCAGTACCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑a                                                                                                    >  2:1184579/1‑150 (MQ=255)
                                                tGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AAc                                                                                                  <  2:1103255/150‑1 (MQ=255)
                                                    cGATGTCGTCGTGCAGTTCCCGGGCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTa                                                                                                                                >  2:1185015/1‑120 (MQ=255)
                                                    cGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTa                                                                                                                                <  1:1185015/120‑1 (MQ=255)
                                                    cGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTgca                                                                                              >  1:1048581/1‑150 (MQ=255)
                                                              gtGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACcagcag                                                                                   >  1:1814768/1‑150 (MQ=255)
                                                               tGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAgg                                                                                  >  1:1929864/1‑150 (MQ=255)
                                                               tGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAgg                                                                                  <  2:869552/150‑1 (MQ=255)
                                                                  aGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGc                                                                                                                  >  1:457302/1‑120 (MQ=255)
                                                                  aGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGc                                                                                                                  <  2:457302/120‑1 (MQ=255)
                                                                    ttCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGc                                                                             >  2:2431187/1‑150 (MQ=255)
                                                                       ccGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCgcccaggccca                                                                                                                                      <  1:821718/95‑1 (MQ=255)
                                                                       ccGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCgcccaggccca                                                                                                                                      >  2:821718/1‑95 (MQ=255)
                                                                         ggCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACca                                                                                       <  2:821376/135‑1 (MQ=255)
                                                                         ggCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACca                                                                                       >  1:821376/1‑135 (MQ=255)
                                                                           cccAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGAtt                                                                      >  1:1080315/1‑150 (MQ=255)
                                                                             cAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGAtt                                                                      >  2:812512/1‑148 (MQ=255)
                                                                             cAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGAtt                                                                      <  1:812512/148‑1 (MQ=255)
                                                                             cAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGc                                                                    >  2:2199400/1‑150 (MQ=255)
                                                                                aTGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTg                                                                 <  2:729578/150‑1 (MQ=255)
                                                                                aTGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTg                                                                 <  1:706231/150‑1 (MQ=255)
                                                                                 tGCTGGCGCTCGTCTTCTTGCACCTTGGGCATGCGCTCGGCC‑AGCGCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCGCAGGCCCAGTTGTACAACAGGACGCAGCGGGGCGACC‑‑‑‑AACTGCA‑GCGCCAGCAGGCTGGCCGGGATTGCCTgg                                                                <  1:2131234/150‑1 (MQ=255)
                                                                                   cTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                              <  2:1814768/150‑1 (MQ=255)
                                                                                   cTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                              <  1:1765563/150‑1 (MQ=255)
                                                                                   cTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                              <  2:507237/150‑1 (MQ=255)
                                                                                      gcgcTCGTCTTATTGCAACTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGc                                                           >  2:524123/1‑150 (MQ=255)
                                                                                       cgcTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGcc                                                          <  1:1010553/150‑1 (MQ=255)
                                                                                         cTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGa                                                                        >  2:1089937/1‑134 (MQ=255)
                                                                                         cTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGa                                                                        <  1:1089937/134‑1 (MQ=255)
                                                                                         cTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAg                                                        >  1:929206/1‑150 (MQ=255)
                                                                                         cTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAg                                                        >  2:2284969/1‑150 (MQ=255)
                                                                                             tcttctTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGTCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTAc                                                    >  2:929963/1‑150 (MQ=255)
                                                                                              cttcttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACa                                                   <  1:358009/150‑1 (MQ=255)
                                                                                               ttcttGCACCTTGAGAATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                              >  2:1177149/1‑138 (MQ=255)
                                                                                               ttcttGCACCTTGAGAATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                              <  1:1177149/138‑1 (MQ=255)
                                                                                               ttcttGCACCTTGAGAATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                              <  1:1177204/138‑1 (MQ=255)
                                                                                               ttcttGCACCTTGAGAATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                              >  2:1177204/1‑138 (MQ=255)
                                                                                                tcttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGc                                                                                                             >  2:1967526/1‑95 (MQ=255)
                                                                                                tcttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGc                                                                                                             <  1:1967526/95‑1 (MQ=255)
                                                                                                tcttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAgc                                                 <  2:1945808/150‑1 (MQ=255)
                                                                                                tcttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAgc                                                 <  1:2431187/150‑1 (MQ=255)
                                                                                                    gCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCa                                             <  1:111160/150‑1 (MQ=255)
                                                                                                       ccTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGt                                          >  1:206316/1‑150 (MQ=255)
                                                                                                        cTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGtt                                         >  2:715772/1‑150 (MQ=255)
                                                                                                        cTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGtt                                         <  1:2199400/150‑1 (MQ=255)
                                                                                                         ttGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTg                                        <  1:2229379/150‑1 (MQ=255)
                                                                                                         ttGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTg                                        <  2:1530093/150‑1 (MQ=255)
                                                                                                         ttGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACAGCC‑GGTTCTGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTg                                        >  2:64137/1‑150 (MQ=255)
                                                                                                         ttGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGACGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTg                                        <  2:713274/150‑1 (MQ=255)
                                                                                                          tGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTgg                                       <  1:417664/150‑1 (MQ=255)
                                                                                                                tGCGCCCGGCCGAGCT‑GGCCGAGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCa                                 >  1:988274/1‑150 (MQ=255)
                                                                                                                tGCGCCCGGCCGAGCT‑GGCCGAGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCa                                 >  1:1701773/1‑150 (MQ=255)
                                                                                                                    cTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGa                             <  2:2237689/150‑1 (MQ=255)
                                                                                                                     tCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAg                            <  2:263212/150‑1 (MQ=255)
                                                                                                                       ggCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGc                          <  2:1427474/150‑1 (MQ=255)
                                                                                                                          c‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCAc                       >  2:456383/1‑150 (MQ=255)
                                                                                                                          c‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCAc                       <  2:761716/150‑1 (MQ=255)
                                                                                                                            aGCTCGGCCGGGCG‑CAGGTGGATCAACCCGGCAGGCACACCCAGTTTTACAATAGCACCCCCCGTGGCGACCCCGTAATCGTACGCACCAAAACGC‑‑G‑CGGG‑‑TGCCTGGTTGGCCAGGTACAGCGCCAAGTTGGCAAGTTGTGAGGCCACg                      >  1:544567/1‑150 (MQ=255)
                                                                                                                              gctggCCGAGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGca                    <  2:988274/147‑1 (MQ=255)
                                                                                                                               tCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGcac                   >  1:608305/1‑150 (MQ=255)
                                                                                                                                 gcctgggCGCCAGGTCGATCCACCT‑GCGCCC‑GGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAgg                                                                                  >  2:550352/5‑84 (MQ=18)
                                                                                                                                 gcctgggCGCCAGGTCGATCCACCT‑GCGCCC‑GGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAgg                                                                                  <  1:550352/80‑1 (MQ=18)
                                                                                                                                      ggcgccAGGTCGATCCACCT‑GCGCCC‑GGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGAc           <  1:1838665/145‑1 (MQ=255)
                                                                                                                                      ggCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGa            >  2:498559/1‑150 (MQ=255)
                                                                                                                                            aGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGaataac       <  2:1005918/150‑6 (MQ=255)
                                                                                                                                             ggTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGTCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGACCGTTg      <  1:929963/150‑1 (MQ=255)
                                                                                                                                              gTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCTAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGACCGTTGc     >  2:1303449/1‑150 (MQ=255)
                                                                                                                                              gTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGACCGTTGc     >  2:1572462/1‑150 (MQ=255)
                                                                                                                                                 gATCGACCTGGCGCTAATGTCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AGCGACA‑TCACCACGACTCGGGCCCGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGCCCACGCACAGCCCGACCGTTGCCCa  <  1:1475216/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGACCGTTGCCCA  >  NC_002947/2426542‑2426834

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