Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,884,682 C→T 7.8% L175L (CTG→TTG)  PP_3429 → two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,884,6820CT 7.8% 264.1 / 8.2 89L175L (CTG→TTG) PP_3429two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (40/42);  new base T (2/5);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGGCGGCGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGTT  >  NC_002947/3884541‑3884819
                                                                                                                                             |                                                                                                                                         
cggcggCGATAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgat                                                                                                                                   <  1:632732/150‑3 (MQ=255)
cggcggCGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgat                                                                                                                                   >  1:629028/1‑148 (MQ=255)
cggcggCGACAGTACGTGACACGTGAATGGCTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCcggc                                                                                                                                   >  1:1643922/1‑150 (MQ=255)
 ggcggcGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgatg                                                                                                                                  <  2:629028/150‑4 (MQ=255)
 ggcggcGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCcggcg                                                                                                                                  <  2:36365/150‑1 (MQ=255)
 ggcggcGACAGTACGGGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGAAAGAAGTGCAGGCCAGCTCTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgatg                                                                                                                                  <  1:903139/150‑4 (MQ=255)
   cggcgACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCggcggc                                                                                                                                <  2:1435041/150‑1 (MQ=255)
   cggcgACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCggcggc                                                                                                                                <  2:1643922/150‑1 (MQ=255)
    ggcgACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGgcggcg                                                                                                                               >  2:2191225/1‑150 (MQ=255)
       gACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCg                                                                                                                            >  1:2088054/1‑150 (MQ=255)
           gTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCa                                                                                                                        <  1:1249617/150‑1 (MQ=255)
              cGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGa                                                                                                                     <  1:2039018/150‑1 (MQ=255)
              cGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGa                                                                                                                     <  1:1247229/150‑1 (MQ=255)
                 gACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGAGCGTGACtt                                                                                                                  >  1:1307161/1‑150 (MQ=255)
                  acacGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTc                                                                                                                 <  2:2088054/150‑1 (MQ=255)
                    acGTGATCATCGCCATGGCAGGCAGGAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCCACATGGCGTTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTcaa                                                                                                               <  1:1955479/150‑1 (MQ=255)
                    acGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTcaa                                                                                                               <  2:482305/150‑1 (MQ=255)
                    acGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCATTCCACGCAGGAGAAAGTGCAGGACCGCGGGGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTcaa                                                                                                               >  1:1499864/1‑150 (MQ=255)
                         aTCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAAcc                                                                                                          >  1:1092193/1‑150 (MQ=255)
                            aaCTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgatg                                                                                                                                  >  1:1585858/1‑120 (MQ=255)
                            aaCTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAGTTGGCcgatg                                                                                                                                  <  2:1585858/123‑4 (MQ=255)
                             aCTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCAt                                                                                                                       <  2:692366/133‑1 (MQ=255)
                             aCTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCAt                                                                                                                       >  1:692366/1‑133 (MQ=255)
                               tCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCa                                                                                                    >  1:1193591/1‑150 (MQ=255)
                                    ggCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGc                                                                                               <  2:1619481/150‑1 (MQ=255)
                                     gCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCa                                                                                              >  2:734967/1‑150 (MQ=255)
                                       aTGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGCCCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGAc                                                                                                                    <  1:900540/126‑1 (MQ=255)
                                       aTGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGAc                                                                                                                    >  2:900540/1‑126 (MQ=255)
                                       aTGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATc                                                                                            >  1:52784/1‑150 (MQ=255)
                                        tGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCa                                                                                           <  1:2181040/150‑1 (MQ=255)
                                         gcagcaGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAg                                                                                          >  1:556172/1‑150 (MQ=255)
                                          cagcagCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGc                                                                                         <  2:2146008/150‑1 (MQ=255)
                                           agcagcGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGct                                                                                        <  1:1071337/150‑2 (MQ=255)
                                           agcagcGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCg                                                                                        >  1:1091765/1‑150 (MQ=255)
                                                   cTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCaacaac                                                                                                           <  2:939119/123‑1 (MQ=255)
                                                   cTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCaacaac                                                                                                           >  1:939119/1‑123 (MQ=255)
                                                     ggACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCGAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGa                                                                              <  1:58519/150‑1 (MQ=255)
                                                      gACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGaa                                                                             <  1:775330/150‑1 (MQ=255)
                                                         aGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCg                                                                                        >  2:2348415/1‑136 (MQ=255)
                                                         aGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCg                                                                                        <  1:2348415/136‑1 (MQ=255)
                                                            cGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATc                                                                       <  2:52784/150‑1 (MQ=255)
                                                              gCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCAc                                                                                                   <  2:930997/120‑1 (MQ=255)
                                                              gCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCAc                                                                                                   >  1:930997/1‑120 (MQ=255)
                                                              gCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCa                                                                     <  1:384850/150‑1 (MQ=255)
                                                               ccTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAg                                                                    >  2:1621708/1‑150 (MQ=255)
                                                                cTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAga                                                                   >  1:238003/1‑149 (MQ=255)
                                                                     aCGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGgcgc                                                              <  2:152758/150‑1 (MQ=255)
                                                                       gCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCt                                                            <  1:1223767/150‑1 (MQ=255)
                                                                          ggAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTgct                                                         >  1:533677/1‑149 (MQ=255)
                                                                            aagaagTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGcc                                                       >  2:2099489/1‑150 (MQ=255)
                                                                               aagTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAg                                                    <  2:2364914/150‑1 (MQ=255)
                                                                                     aGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCa                                              <  2:1091765/150‑1 (MQ=255)
                                                                                          cgcgGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGAc                                         <  2:262016/150‑1 (MQ=255)
                                                                                          cgcgGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGAc                                         >  1:1796409/1‑150 (MQ=255)
                                                                                          cgcgGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTAAACAGCATCAGCGTCTGCTTTGAATTCATCTATTTTCGCCTGCGCCAAGGCCGCAGAGAc                                         >  2:2202116/1‑150 (MQ=255)
                                                                                                  gTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGa                                 >  1:890489/1‑150 (MQ=255)
                                                                                                   ttCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAgcggc                                                                >  1:1930573/1‑118 (MQ=255)
                                                                                                   ttCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAgcggc                                                                <  2:1930573/118‑1 (MQ=255)
                                                                                                   ttCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGAt                                >  2:2167517/1‑150 (MQ=255)
                                                                                                   ttCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGAt                                >  1:1485929/1‑150 (MQ=255)
                                                                                                   ttCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTACGCCAAGGCCGCAGCGAAGGCCTGGAt                                >  1:500598/1‑150 (MQ=255)
                                                                                                       cccGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCg                            <  2:2114383/150‑1 (MQ=255)
                                                                                                        ccGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGt                           <  1:1526748/150‑1 (MQ=255)
                                                                                                         cGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGtt                          >  2:2328501/1‑150 (MQ=255)
                                                                                                            tGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCtg                       >  2:819494/1‑150 (MQ=255)
                                                                                                                aGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCa                   >  2:493442/1‑150 (MQ=255)
                                                                                                                aGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCa                   >  1:167770/1‑150 (MQ=255)
                                                                                                                aGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCa                   >  1:1191672/1‑150 (MQ=255)
                                                                                                                 ggCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCAt                  <  2:1927201/150‑1 (MQ=255)
                                                                                                                 ggCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCAt                  <  2:1499864/150‑1 (MQ=255)
                                                                                                                   cACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATgg                <  1:2191225/150‑1 (MQ=255)
                                                                                                                   cACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATgg                <  2:1191672/150‑1 (MQ=255)
                                                                                                                      gCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGcc             >  1:1693523/1‑150 (MQ=255)
                                                                                                                        aaGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACgg                                       <  1:235700/122‑1 (MQ=255)
                                                                                                                        aaGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACgg                                       >  2:235700/1‑122 (MQ=255)
                                                                                                                          gATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGa         >  2:1377866/1‑150 (MQ=255)
                                                                                                                          gATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGa         >  2:2049802/1‑150 (MQ=255)
                                                                                                                           aTGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAg        <  2:1796409/150‑1 (MQ=255)
                                                                                                                           aTGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAg        <  2:1307161/150‑1 (MQ=255)
                                                                                                                            tGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGAACTTTTTAACTCATCCATCTGGGCCGGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAgg       <  1:2099489/150‑1 (MQ=255)
                                                                                                                              gAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGtt                          >  2:321698/1‑129 (MQ=255)
                                                                                                                              gAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGtt                          <  1:321698/129‑1 (MQ=255)
                                                                                                                              gAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAg        >  2:2314517/1‑147 (MQ=255)
                                                                                                                              gAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAg        <  1:2314517/147‑1 (MQ=255)
                                                                                                                              gAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGcc     >  1:1949803/1‑150 (MQ=255)
                                                                                                                                ggCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGtt                          >  1:1567288/1‑127 (MQ=255)
                                                                                                                                ggCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGtt                          <  2:1567288/127‑1 (MQ=255)
                                                                                                                                 gCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGtt  <  2:401171/150‑1 (MQ=255)
                                                                                                                                   cgccGGCCAGTTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGaa                                                                             <  1:530325/70‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                         
CGGCGGCGACAGTACGTGACACGTGATCAACTCCATGGCATGCAGCAGCGCCTGGACAGCCGGCCTTCAACGCCGGAAGAAGTGCAGGACCGCGGTGAGTTCGCCCGCTGCCAGGCACGCAAGATGGAGGCCATCGGCCAACTGGCCGGCGGCGTCGCCCATGACTTCAACAACCTGCTCACCAGCATCAGCGGCAGCTTCGAACTCATCCAGCGGCGCCTGCGCCAAGGCCGCAGCGACGGCCTGGATGGCGTTCTGTGCATGGGCCGCGAGGCCGTT  >  NC_002947/3884541‑3884819

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