Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 35,012 C→T 8.4% intergenic (+30/+28) czrSA → / ← PP_0031 metal‑sensor histidine kinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294735,0120CT 8.4% 263.1 / 11.8 96intergenic (+30/+28)czrSA/PP_0031metal‑sensor histidine kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (45/42);  new base T (4/4);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

GATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTCTCTT  >  NC_002947/34865‑35161
                                                                                                                                                   |                                                                                                                                                     
gATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGa                                                                                                                                                     >  2:1054282/1‑150 (MQ=255)
gATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGa                                                                                                                                                     >  1:1936203/1‑150 (MQ=255)
   cATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                                  >  1:948330/1‑150 (MQ=255)
         tCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCttt                                                                                                                                            >  2:1203267/1‑150 (MQ=255)
             cACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttgtt                                                                                                                                        <  1:1005789/150‑1 (MQ=255)
              aCCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                                  >  2:1901051/1‑139 (MQ=255)
              aCCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                                  <  1:1901051/139‑1 (MQ=255)
              aCCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTtttgttt                                                                                                                                       <  1:2141301/150‑1 (MQ=255)
                gaGGGATCGGCGAGCGAGGAGGGGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCACCAACTGCAGGCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAAGCCCGACAATAGACGGGCTATTTGTTTGc                                                                                                                                     >  2:1338028/2‑150 (MQ=255)
                cAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGc                                                                                                                                     <  2:170029/150‑1 (MQ=255)
                 aGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGCCtaatcgt                                                                                                                                        >  1:959619/1‑140 (MQ=255)
                 aGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGCCtaatcgt                                                                                                                                        <  2:959619/146‑7 (MQ=255)
                 aGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGcc                                                                                                                                    <  2:79410/150‑1 (MQ=255)
                        ggCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGc                                                                                                                             >  1:994532/1‑150 (MQ=255)
                           gAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATc                                                                                                                          >  2:482208/1‑150 (MQ=255)
                                tGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGc                                                                                                                     <  1:1227861/150‑1 (MQ=255)
                                  gAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCTCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGc                                                                                                                             <  1:2154730/140‑1 (MQ=255)
                                  gAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCTCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGc                                                                                                                             >  2:2154730/1‑140 (MQ=255)
                                     gTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAg                                                                                                                <  1:1420507/150‑1 (MQ=255)
                                          cGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGGCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACca                                                                                                           >  1:1831147/1‑150 (MQ=255)
                                            aagaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGt                                                                                                                       >  2:326931/1‑136 (MQ=255)
                                            aagaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACcaca                                                                                                         <  1:1203267/150‑1 (MQ=255)
                                            aagaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACcaca                                                                                                         <  1:410761/150‑1 (MQ=255)
                                            aagaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGCGCCTTTTGTTTGCCTCAATGCATCAGt                                                                                                                       <  1:326931/136‑1 (MQ=255)
                                              gaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAgg                                                                                                       >  1:1465346/1‑150 (MQ=255)
                                              gaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAgg                                                                                                       <  2:1219072/150‑1 (MQ=255)
                                              gaaAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAgg                                                                                                       >  1:1465368/1‑150 (MQ=255)
                                               aaaaCTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGAGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAggg                                                                                                      >  2:457676/1‑150 (MQ=255)
                                                aaaCTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGt                                                                                                     <  1:1618159/150‑1 (MQ=255)
                                                   cTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc                                                                                                                                                >  2:2049180/1‑104 (MQ=255)
                                                   cTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCGATTGACGGGc                                                                                                                                                <  1:2049180/104‑1 (MQ=255)
                                                   cTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGaa                                                                                                                 <  1:431700/135‑1 (MQ=255)
                                                   cTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGaa                                                                                                                 >  2:431700/1‑135 (MQ=255)
                                                    ttCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAAt                                                                                                                               <  1:1113650/120‑1 (MQ=255)
                                                    ttCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAAt                                                                                                                               >  2:1113650/1‑120 (MQ=255)
                                                    ttCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCAt                                                                                                 <  1:1186614/150‑1 (MQ=255)
                                                        ccTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGc                                                                                             <  2:1021577/150‑1 (MQ=255)
                                                         cTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCg                                                                                            >  1:1209965/1‑150 (MQ=255)
                                                                cACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAAt                                                                                                                               <  2:1916185/108‑1 (MQ=255)
                                                                cACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAAt                                                                                                                               >  1:1916185/1‑108 (MQ=255)
                                                                cACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAAAGACCGACGATAAAATTGCGCAGACGCGATTAAGCCCATCAATCTACGGTCCTTGTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAg                                                                                     >  2:883064/1‑150 (MQ=255)
                                                                  ccTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCaa                                                                                                                                >  2:1367168/1‑105 (MQ=255)
                                                                  ccTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCaa                                                                                                                                <  1:1367168/105‑1 (MQ=255)
                                                                   cTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAAt                                                                                  >  2:2306227/1‑150 (MQ=255)
                                                                      cACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCAc                                                                               <  2:948330/150‑1 (MQ=255)
                                                                              ggAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                                  >  1:784116/1‑75 (MQ=255)
                                                                              ggAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACgg                                                                                                                                                  <  2:784116/75‑1 (MQ=255)
                                                                              ggAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCCCCGGGTCCATGGGCGACCGCAGAATCACCGACTtaa                                                                       >  1:243668/1‑150 (MQ=255)
                                                                                 aGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCttt                                                                                                                                            >  2:1713073/1‑78 (MQ=255)
                                                                                 aGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCttt                                                                                                                                            <  1:1713073/78‑1 (MQ=255)
                                                                                    ccAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtct                                                                 >  2:1381251/1‑150 (MQ=255)
                                                                                     cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCAt                                                                                                                           <  1:959085/91‑1 (MQ=255)
                                                                                     cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCAt                                                                                                                           >  2:959085/1‑91 (MQ=255)
                                                                                     cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtctt                                                                <  2:1936203/150‑1 (MQ=255)
                                                                                     cAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCATCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtctt                                                                >  2:1556800/1‑150 (MQ=255)
                                                                                                  cGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCa                                                                                      <  1:740428/115‑1 (MQ=255)
                                                                                                  cGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCa                                                                                      >  2:740428/1‑115 (MQ=255)
                                                                                                    gCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGAtt                                                 >  1:399853/1‑150 (MQ=255)
                                                                                                     cAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTg                                                <  2:1465346/150‑1 (MQ=255)
                                                                                                     cAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTg                                                <  2:1465368/150‑1 (MQ=255)
                                                                                                     cAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTg                                                <  1:271202/150‑1 (MQ=255)
                                                                                                       aTGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGc                                              <  2:1209965/150‑1 (MQ=255)
                                                                                                           ccgccgCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGttt                                          >  1:91736/1‑150 (MQ=255)
                                                                                                            cgccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCGTGAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGACAGGATTGGCGtttt                                         >  2:857742/1‑150 (MQ=255)
                                                                                                             gccgcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGttttt                                        >  1:1094726/1‑150 (MQ=255)
                                                                                                                gcTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTGTTTCca                                     >  1:714155/1‑150 (MQ=255)
                                                                                                                  ttAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCcaca                                   <  1:1054282/150‑1 (MQ=255)
                                                                                                                  ttAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCcaca                                   >  1:2287714/1‑150 (MQ=255)
                                                                                                                  tgaaaTTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCcaca                                   <  1:2306227/148‑1 (MQ=255)
                                                                                                                       ttGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGgtgt                              >  2:2292801/1‑150 (MQ=255)
                                                                                                                          cgcAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGTACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGCCAGGATTGGCGGTTTTCCACCGGTGTgag                           >  2:225152/1‑150 (MQ=255)
                                                                                                                          cgcAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCCCAGGTGTgag                           <  2:91736/150‑1 (MQ=255)
                                                                                                                           gcAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCcaca                                   >  1:1568068/1‑141 (MQ=255)
                                                                                                                           gcAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCcaca                                   <  2:1568068/141‑1 (MQ=255)
                                                                                                                           gcAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTgaga                          >  2:511021/1‑150 (MQ=255)
                                                                                                                            cAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGaa                         <  1:261942/150‑1 (MQ=255)
                                                                                                                               acacGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCAcc                                                                              >  1:840453/1‑94 (MQ=255)
                                                                                                                               acacGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCAcc                                                                              <  2:840453/94‑1 (MQ=255)
                                                                                                                                   gATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGa                  <  1:2275874/150‑1 (MQ=255)
                                                                                                                                   gATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGa                  <  2:714155/150‑1 (MQ=255)
                                                                                                                                    aaaaGGCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGAcc                                                                                         >  1:1437218/4‑78 (MQ=38)
                                                                                                                                    aaaaGGCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGAcc                                                                                         <  2:1437218/75‑1 (MQ=38)
                                                                                                                                      tAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCa               <  1:225152/150‑1 (MQ=255)
                                                                                                                                        ggCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGc                                              >  1:494656/1‑117 (MQ=255)
                                                                                                                                        ggCCCGTCGATTGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGc                                              <  2:494656/117‑1 (MQ=255)
                                                                                                                                         gCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATgg                                                                                               <  2:210746/67‑1 (MQ=255)
                                                                                                                                         gCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATgg                                                                                               >  1:210746/1‑67 (MQ=255)
                                                                                                                                             gTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtctt                                                                <  2:1633478/94‑1 (MQ=255)
                                                                                                                                             gTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAAtctt                                                                >  1:1633478/1‑94 (MQ=255)
                                                                                                                                             gTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACtt        >  1:423044/1‑150 (MQ=255)
                                                                                                                                                aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTAc                     <  1:187815/134‑1 (MQ=255)
                                                                                                                                                aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTAc                     >  2:187815/1‑134 (MQ=255)
                                                                                                                                                aaTCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTtct     >  2:227458/1‑150 (MQ=255)
                                                                                                                                                  tCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTctct   <  1:227458/150‑1 (MQ=255)
                                                                                                                                                  tCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTctct   >  2:281146/1‑150 (MQ=255)
                                                                                                                                                  tCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTctct   >  1:2020877/1‑150 (MQ=255)
                                                                                                                                                   cGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTCTCtt  <  2:399853/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATCATGAGTCTGCACCAGGGATCGGCGAGCGTGGAGGTGACCGAAGAAAACTTCACCTGGTTCCACCTTCACTTTCCGGAAGCCCAGCAAATGCAGCCGGCAATGACCGCCGCTTAAATTGCGCAGACACGATTAGGCCCGTCAATCGACGGGCCTTTTGTTTGCCTCAATGCATCAGTGCCGAAGGACCACAGGGTCCATGGGCGACCGCAGAATCACCGACTTAATAATCTTCTGTTCGTCAGGATTGGCGTTTTTCCACAGGTGTGAGAAGTACTGAGCATTCACTTTCTCTT  >  NC_002947/34865‑35161

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: