Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,698,213 G→T 7.1% R264L (CGT→CTT)  hslU → protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,698,2130GT 7.1% 312.1 / 9.2 99R264L (CGT→CTT) hslUprotease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base G (49/43);  new base T (3/4);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

AAGCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCTGA  >  NC_002947/5698065‑5698359
                                                                                                                                                    |                                                                                                                                                   
atcttcaaGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAACGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGt                                                                                                                                                    <  1:919829/144‑1 (MQ=255)
aaGCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGt                                                                                                                                                    <  1:2186642/150‑1 (MQ=255)
 aGCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTg                                                                                                                                                   <  1:301031/150‑1 (MQ=255)
  acTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTgg                                                                                                                                                  <  1:2217329/149‑1 (MQ=255)
     gAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGATATGGTGTGCATAGAGGAAATCGACAAGGTTGCCAAGCGTgagat                                                                                                                                               >  2:919829/1‑146 (MQ=255)
         gTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCACGCTTGGCAACCtt                                                                                                                                           <  1:1070137/150‑1 (MQ=255)
         gTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGtt                                                                                                                                           >  1:276823/1‑150 (MQ=255)
         gTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAAGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAACGTTGCCACGCTTGGCAACCtt                                                                                                                                           <  1:163450/150‑1 (MQ=255)
            aaaGAAGCGCTGAAGATGGTTTGCGTTTAAGAAGCGAGCCTCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTggc                                                                                                                                        >  1:1261494/1‑150 (MQ=255)
            aaaGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTCggc                                                                                                                                        >  2:1628805/1‑150 (MQ=255)
                     cTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTgg                                                                                                                                                  <  1:78459/131‑1 (MQ=255)
                     cTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTgg                                                                                                                                                  >  2:78459/1‑131 (MQ=255)
                       gAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACgt                                                                                                                             >  2:41697/1‑150 (MQ=255)
                            tGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGTAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTccc                                                                                                                        <  2:141521/150‑1 (MQ=255)
                              gTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGt                                                                                                                      <  1:521106/150‑1 (MQ=255)
                               ttCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTg                                                                                                                                                   >  1:266934/1‑120 (MQ=255)
                               ttCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTg                                                                                                                                                   <  2:266934/120‑1 (MQ=255)
                                 cgcgATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGaa                                                                                                                   <  2:2338386/150‑1 (MQ=255)
                                    gATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGAc                                                                                                                               <  1:10666/135‑1 (MQ=255)
                                    gATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGAc                                                                                                                               >  2:10666/1‑135 (MQ=255)
                                    gATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGc                                                                                                                <  1:1682435/150‑1 (MQ=255)
                                         agaagCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACa                                                                                                           <  2:1828277/150‑1 (MQ=255)
                                         agaagCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACa                                                                                                           >  2:314603/1‑150 (MQ=255)
                                         agaagCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACa                                                                                                           <  2:304236/150‑1 (MQ=255)
                                              cGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGcgcg                                                                                                     >  1:709111/1‑150 (MQ=255)
                                               gAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGc                                                                                                                >  2:1977639/1‑139 (MQ=255)
                                               gAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGa                                                                                                    >  2:1065248/1‑150 (MQ=255)
                                               gAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCAAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGc                                                                                                                <  1:1977639/139‑1 (MQ=255)
                                                 gccgccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGAcc                                                                                                  <  2:1144932/150‑1 (MQ=255)
                                                 gccgccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGAcc                                                                                                  <  2:2380851/150‑1 (MQ=255)
                                                    gccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACctgc                                                                                               >  1:920638/1‑150 (MQ=255)
                                                     ccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTgg                                                                                                                                         >  2:1881744/1‑108 (MQ=255)
                                                     ccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTgg                                                                                                                                         <  1:1881744/108‑1 (MQ=255)
                                                          tCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGc                                                                                         <  1:314603/150‑1 (MQ=255)
                                                          tCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGc                                                                                         <  2:1261494/150‑1 (MQ=255)
                                                           cAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCt                                                                                        <  1:1882632/150‑1 (MQ=255)
                                                            aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGaa                                                                                                                   <  1:298726/123‑1 (MQ=255)
                                                            aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGaa                                                                                                                   >  2:298726/1‑123 (MQ=255)
                                                            aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTg                                                                                       >  2:2382365/1‑150 (MQ=255)
                                                            aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTg                                                                                       >  1:2207240/1‑150 (MQ=255)
                                                             aCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                      <  1:1414752/150‑1 (MQ=255)
                                                                aGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCg                                                                                   >  2:1942422/1‑150 (MQ=255)
                                                                  gAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGaa                                                                                 <  2:797471/150‑1 (MQ=255)
                                                                  gAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGaa                                                                                 >  1:482737/1‑150 (MQ=255)
                                                                      agCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTg                                                                                                                                                   <  1:1782613/81‑1 (MQ=255)
                                                                      agCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTg                                                                                                                                                   >  2:1782613/1‑81 (MQ=255)
                                                                      agCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCt                                                                             <  2:1975069/150‑1 (MQ=255)
                                                                         tCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCa                                                                          >  1:1664253/1‑150 (MQ=255)
                                                                          caaggccaaggcCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCAc                                                                         >  1:1065072/1‑150 (MQ=255)
                                                                           aaggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTg                                                                            >  2:1284067/1‑146 (MQ=255)
                                                                           aaggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTg                                                                            <  1:1284067/146‑1 (MQ=255)
                                                                           aaggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGCCGCAACTTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCAcc                                                                        >  1:1497102/1‑150 (MQ=255)
                                                                            aggccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCAcca                                                                       >  1:481054/1‑149 (MQ=255)
                                                                                caaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCaa                                                                   >  1:1198206/1‑150 (MQ=255)
                                                                                   ggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAactc                                                                <  2:655169/150‑3 (MQ=255)
                                                                                       cTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGTCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAg                                                            >  1:914198/1‑150 (MQ=255)
                                                                                       cTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACctgctg                                                                                             <  2:2163762/117‑1 (MQ=255)
                                                                                       cTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACctgctg                                                                                             >  1:2163762/1‑117 (MQ=255)
                                                                                       cTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                      <  1:2364096/124‑1 (MQ=255)
                                                                                       cTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                      >  2:2364096/1‑124 (MQ=255)
                                                                                       cTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAg                                                            >  2:1114777/1‑150 (MQ=255)
                                                                                          gAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTg                                                         <  2:1065072/150‑1 (MQ=255)
                                                                                                   gAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTggc                                                                                                                                        >  2:1001136/1‑63 (MQ=255)
                                                                                                   gAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTggc                                                                                                                                        <  1:1001136/63‑1 (MQ=255)
                                                                                                   gAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTc                                                <  1:1942422/150‑1 (MQ=255)
                                                                                                   gAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTc                                                >  2:2147752/1‑150 (MQ=255)
                                                                                                    aaCAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCa                                               >  2:1847987/1‑150 (MQ=255)
                                                                                                     aCAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCaa                                              >  2:1659885/1‑150 (MQ=255)
                                                                                                        gCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgac                                           <  1:1114777/150‑1 (MQ=255)
                                                                                                         cACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgacc                                          >  1:986977/1‑150 (MQ=255)
                                                                                                            ggTATCGTGTTCATCGACGAAATCGACAAGGTTGCCCAGGGTGGGAACTTTGTCGTGGTGGACTTTTCCCATGTAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgaccgac                                       >  1:1855516/1‑150 (MQ=255)
                                                                                                               aTCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACcac                                    >  2:1616458/1‑150 (MQ=255)
                                                                                                                   tgtTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTc                                >  1:2191123/1‑150 (MQ=255)
                                                                                                                   tgtTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTc                                <  2:2191123/150‑1 (MQ=255)
                                                                                                                      tCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGt                             <  1:2280990/150‑1 (MQ=255)
                                                                                                                        aTCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGttg                           >  2:1117038/1‑149 (MQ=255)
                                                                                                                        aTCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTc                           <  2:482737/150‑1 (MQ=255)
                                                                                                                         tCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCAGCGCCGACGGGTCCCGTCAAGGCGTACAGAGCGGAACGGTCAGTACTATATAAGCCGGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGAACGCAGTCTGTTCa                          >  1:1151102/1‑150 (MQ=255)
                                                                                                                             cgaACTCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGc                      >  2:1034380/1‑150 (MQ=255)
                                                                                                                                  tCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCgg                 >  2:797314/1‑150 (MQ=255)
                                                                                                                                  tCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCgg                 <  2:1855516/150‑1 (MQ=255)
                                                                                                                                    gACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCCTTGCCTCGGGt               >  1:123320/1‑150 (MQ=255)
                                                                                                                                    gACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGt               <  2:1151102/150‑1 (MQ=255)
                                                                                                                                     aCAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGt               <  1:953701/149‑1 (MQ=255)
                                                                                                                                     aCAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGt               >  2:953701/1‑149 (MQ=255)
                                                                                                                                      cAATGTTACCAAGCGTGGCAACGTTGTCGGCGCCGACGTGTCCCGTGAAGGCGGACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCGTGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGc             <  1:2299971/150‑1 (MQ=255)
                                                                                                                                       aacgTTGCCACGCTTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGt                             <  1:842231/130‑1 (MQ=255)
                                                                                                                                       aacgTTGCCACGCTTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGt                             >  2:842231/4‑133 (MQ=255)
                                                                                                                                       aacgTTGCCACGCTTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCa            >  1:2230988/4‑150 (MQ=255)
                                                                                                                                       aaGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCAGGTgga            >  2:939813/1‑148 (MQ=255)
                                                                                                                                          gTTGCCACGCTTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAga                                            >  1:1964759/1‑115 (MQ=255)
                                                                                                                                          gTTGCCACGCTTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAga                                            <  2:1964759/115‑1 (MQ=255)
                                                                                                                                          gTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTc         >  2:1328397/1‑150 (MQ=255)
                                                                                                                                             gCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCAc      >  1:1787637/1‑150 (MQ=255)
                                                                                                                                               cAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACCTTCTGTTCATTTCCTCGGGTGCATTCCACCt    >  2:492743/1‑150 (MQ=255)
                                                                                                                                                aaGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATgg                                                  <  1:2293730/103‑1 (MQ=255)
                                                                                                                                                aaGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATgg                                                  >  2:2293730/1‑103 (MQ=255)
                                                                                                                                                 aGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCTGa  <  2:1664253/150‑1 (MQ=255)
                                                                                                                                                 aGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCTGa  <  2:1497102/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAGCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCTGA  >  NC_002947/5698065‑5698359

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