Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,632,096 G→C intergenic (+658/‑172) emrB → / → RPPX_20645 multidrug resistance protein B/oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,632,0960GC99.4% 566.1 / NA 169intergenic (+658/‑172)emrB/RPPX_20645multidrug resistance protein B/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (90/78);  minor base A (0/1);  total (90/79)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

TGCGCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCC  >  CP009974/4631962‑4632243
                                                                                                                                      |                                                                                                                                                   
                                           cgcACTCCGGCCTTGCGCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCATCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGACCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCa                                                                                          <  1:3442907/149‑1 (MQ=11)
                                           cTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCa                                                                                          <  2:2923435/151‑1 (MQ=2)
                                           cTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAACATCAAAAGCCa                                                                                          >  1:5555610/1‑151 (MQ=9)
                                            tCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAg                                                                                         <  1:3578627/151‑1 (MQ=9)
                                                tCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCa                                                                                     >  1:3038556/1‑151 (MQ=16)
                                                  cGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAAc                                                                                   <  1:1220011/151‑1 (MQ=25)
                                                  cGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAAc                                                                                   <  2:6832863/151‑1 (MQ=255)
                                                      ctggCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAACATCAAGATCAAAAGCCAGATCAACAGcc                                                                               <  1:7233350/148‑1 (MQ=255)
                                                       ttGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCca                                                                              <  2:3253735/151‑1 (MQ=255)
                                                        tGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCACGATCAAAAGCCAGATCAACAGCcac                                                                             >  1:6033989/1‑151 (MQ=255)
                                                        tGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCcac                                                                             <  2:2352084/151‑1 (MQ=255)
                                                         gCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCcaca                                                                            >  2:208330/1‑151 (MQ=255)
                                                         gCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCcaca                                                                            >  1:5655717/1‑151 (MQ=255)
                                                           tCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATc                                                                          <  1:1035482/151‑1 (MQ=255)
                                                           tCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATc                                                                          >  2:1256157/1‑151 (MQ=255)
                                                           tCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATc                                                                          >  1:1391954/1‑151 (MQ=255)
                                                                gCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcag                                                                     <  1:2681883/151‑1 (MQ=255)
                                                                  aGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagca                                                                   <  1:6207778/151‑1 (MQ=255)
                                                                    gACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCTAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCaga                                                                 <  1:5461601/151‑1 (MQ=255)
                                                                    gACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCaga                                                                 >  1:1800871/1‑151 (MQ=255)
                                                                     aCCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCagag                                                                <  2:6865292/151‑1 (MQ=255)
                                                                       cgcgCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAggg                                                              <  1:16943/151‑1 (MQ=255)
                                                                       cgcgCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAggg                                                              >  2:3118326/1‑151 (MQ=255)
                                                                       cgcgCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAggg                                                              <  1:5539135/151‑1 (MQ=255)
                                                                        gcgcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGt                                                             <  1:128805/151‑1 (MQ=255)
                                                                         cgcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCGGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGtt                                                            >  2:6996602/1‑151 (MQ=255)
                                                                          gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           >  1:3243552/1‑151 (MQ=255)
                                                                          gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           <  1:6476793/151‑1 (MQ=255)
                                                                          gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           >  1:236343/1‑151 (MQ=255)
                                                                          gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           <  1:4825397/151‑1 (MQ=255)
                                                                          gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           >  1:382317/1‑151 (MQ=255)
                                                                          gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           >  1:382785/1‑151 (MQ=255)
                                                                          gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           <  2:4076938/151‑1 (MQ=255)
                                                                           ccGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCCGAGGGTTgg                                                          >  1:4412415/1‑151 (MQ=255)
                                                                           ccGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTgg                                                          <  1:3712376/151‑1 (MQ=255)
                                                                            cGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGc                                                         <  2:3953078/151‑1 (MQ=255)
                                                                             gAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCt                                                        <  1:6776415/151‑1 (MQ=255)
                                                                             gAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCt                                                        >  1:4664084/1‑151 (MQ=255)
                                                                              aaGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtg                                                       >  2:1586802/1‑151 (MQ=255)
                                                                               aGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtgt                                                      <  1:1821711/151‑1 (MQ=255)
                                                                                 gCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGc                                                    >  1:862381/1‑151 (MQ=255)
                                                                                  ccccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCt                                                   >  1:2815386/1‑151 (MQ=255)
                                                                                  ccccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCt                                                   <  2:3861846/151‑1 (MQ=255)
                                                                                  ccccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCt                                                   <  2:1022871/151‑1 (MQ=255)
                                                                                   cccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTg                                                  <  2:4084998/151‑1 (MQ=255)
                                                                                    ccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCCGa                                                 >  2:3983716/1‑151 (MQ=255)
                                                                                     cATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGCGTTGGCTATGCTGAt                                                >  2:4441806/1‑151 (MQ=255)
                                                                                      aTCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATc                                               >  2:6006403/1‑151 (MQ=255)
                                                                                       tCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCt                                              <  1:3704426/151‑1 (MQ=255)
                                                                                        ccTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCtt                                             <  1:5335078/151‑1 (MQ=255)
                                                                                            gggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTc                                         >  2:2685315/1‑151 (MQ=255)
                                                                                             ggCCTCAGCGCTTGCCGGCCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCa                                        <  2:6909956/151‑1 (MQ=255)
                                                                                             ggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGCAGTCGCGAAGATCGAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCa                                        <  2:3764473/151‑1 (MQ=255)
                                                                                             ggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCa                                        <  2:2031584/151‑1 (MQ=255)
                                                                                              gCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCAt                                       <  1:845055/151‑1 (MQ=255)
                                                                                              gCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCAt                                       <  1:3494797/151‑1 (MQ=255)
                                                                                                cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTc                                     <  2:6712534/151‑1 (MQ=255)
                                                                                                cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTc                                     <  2:3675225/151‑1 (MQ=255)
                                                                                                 tCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTcc                                    >  1:4189247/1‑151 (MQ=255)
                                                                                                   aGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCaa                                  >  2:4847514/1‑151 (MQ=255)
                                                                                                    gcgcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGATGGCTGTGCTGATCTTCGTCATTCCAAc                                 >  1:2335948/1‑151 (MQ=255)
                                                                                                       cTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCa                              >  1:4989013/1‑151 (MQ=255)
                                                                                                        ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCACAGGGTTGGCTGTGCTGATCCTCGTCATTCCAACGCac                             >  2:4394571/1‑150 (MQ=255)
                                                                                                        ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCCGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTCACCCTCGTCCCTCCAACGCAg                             >  2:4394584/1‑151 (MQ=18)
                                                                                                         tGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgc                            >  2:1307819/1‑151 (MQ=255)
                                                                                                          gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcg                           <  2:2461179/151‑1 (MQ=255)
                                                                                                          gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcc                           >  2:2546561/1‑150 (MQ=255)
                                                                                                           gcGGGCATCCAGGCCCGTCACCTGCCGCCGCAAGGCCTGCGTGCGGCCTCCTGAAGTCGCGAAGATCAAGAGCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcgc                          <  2:6412906/150‑1 (MQ=25)
                                                                                                           aCGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcgc                          >  1:619678/1‑151 (MQ=255)
                                                                                                           aCGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcgc                          <  1:3337701/151‑1 (MQ=255)
                                                                                                            cGGGCATCCATGCCCGGCACCTGCCGCCGGCAGGCCTGCGTTCGGCCCCCTGAAGGCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGcc                         <  1:3750094/151‑1 (MQ=21)
                                                                                                             gggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCa                        >  2:4220712/1‑151 (MQ=255)
                                                                                                              ggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCAt                       <  2:5737992/151‑1 (MQ=255)
                                                                                                              ggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCAt                       <  2:544906/151‑1 (MQ=255)
                                                                                                              ggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCAt                       <  2:5863240/151‑1 (MQ=255)
                                                                                                               gCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCATCGCCAtt                      >  2:4472609/1‑151 (MQ=255)
                                                                                                                cATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCAttt                     >  1:5481304/1‑151 (MQ=255)
                                                                                                                 aTCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTc                    <  2:2982524/151‑1 (MQ=255)
                                                                                                                    cATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGcc                 <  1:4241526/151‑1 (MQ=255)
                                                                                                                     aTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCt                <  2:49071/151‑1 (MQ=255)
                                                                                                                     aTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGAGCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCTATTTCGGCt                >  1:6399387/1‑151 (MQ=21)
                                                                                                                      tGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCCACGCAGCGCCATTTCGCCtc               >  2:7215684/1‑150 (MQ=255)
                                                                                                                      tGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTa               <  1:4000990/151‑1 (MQ=255)
                                                                                                                      tGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTa               >  2:506850/1‑151 (MQ=255)
                                                                                                                         ccGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTg            >  2:447450/1‑151 (MQ=255)
                                                                                                                          cGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGt           <  2:1852423/151‑1 (MQ=255)
                                                                                                                          cGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGt           <  1:5920329/151‑1 (MQ=255)
                                                                                                                           gTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTc          >  1:778691/1‑151 (MQ=255)
                                                                                                                            tCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGcc                                                                               >  1:677771/1‑81 (MQ=255)
                                                                                                                            tCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCg         <  1:1971549/151‑1 (MQ=255)
                                                                                                                            tCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCg         >  2:6639959/1‑151 (MQ=255)
                                                                                                                             cACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGa        <  1:5925537/151‑1 (MQ=255)
                                                                                                                             cACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGa        <  1:726988/151‑1 (MQ=255)
                                                                                                                             cACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGa        >  1:5182328/1‑151 (MQ=255)
                                                                                                                              aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGTCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       <  2:4654880/151‑1 (MQ=255)
                                                                                                                              aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       >  1:1018091/1‑151 (MQ=255)
                                                                                                                              aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       >  1:933804/1‑151 (MQ=255)
                                                                                                                               ccTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATg      <  2:3452910/151‑1 (MQ=255)
                                                                                                                                  gccTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTc   >  2:5564132/2‑151 (MQ=255)
                                                                                                                                  gccTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCCTCATTCCAACGCAGCGCCACTTCGCCTATTCTCgccgccc   >  2:5564152/2‑145 (MQ=25)
                                                                                                                                      |                                                                                                                                                   
TGCGCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCC  >  CP009974/4631962‑4632243

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