Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 583,341 T→G 93.6% S33R (AGT→CGT)  ybcH ← PRK09936 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273583,3410TG93.6% 1195.0 / 42.7 312S33R (AGT→CGT) ybcHPRK09936 family protein
Reads supporting (aligned to +/- strand):  ref base T (7/13);  new base G (137/155);  total (144/168)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 312 total aligned reads displayed.

AGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACG  >  CP009273/583195‑583488
                                                                                                                                                  |                                                                                                                                                   
 gCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtc                                                                                                                                                  >  2:4782196/1‑149 (MQ=255)
 gCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtc                                                                                                                                                  <  2:2918092/149‑1 (MQ=255)
 gCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGAtc                                                                                                                                                  >  2:3394995/1‑149 (MQ=255)
       gCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGtt                                                                                                                                            >  2:2722155/1‑149 (MQ=255)
       gCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGtt                                                                                                                                            >  2:382237/1‑149 (MQ=255)
       gCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGtt                                                                                                                                            >  2:3921988/1‑149 (MQ=255)
           gCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCCGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTtgt                                                                                                                                        >  1:2349531/1‑149 (MQ=255)
             ttAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTgg                                                                                                                                      <  2:1748831/149‑1 (MQ=255)
               aaaCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGtt                                                                                                                                    <  1:1675821/149‑1 (MQ=255)
                 aCAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGc                                                                                                                                  <  1:2012202/149‑1 (MQ=255)
                     tAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGa                                                                                                                              >  2:2472877/1‑149 (MQ=255)
                        cGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataa                                                                                                                           <  1:3823302/149‑1 (MQ=255)
                        cGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataa                                                                                                                           <  1:2881231/149‑1 (MQ=255)
                         gTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataat                                                                                                                          <  2:4059366/149‑1 (MQ=255)
                         gTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataat                                                                                                                          <  1:2684856/149‑1 (MQ=255)
                         gTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataat                                                                                                                          <  1:57435/149‑1 (MQ=255)
                         gTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGataat                                                                                                                          <  1:876570/149‑1 (MQ=255)
                           gcgcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATAc                                                                                                                        <  1:679420/149‑1 (MQ=255)
                             gcTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCt                                                                                                                      >  1:829298/1‑149 (MQ=255)
                               tGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCttt                                                                                                                    >  2:3233030/1‑149 (MQ=255)
                               tGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCttt                                                                                                                    <  1:1001319/149‑1 (MQ=255)
                                  tCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCAt                                                                                                                 >  2:2593083/1‑149 (MQ=255)
                                   ctggctggGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATc                                                                                                                <  1:4090113/149‑1 (MQ=255)
                                         gggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                          <  2:4420643/149‑1 (MQ=255)
                                         gggAAAATGCATCGCCGTAACGTGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAg                                                                                                          <  1:4004973/149‑1 (MQ=255)
                                          ggTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAgg                                                                                                         <  1:2821976/149‑1 (MQ=255)
                                             aaaTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaa                                                                                                      >  1:1868355/1‑149 (MQ=255)
                                             aaaTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaa                                                                                                      <  2:4320088/149‑1 (MQ=255)
                                               aTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaaaa                                                                                                    >  2:5223924/1‑149 (MQ=255)
                                               aTGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGaaaaa                                                                                                    <  2:1868603/149‑1 (MQ=255)
                                                  cATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAggg                                                                                                 >  1:1616169/1‑149 (MQ=255)
                                                    tCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCt                                                                                               <  1:2392395/149‑1 (MQ=255)
                                                     cGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTg                                                                                              <  1:4782218/149‑1 (MQ=255)
                                                     cGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTg                                                                                              <  1:1016611/149‑1 (MQ=255)
                                                     cGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTg                                                                                              <  1:3394995/149‑1 (MQ=255)
                                                      gCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGa                                                                                             <  1:4675736/149‑1 (MQ=255)
                                                      gCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGa                                                                                             >  1:2080747/1‑149 (MQ=255)
                                                       ccGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGAc                                                                                            <  2:4705020/149‑1 (MQ=255)
                                                         gTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCa                                                                                          >  2:4528788/1‑149 (MQ=255)
                                                            aCGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaaaa                                                                                       <  2:2763009/149‑1 (MQ=255)
                                                            aCGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaaaa                                                                                       >  1:4590140/1‑149 (MQ=255)
                                                            aCGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCaaaa                                                                                       >  1:3461731/1‑149 (MQ=255)
                                                             cGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAg                                                                                      <  2:1331511/149‑1 (MQ=255)
                                                              gggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGc                                                                                     >  2:664470/1‑149 (MQ=255)
                                                               ggTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCa                                                                                    <  1:4235033/149‑1 (MQ=255)
                                                                          aaaCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAAt                                                                         >  1:4436906/1‑149 (MQ=255)
                                                                          aaaCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAAt                                                                         <  2:865799/149‑1 (MQ=255)
                                                                           aaCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATa                                                                        <  1:3362762/149‑1 (MQ=255)
                                                                            aCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATAc                                                                       <  2:31598/149‑1 (MQ=255)
                                                                            aCAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATAc                                                                       <  2:671814/149‑1 (MQ=255)
                                                                               aGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCACTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGaa                                                                    >  1:259693/1‑149 (MQ=255)
                                                                                 ggTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGaaaa                                                                  <  2:1137067/149‑1 (MQ=255)
                                                                                 ggTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGaaaa                                                                  <  2:1631590/149‑1 (MQ=255)
                                                                                 ggTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGaaaa                                                                  <  2:4799802/149‑1 (MQ=255)
                                                                                  gTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAAt                                                                 >  2:4287411/1‑149 (MQ=255)
                                                                                  gTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAAt                                                                 <  1:44825/149‑1 (MQ=255)
                                                                                  gTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAAt                                                                 <  2:3693657/149‑1 (MQ=255)
                                                                                    aTCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGa                                                               <  2:1002776/149‑1 (MQ=255)
                                                                                      cGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAAc                                                             >  1:2534800/1‑149 (MQ=255)
                                                                                      cGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAAc                                                             >  1:4896608/1‑149 (MQ=255)
                                                                                       gAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACt                                                            >  2:3540952/1‑149 (MQ=255)
                                                                                       gAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACt                                                            <  2:4348131/149‑1 (MQ=255)
                                                                                        aaGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACtt                                                           >  2:545341/1‑149 (MQ=255)
                                                                                            cTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGc                                                       >  2:1425047/1‑149 (MQ=255)
                                                                                             ttGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCa                                                      >  2:2620691/1‑149 (MQ=255)
                                                                                              tGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAt                                                     <  1:3667627/149‑1 (MQ=255)
                                                                                               gCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAtt                                                    <  2:631509/149‑1 (MQ=255)
                                                                                                cAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCAttt                                                   <  1:927699/149‑1 (MQ=255)
                                                                                                 aaaCGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTa                                                  <  1:4333991/149‑1 (MQ=255)
                                                                                                  aaCGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTAc                                                 >  1:2284547/1‑149 (MQ=255)
                                                                                                        aCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGc                                           <  1:614266/149‑1 (MQ=255)
                                                                                                        aCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGc                                           >  2:3896271/1‑149 (MQ=255)
                                                                                                        aCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGc                                           >  2:4178976/1‑149 (MQ=255)
                                                                                                        aCTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGc                                           <  1:264611/149‑1 (MQ=255)
                                                                                                            aCTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTc                                       >  1:327850/1‑149 (MQ=255)
                                                                                                             cTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGTTTTTGTGGTTGCCAGATAATACCTTTCATCTCAAAGGAAAAATGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTcc                                      <  2:206545/149‑1 (MQ=255)
                                                                                                               catcaGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGCTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCaa                                    <  1:5100459/149‑1 (MQ=255)
                                                                                                                  caGCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATgg                                 <  1:1173091/149‑1 (MQ=255)
                                                                                                                    gCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTg                               >  1:5043733/1‑149 (MQ=255)
                                                                                                                    gCCCCTGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTg                               >  1:4931354/1‑149 (MQ=255)
                                                                                                                     ccccTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTga                              <  1:3685487/149‑1 (MQ=255)
                                                                                                                         tGCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGaaa                          <  2:1703545/149‑1 (MQ=255)
                                                                                                                          gCCACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAc                         >  2:2217685/1‑149 (MQ=255)
                                                                                                                           ccACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAcg                        <  1:3351192/149‑1 (MQ=255)
                                                                                                                            cACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAcgc                       <  1:3707230/149‑1 (MQ=255)
                                                                                                                            cACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAcgc                       >  1:357163/1‑149 (MQ=255)
                                                                                                                            cACTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAAcgc                       >  1:1381119/1‑149 (MQ=255)
                                                                                                                              cTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgt                     <  1:1122182/149‑1 (MQ=255)
                                                                                                                              cTGGGTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgt                     >  1:916225/1‑149 (MQ=255)
                                                                                                                                 ggTATCGGTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCgttgt                  <  1:1520702/149‑1 (MQ=255)
                                                                                                                                     tCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTgcg              >  2:3292834/1‑149 (MQ=255)
                                                                                                                                        gTAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGctc           <  1:3243095/149‑1 (MQ=255)
                                                                                                                                         tAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGctct          <  2:123444/149‑1 (MQ=255)
                                                                                                                                         tAACCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGctct          <  1:2472877/149‑1 (MQ=255)
                                                                                                                                          aaCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCtctc         <  2:1845530/149‑1 (MQ=255)
                                                                                                                                           aCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTcc        <  2:4993682/149‑1 (MQ=255)
                                                                                                                                           aCCTGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTcc        >  2:5160793/1‑149 (MQ=255)
                                                                                                                                              tGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGtt     <  1:3233030/149‑1 (MQ=255)
                                                                                                                                              tGACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGtt     <  1:2571537/149‑1 (MQ=255)
                                                                                                                                               gACGATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTa    >  1:4980474/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCTGCGGCCCGCTTAAACAATAACGTGCGCTGTTCTGGCTGGGTAAATGCATCGCCGTAACGGGTCCATTGCAAAACAAGGGTATCGAAGCCTTGCAAACGTAACTGACTCATCAGCCCCTGCCACTGGGTATCGGTAACCTGACTATCTCGGTTTTGTGGTTGCCAGATAATACCTTTCATCGCAAAGGAAAAAGGGCTGACCAAAAGCAGTGTCAGCAATACGAAAATGAACTTACGCATTTACCAGTGCACTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACG  >  CP009273/583195‑583488

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