Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,173,179 G→T 100% P669P (CCG→CCT rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,173,1790GT95.5% 722.4 / 22.8 243P669P (CCG→CCTrpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (6/5);  new base T (117/115);  total (123/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 243 total aligned reads displayed.

AACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAG  >  CP009273/4173036‑4173324
                                                                                                                                               |                                                                                                                                                 
aaCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTcc                                                                                                                                               <  1:3278041/148‑1 (MQ=255)
  cTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTg                                                                                                                                             >  1:3163784/1‑148 (MQ=255)
     ggATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGaa                                                                                                                                          <  1:3321813/148‑1 (MQ=255)
     ggATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGaa                                                                                                                                          <  2:1810397/148‑1 (MQ=255)
     ggATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGaa                                                                                                                                          <  2:1810211/148‑1 (MQ=255)
      gATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAac                                                                                                                                         <  1:1653209/148‑1 (MQ=255)
        tGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAacac                                                                                                                                       >  1:1846153/1‑148 (MQ=255)
         gaagaaGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACg                                                                                                                                      >  2:3306516/1‑148 (MQ=255)
           agaagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGAt                                                                                                                                    <  1:168204/148‑1 (MQ=255)
             aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGa                                                                                                                                  <  1:1567055/148‑1 (MQ=255)
             aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGa                                                                                                                                  >  2:74168/1‑148 (MQ=255)
             aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGa                                                                                                                                  >  1:2394631/1‑148 (MQ=255)
              agGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGAc                                                                                                                                 >  1:13137/1‑148 (MQ=255)
                         tAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGc                                                                                                                      >  2:2603905/1‑148 (MQ=255)
                          agaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCa                                                                                                                     <  1:1888806/148‑1 (MQ=255)
                          agaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCa                                                                                                                     <  2:1574260/148‑1 (MQ=255)
                          agaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCa                                                                                                                     >  2:1501183/1‑148 (MQ=255)
                          agaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGAAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCa                                                                                                                     <  2:22872/148‑1 (MQ=255)
                           gaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAt                                                                                                                    >  2:249073/1‑148 (MQ=255)
                           gaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAt                                                                                                                    >  2:248792/1‑148 (MQ=255)
                           gaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAt                                                                                                                    >  2:337569/1‑148 (MQ=255)
                           gaagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAt                                                                                                                    >  2:881455/1‑148 (MQ=255)
                            aagaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCAtt                                                                                                                   >  1:2424493/1‑148 (MQ=255)
                              gaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGa                                                                                                                 <  1:561768/148‑1 (MQ=255)
                               aCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGAt                                                                                                                >  2:841470/1‑148 (MQ=255)
                                ccTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATg                                                                                                               <  1:291337/148‑1 (MQ=255)
                                 cTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATgg                                                                                                              >  2:3431022/1‑148 (MQ=255)
                                   ggTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGt                                                                                                            >  1:251915/1‑148 (MQ=255)
                                    gTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTg                                                                                                           <  2:562014/148‑1 (MQ=255)
                                     tAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGc                                                                                                          <  1:986227/148‑1 (MQ=255)
                                       aCTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGa                                                                                                        <  2:1895413/148‑1 (MQ=255)
                                       aCTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGa                                                                                                        <  2:3112425/148‑1 (MQ=255)
                                         ttGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                                                                      >  2:480601/1‑148 (MQ=255)
                                         ttGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAc                                                                                                      <  1:914538/148‑1 (MQ=255)
                                          tGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCa                                                                                                     >  2:1076020/1‑148 (MQ=255)
                                              gTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCa                                                                                                 <  2:581325/148‑1 (MQ=255)
                                                aGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAAc                                                                                               >  1:2108832/1‑148 (MQ=255)
                                                 gCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACg                                                                                              >  2:3384080/1‑148 (MQ=255)
                                                  cAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGt                                                                                             >  1:379166/1‑148 (MQ=255)
                                                   aaaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTc                                                                                            >  2:3035374/1‑148 (MQ=255)
                                                   aaaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTc                                                                                            <  2:1292141/148‑1 (MQ=255)
                                                    aaGGAGAATCCATCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCAGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCa                                                                                           <  2:3095756/148‑1 (MQ=255)
                                                     aGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAg                                                                                          <  1:2174705/148‑1 (MQ=255)
                                                     aGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAg                                                                                          >  1:842105/1‑148 (MQ=255)
                                                      ggCGAATCCAGCTTGTTCAGCCGTGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                                                         >  2:2398202/1‑148 (MQ=255)
                                                      ggCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAgg                                                                                         >  2:757788/1‑148 (MQ=255)
                                                       gCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGc                                                                                        >  1:486990/1‑148 (MQ=255)
                                                              cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCg                                                                                 >  1:2153448/1‑148 (MQ=255)
                                                               aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                                                                <  1:2159809/148‑1 (MQ=255)
                                                                 cTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAct                                                                              <  1:1035628/148‑1 (MQ=255)
                                                                         gCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                                                      <  1:2828459/148‑1 (MQ=255)
                                                                         gCCGCGACCAGGTTGACTACATGGACGTATCCACCCTGCAGGGGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTgcgcgc                                                                      <  2:469022/148‑1 (MQ=255)
                                                                           cgcgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTg                                                                    <  2:2400422/148‑1 (MQ=255)
                                                                            gcgACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGa                                                                   <  1:2327891/148‑1 (MQ=255)
                                                                             cgACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGAt                                                                  >  2:667798/1‑148 (MQ=255)
                                                                             cgACCAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGAt                                                                  >  2:2450863/1‑148 (MQ=255)
                                                                                ccAGGTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAg                                                               >  2:3492069/1‑148 (MQ=255)
                                                                                    gTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCTTTCCGACTCTGCGCGCTGATAAGCCGc                                                           <  2:1004354/148‑1 (MQ=255)
                                                                                    gTTGACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCTTTCCGACTCTGCGCGCTGATAAGCCGc                                                           <  2:2780611/148‑1 (MQ=255)
                                                                                       gACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtgg                                                        <  2:2903715/148‑1 (MQ=255)
                                                                                       gACTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtgg                                                        <  2:961573/148‑1 (MQ=255)
                                                                                        aCTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                                                       >  2:3230604/1‑148 (MQ=255)
                                                                                        aCTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggt                                                       <  1:2906319/148‑1 (MQ=255)
                                                                                         cTACATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggtt                                                      >  2:1593434/1‑148 (MQ=255)
                                                                                           aCATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg                                                    >  2:2436361/1‑148 (MQ=255)
                                                                                            cATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt                                                   <  1:820383/148‑1 (MQ=255)
                                                                                            cATGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttggt                                                   <  1:2603905/148‑1 (MQ=255)
                                                                                             aTGGACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                                                  >  1:2640023/1‑148 (MQ=255)
                                                                                               ggACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt                                                >  1:2753010/1‑148 (MQ=255)
                                                                                                gACGTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg                                               <  1:1498030/148‑1 (MQ=255)
                                                                                                   gTATCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTa                                            <  2:943263/148‑1 (MQ=255)
                                                                                                     aTCCACCCTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATg                                          <  1:337601/148‑1 (MQ=255)
                                                                                                          cccTGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg                                     >  2:1264039/1‑148 (MQ=255)
                                                                                                            ctgcagGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTg                                   <  2:2640609/146‑1 (MQ=255)
                                                                                                            ctgcagGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTg                                   >  2:1625807/3‑148 (MQ=255)
                                                                                                              gcagGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCGGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCt                                 <  1:1959738/148‑1 (MQ=255)
                                                                                                                agGTGGTATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGt                               <  2:1329142/148‑1 (MQ=255)
                                                                                                                 ggtggtATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGtt                              >  2:2343173/1‑148 (MQ=255)
                                                                                                                 ggtggtATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGtt                              >  2:257575/1‑148 (MQ=255)
                                                                                                                 ggtggtATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGtt                              >  2:3530868/1‑148 (MQ=255)
                                                                                                                   tggtATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGc                            >  2:272000/1‑148 (MQ=255)
                                                                                                                     gtATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCg                          <  2:40215/148‑1 (MQ=255)
                                                                                                                      tATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGt                         <  1:249144/148‑1 (MQ=255)
                                                                                                                      tATCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGt                         <  1:248792/148‑1 (MQ=255)
                                                                                                                       aTCCGTCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGtt                        <  1:3086276/148‑1 (MQ=255)
                                                                                                                            tCGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTc                   <  1:3091717/148‑1 (MQ=255)
                                                                                                                             cGGTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTcc                  >  1:3002467/1‑148 (MQ=255)
                                                                                                                               gTGCGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCTGGCGGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCgg                <  1:74168/148‑1 (MQ=255)
                                                                                                                                  cGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGgtgt             <  1:718065/148‑1 (MQ=255)
                                                                                                                                  cGTCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGgtgt             <  1:1778206/148‑1 (MQ=255)
                                                                                                                                    tCCCTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTaa           <  2:1863119/148‑1 (MQ=255)
                                                                                                                                     cccTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAAc          >  1:1078647/1‑148 (MQ=255)
                                                                                                                                       cTGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTg        >  2:2239936/1‑148 (MQ=255)
                                                                                                                                        tGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGc       <  2:637374/148‑1 (MQ=255)
                                                                                                                                        tGATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGc       >  1:1006046/1‑148 (MQ=255)
                                                                                                                                         gATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCg      <  2:2029390/148‑1 (MQ=255)
                                                                                                                                         gATCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCg      <  2:37273/148‑1 (MQ=255)
                                                                                                                                          aTCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCgg     <  1:2672687/148‑1 (MQ=255)
                                                                                                                                          aTCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCgg     <  1:3404718/148‑1 (MQ=255)
                                                                                                                                           tCCCTTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGt    >  1:1215532/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
AACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTTGACTCCGGTGTAACTGCGGTAG  >  CP009273/4173036‑4173324

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: