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.1% 360.0 / ‑4.8 105intergenic (+658/‑172)emrB/RPPX_20645multidrug resistance protein B/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (59/45);  minor base A (0/1);  total (59/46)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 105 total aligned reads displayed.

CACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCC  >  CP009974/4632007‑4632243
                                                                                         |                                                                                                                                                   
cACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGa                                                                                        <  2:2071464/151‑1 (MQ=12)
cACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGa                                                                                        <  2:2069710/151‑1 (MQ=12)
cACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGa                                                                                        >  1:2210295/1‑151 (MQ=9)
 aCTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGAt                                                                                       >  2:148870/1‑151 (MQ=9)
 aCTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGAt                                                                                       <  1:1795435/151‑1 (MQ=12)
  cTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAATGCCTGCGTTCGGCCTCCTTAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATc                                                                                      >  1:4109753/1‑151 (MQ=255)
  cTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATc                                                                                      <  1:2088306/151‑1 (MQ=18)
  cTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATc                                                                                      <  1:2092973/151‑1 (MQ=16)
    ccGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCaa                                                                                    >  1:3160507/1‑151 (MQ=18)
       gCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAg                                                                                 <  1:5066247/151‑1 (MQ=255)
       gCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTCGCCGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGCCCGGCCTCCTGCAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAg                                                                                 >  2:4938357/1‑151 (MQ=11)
        ccTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGc                                                                                >  2:1559745/1‑151 (MQ=255)
         cTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGcc                                                                               >  1:3588148/1‑151 (MQ=255)
          ttGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCca                                                                              >  1:2813537/1‑151 (MQ=255)
            gCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGGAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCcaca                                                                            >  2:254758/1‑151 (MQ=255)
            gCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCcaca                                                                            >  2:3888103/1‑151 (MQ=255)
               cccGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGATCGGCCTCCTGAAGTCGCGAAGATCAAGATAAAGATCAAAAGCCAGATCAACAGCCacaaca                                                                         >  2:1876946/1‑151 (MQ=255)
                 cGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAc                                                                       >  1:140166/1‑151 (MQ=255)
                  ggCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCACGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAca                                                                      >  1:1048955/1‑151 (MQ=255)
                  ggCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCGTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAca                                                                      <  2:5046367/151‑1 (MQ=255)
                    cAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagc                                                                    <  1:2472645/151‑1 (MQ=255)
                     aGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagca                                                                   >  1:5100455/1‑151 (MQ=255)
                     aGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagca                                                                   >  1:525651/1‑151 (MQ=255)
                      ggACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAAcagcag                                                                  >  2:4468581/1‑151 (MQ=255)
                       gACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCaga                                                                 <  1:3322712/151‑1 (MQ=255)
                       gACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCaga                                                                 >  1:1508102/1‑151 (MQ=255)
                          cgcgCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAggg                                                              >  1:4007662/1‑151 (MQ=255)
                           gcgcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGt                                                             >  2:4172342/1‑151 (MQ=255)
                             gcCGAAGGCCCCATCCTGGGGCGTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           >  2:4653278/1‑151 (MQ=255)
                             gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           <  2:1944261/151‑1 (MQ=255)
                             gcCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTg                                                           >  1:3962324/1‑151 (MQ=255)
                              ccGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTgg                                                          <  1:3563638/151‑1 (MQ=255)
                               cGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGc                                                         <  1:3129743/151‑1 (MQ=255)
                                gAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCt                                                        >  1:5118537/1‑151 (MQ=255)
                                gAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCt                                                        >  1:3526191/1‑151 (MQ=255)
                                 aaGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtg                                                       <  2:4367345/151‑1 (MQ=255)
                                 aaGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtg                                                       <  2:317184/151‑1 (MQ=255)
                                  aGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtgt                                                      <  2:1807477/151‑1 (MQ=255)
                                  aGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtgt                                                      >  1:4155232/1‑151 (MQ=255)
                                  aGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCtgt                                                      <  1:2675570/151‑1 (MQ=255)
                                    gCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGc                                                    >  2:1999599/1‑151 (MQ=255)
                                    gCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGc                                                    <  2:2780505/151‑1 (MQ=255)
                                    gCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGACTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGc                                                    <  2:4085522/151‑1 (MQ=255)
                                     ccccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCt                                                   <  1:3881590/151‑1 (MQ=255)
                                     ccccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCt                                                   <  1:2860121/151‑1 (MQ=255)
                                      cccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTg                                                  >  1:4549029/1‑151 (MQ=255)
                                      cccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGCGTTGGCTGTGCTg                                                  >  1:5043365/1‑151 (MQ=255)
                                       ccATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGa                                                 >  2:1043490/1‑151 (MQ=255)
                                         aTCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATc                                               <  1:2471276/151‑1 (MQ=255)
                                         aTCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATc                                               <  1:1380664/151‑1 (MQ=255)
                                          tCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCt                                              >  2:4639639/1‑151 (MQ=255)
                                           ccTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCtt                                             <  1:3758634/151‑1 (MQ=255)
                                           ccTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCtt                                             <  2:3262749/151‑1 (MQ=255)
                                            cTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTc                                            <  1:3851457/151‑1 (MQ=255)
                                            cTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTc                                            >  1:3759767/1‑151 (MQ=255)
                                             gggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCg                                           <  2:1218351/150‑1 (MQ=255)
                                              ggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGt                                          <  2:2963515/151‑1 (MQ=255)
                                              ggggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAAAAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGt                                          <  1:566569/151‑1 (MQ=255)
                                               gggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTc                                         <  2:1237357/151‑1 (MQ=255)
                                                ggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGCTCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCACAGGGTTGGCTGTGCTGATCTTCGTCa                                        >  2:5099876/1‑151 (MQ=255)
                                                ggCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCa                                        >  2:4171906/1‑151 (MQ=255)
                                                 gCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCAt                                       <  2:1918192/151‑1 (MQ=255)
                                                  ccTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCAtt                                      >  2:3018026/1‑151 (MQ=255)
                                                   cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTc                                     >  1:4913475/1‑151 (MQ=255)
                                                   cTCAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTc                                     >  2:1647948/1‑151 (MQ=255)
                                                     cAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCa                                   <  2:3966038/151‑1 (MQ=255)
                                                     cAGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCa                                   <  1:2583878/151‑1 (MQ=255)
                                                      ggcgcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCaa                                  <  2:214580/150‑1 (MQ=255)
                                                      aGCGCTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCaa                                  <  1:3849624/151‑1 (MQ=255)
                                                      aGCGCTTGACGGGCATCCATGCCCGCCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCACGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCAGTCCaa                                  >  2:3564696/1‑151 (MQ=18)
                                                       gcgcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAAc                                 >  1:2675628/1‑151 (MQ=255)
                                                        cgcTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGCGTTGGCTGTGCTGATCTTCGTCATTCCTACg                                >  2:4340517/1‑151 (MQ=255)
                                                          cTTGACGGGTATCCATGCCCGTCACCTGCCTACGCAATGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCACGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCa                              <  2:4092000/151‑1 (MQ=25)
                                                          cTTGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCa                              <  1:2401209/151‑1 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCac                             >  2:501157/1‑150 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             >  1:1183061/1‑151 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             <  1:3209819/151‑1 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             >  1:477439/1‑151 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             >  1:1443025/1‑151 (MQ=255)
                                                           ttGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAg                             >  1:5061382/1‑151 (MQ=255)
                                                            tGACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgc                            >  1:1554659/1‑151 (MQ=255)
                                                             gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcg                           >  1:2456041/1‑151 (MQ=255)
                                                             gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcg                           >  1:3096590/1‑151 (MQ=255)
                                                             gACGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcg                           <  1:4146292/151‑1 (MQ=255)
                                                              aCGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAgcgc                          >  1:1125540/1‑151 (MQ=255)
                                                               cGGGCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGcc                         >  1:231518/1‑151 (MQ=255)
                                                                gggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCa                        >  2:3907847/1‑151 (MQ=255)
                                                                gggCATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCa                        >  2:819589/1‑151 (MQ=255)
                                                                  gcatggaTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCAtt                      >  1:2007793/7‑151 (MQ=255)
                                                                   cATCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCAttt                     >  1:1924779/1‑151 (MQ=255)
                                                                    aTCCATGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTc                    >  2:3309056/1‑151 (MQ=255)
                                                                       gaTGCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGcc                 >  2:5072630/2‑151 (MQ=255)
                                                                          gCCCGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTAt              <  2:525651/151‑1 (MQ=255)
                                                                           cccGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGCCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCACCGCCATTTCGCCTAtt             >  2:2471888/1‑151 (MQ=255)
                                                                             cGTCACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGt           <  2:610928/151‑1 (MQ=255)
                                                                                cACCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGa        <  2:217213/151‑1 (MQ=255)
                                                                                 aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       <  1:4468581/151‑1 (MQ=255)
                                                                                 aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       >  2:4685912/1‑151 (MQ=255)
                                                                                 aCCTGCCTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGAt       <  1:1912666/151‑1 (MQ=255)
                                                                                      ccTCCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTcc  <  1:231418/151‑1 (MQ=255)
                                                                                         |                                                                                                                                                   
CACTCCGGCCTTGCTCCCGGCAGGACCGCGCCGAAGGCCCCATCCTGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAAGATCAAAAGCCAGATCAACAGCCACATCAACAGCAGAGGGTTGGCTGTGCTGATCTTCGTCATTCCAACGCAGCGCCATTTCGCCTATTGTCGATGCTCC  >  CP009974/4632007‑4632243

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