Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,426,684 G→C 5.3% H85D (CAC→GAC)  PP_2127 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,426,6840GC 5.3% 302.2 / 6.8 95H85D (CAC→GAC) PP_2127sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (45/45);  new base C (2/3);  total (47/48)
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.81e-01

GGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGACCGTTGC  >  NC_002947/2426538‑2426831
                                                                                                                                                    |                                                                                                                                                        
gggCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CCGGTGGAt                                                                                                                                                       >  2:807260/1‑150 (MQ=255)
    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                                                                                    >  1:821376/1‑135 (MQ=255)
                                                                             ggCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACca                                                                                    <  2:821376/135‑1 (MQ=255)
                                                                               cccAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGAtt                                                                   >  1:1080315/1‑150 (MQ=255)
                                                                                 cAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGAtt                                                                   <  1:812512/148‑1 (MQ=255)
                                                                                 cAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGAtt                                                                   >  2:812512/1‑148 (MQ=255)
                                                                                 cAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGc                                                                 >  2:2199400/1‑150 (MQ=255)
                                                                                    aTGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTg                                                              <  1:706231/150‑1 (MQ=255)
                                                                                    aTGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTg                                                              <  2:729578/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                                                                     <  1:1089937/134‑1 (MQ=255)
                                                                                             cTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGa                                                                     >  2:1089937/1‑134 (MQ=255)
                                                                                             cTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAg                                                     >  2:2284969/1‑150 (MQ=255)
                                                                                             cTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAg                                                     >  1:929206/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:1177204/1‑138 (MQ=255)
                                                                                                   ttcttGCACCTTGAGAATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                           >  2:1177149/1‑138 (MQ=255)
                                                                                                   ttcttGCACCTTGAGAATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                           <  1:1177204/138‑1 (MQ=255)
                                                                                                   ttcttGCACCTTGAGAATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGtt                                                           <  1:1177149/138‑1 (MQ=255)
                                                                                                    tcttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGc                                                                                                          <  1:1967526/95‑1 (MQ=255)
                                                                                                    tcttGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGc                                                                                                          >  2:1967526/1‑95 (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                                      <  1:2199400/150‑1 (MQ=255)
                                                                                                            cTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGtt                                      >  2:715772/1‑150 (MQ=255)
                                                                                                             ttGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTg                                     <  2:1530093/150‑1 (MQ=255)
                                                                                                             ttGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTg                                     <  1:2229379/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:1701773/1‑150 (MQ=255)
                                                                                                                    tGCGCCCGGCCGAGCT‑GGCCGAGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCa                              >  1:988274/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)
                                                                                                                                                    |                                                                                                                                                        
GGGCATGCGCCGTTGCAGCCACTGCAGTTGCAGCTTGGCAGCGGTGAGCAGTTGACCGATGTCGTCGTGCAGTTCCCGGCCCAGATGCTGGCGCTCGTCTTCTTGCACCTTGAGCATGCGCTCGGCC‑AGCTCGGCCGGGCG‑CAGGTGGATCGACCTGGCGCCCAGGCCCAGTTGTACAACAGCACCCAGCGTGGCGACC‑‑‑‑AACTGCA‑GCACCAGCAGGCTGGCCGGGATTGCCTGGTTGGCCAGGTACAGCGCCATGTTGGCAATCAGTGAGGCCACGCACAGCCCGACCGTTGC  >  NC_002947/2426538‑2426831

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: ‑