Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,756,930 A→T L392H (CTT→CAT)  RPPX_16780 ← phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,756,9300AT99.3% 1039.5 / ‑1.6 278L392H (CTT→CAT) RPPX_16780phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (142/134);  total (143/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 278 total aligned reads displayed.

TCGGCAATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGCAC  >  CP009974/3756783‑3757075
                                                                                                                                                   |                                                                                                                                                 
 cGGCAATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGat                                                                                                                                               <  2:2058459/151‑1 (MQ=255)
  ggCAATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGata                                                                                                                                              <  2:761234/151‑1 (MQ=255)
  ggCAATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGata                                                                                                                                              >  1:4801325/1‑151 (MQ=255)
   gCAATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGatac                                                                                                                                             <  2:2514870/151‑2 (MQ=255)
   gCAATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGatac                                                                                                                                             <  1:3298378/151‑2 (MQ=255)
     aaTCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCg                                                                                                                                           <  2:3284979/151‑1 (MQ=255)
     aaTCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCg                                                                                                                                           <  1:2110508/151‑1 (MQ=255)
            aaTGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCAt                                                                                                                                    >  2:2152980/1‑151 (MQ=255)
                  aaCTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGc                                                                                                                              <  1:5063527/151‑1 (MQ=255)
                   aCTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGcc                                                                                                                             >  1:5368904/1‑151 (MQ=255)
                   aCTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGcc                                                                                                                             <  1:729460/151‑1 (MQ=255)
                   aCTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGcc                                                                                                                             >  2:1418198/1‑151 (MQ=255)
                    cTCGAGCGACTGATCGGGATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGccc                                                                                                                            <  2:5907966/151‑1 (MQ=255)
                    cTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGccc                                                                                                                            >  1:5979769/1‑151 (MQ=255)
                    cTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGccc                                                                                                                            <  2:3317892/151‑1 (MQ=255)
                     tCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTTCCCt                                                                                                                           >  1:465215/1‑151 (MQ=255)
                     tCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGcccc                                                                                                                           >  2:445901/1‑150 (MQ=255)
                     tCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGGGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCt                                                                                                                           <  2:5007383/151‑1 (MQ=255)
                            aCTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTg                                                                                                                    <  2:6263421/151‑1 (MQ=255)
                             cTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTgg                                                                                                                   >  1:5256315/1‑151 (MQ=255)
                                aTCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGgtgt                                                                                                                >  2:465199/1‑151 (MQ=255)
                                aTCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGgtgt                                                                                                                >  2:1545379/1‑151 (MQ=255)
                                 tCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTa                                                                                                               <  1:6530539/151‑1 (MQ=255)
                                    gCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACtt                                                                                                            <  2:2865146/151‑1 (MQ=255)
                                           tCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGgaagaa                                                                                                     <  2:1127781/151‑1 (MQ=255)
                                             cgcgGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACt                                                                                                   <  2:5588880/151‑1 (MQ=255)
                                             cgcgGAACACAATGGGTGTGAAGCCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACt                                                                                                   <  2:5083950/151‑1 (MQ=255)
                                              gcgGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTg                                                                                                  >  1:1887500/1‑151 (MQ=255)
                                               cgGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGGTGCATTCGGTGACGTTCTGCCCAGGCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGc                                                                                                 <  2:5140588/151‑1 (MQ=255)
                                                ggATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGCAGAACTGCt                                                                                                >  2:6128246/1‑151 (MQ=255)
                                                  aTCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTc                                                                                              <  2:1557715/151‑1 (MQ=255)
                                                     acaATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGTGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCAcc                                                                                           <  1:179012/151‑1 (MQ=255)
                                                      caATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCt                                                                                          <  1:702236/151‑1 (MQ=255)
                                                      caATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAACAACTGCTTCAccc                                                                                          >  2:3508558/1‑150 (MQ=255)
                                                          gggTGTGCTACCGGTAAGACACCGCACCGCCGGTGGTGGGCCGGGCACCGCCGATGCAGTCGGTGACGTTCTGCCGCGTCCACTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATc                                                                                      <  2:4067643/151‑1 (MQ=255)
                                                          gggTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATc                                                                                      >  2:928463/1‑151 (MQ=255)
                                                          gggTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATc                                                                                      >  1:1242482/1‑151 (MQ=255)
                                                            gtgtGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCa                                                                                    <  2:343716/151‑1 (MQ=255)
                                                              gtgATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGa                                                                                  <  2:3744330/151‑1 (MQ=255)
                                                                gATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAAt                                                                                <  1:4415021/151‑1 (MQ=255)
                                                                  tACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCt                                                                              >  2:3774466/1‑151 (MQ=255)
                                                                    ccGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTgc                                                                            >  1:1031965/1‑151 (MQ=255)
                                                                    ccGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTgc                                                                            <  2:3595866/151‑1 (MQ=255)
                                                                     cGGTCAGACAACGCATCCTCGGTGACCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTGCTTGGAAGAATTGCTTTACCTCATCCATAATCTtcg                                                                           >  2:616255/1‑151 (MQ=255)
                                                                      ggTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTgcgc                                                                          <  1:2054652/151‑1 (MQ=255)
                                                                      ggTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTgcgc                                                                          <  1:4062206/151‑1 (MQ=255)
                                                                       gTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCa                                                                         >  1:5927586/1‑151 (MQ=255)
                                                                       gTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCa                                                                         >  2:3544160/1‑151 (MQ=255)
                                                                        tCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCaa                                                                        >  2:1921302/1‑151 (MQ=255)
                                                                          aGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCaaaa                                                                      >  1:6098281/1‑151 (MQ=255)
                                                                          aGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCaaaa                                                                      >  1:3633709/1‑151 (MQ=255)
                                                                           gACAACGCATCCTGGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAAt                                                                     <  2:5699640/151‑1 (MQ=255)
                                                                              aaCGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGGACATGGAAGAACTGCTTCACCTCATCCACAGTCTGCGCAAAACcgc                                                                  >  2:4293879/1‑151 (MQ=255)
                                                                                 gCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCggg                                                               >  2:6420594/1‑151 (MQ=255)
                                                                                  cATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGt                                                              >  2:2037621/1‑151 (MQ=255)
                                                                                   aTCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTc                                                             <  1:4542952/151‑1 (MQ=255)
                                                                                     ccTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCtt                                                           >  1:6402776/1‑151 (MQ=255)
                                                                                     ccTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCtt                                                           >  1:3094321/1‑151 (MQ=255)
                                                                                     ccTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCtt                                                           >  2:6330579/1‑151 (MQ=255)
                                                                                      ttCGGTGATCGGCCGGCACCGCCGATGCATTCGGTGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAAcgc                                                                                                  >  2:4715174/2‑108 (MQ=255)
                                                                                       tCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGt                                                         >  1:1466647/1‑151 (MQ=255)
                                                                                          gTGATCGGCCGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGTAAGAACTGCTTCACCTCAGCCAGAATCTGCGCAAAATCGCGGGTCTTGTAAc                                                      >  2:1725599/1‑151 (MQ=255)
                                                                                                  ccGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGTCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGctgctg                                              >  2:114391/1‑151 (MQ=255)
                                                                                                   cGGGCACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTgg                                             <  1:4180341/151‑1 (MQ=255)
                                                                                                    gggCCCCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGc                                            <  2:6536526/151‑1 (MQ=255)
                                                                                                       cACCGCCGATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCtt                                         >  1:3888277/1‑151 (MQ=255)
                                                                                                         ccgccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGa                                       >  1:2103366/1‑151 (MQ=255)
                                                                                                         ccgccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGa                                       <  1:3727303/151‑1 (MQ=255)
                                                                                                         ccgccgATGCATTCGGTGACGTACTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGa                                       >  2:5812953/1‑151 (MQ=255)
                                                                                                          cgccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGAt                                      >  1:2415909/1‑151 (MQ=255)
                                                                                                           gccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATg                                     <  2:4921092/151‑1 (MQ=255)
                                                                                                           gccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATg                                     >  2:2825216/1‑151 (MQ=255)
                                                                                                            ccgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATgg                                    >  1:5495807/1‑151 (MQ=255)
                                                                                                             cgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGgt                                   <  1:6196538/151‑1 (MQ=255)
                                                                                                             cgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGgt                                   <  1:4938297/151‑1 (MQ=255)
                                                                                                             cgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGgt                                   >  2:4483148/1‑151 (MQ=255)
                                                                                                             cgATGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGgt                                   >  1:1751625/1‑151 (MQ=255)
                                                                                                                tGCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGtt                                <  2:2706257/151‑1 (MQ=255)
                                                                                                                 gCATTCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTg                               <  2:545338/151‑1 (MQ=255)
                                                                                                                    ttCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCg                            <  2:249092/151‑1 (MQ=255)
                                                                                                                     tCGGTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGt                           >  1:2989193/1‑151 (MQ=255)
                                                                                                                       ggTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGc                         >  1:2647043/1‑151 (MQ=255)
                                                                                                                        gTGACGTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCa                        >  1:5356601/1‑151 (MQ=255)
                                                                                                                             gTTCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcgga                   >  1:1957112/1‑151 (MQ=255)
                                                                                                                              ttCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggat                  >  1:1659723/1‑151 (MQ=255)
                                                                                                                              ttCTGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggat                  >  1:3787886/1‑151 (MQ=255)
                                                                                                                                 tGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggatcgg               >  2:4828469/1‑151 (MQ=255)
                                                                                                                                 tGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCatcggatcgg               >  2:129860/1‑151 (MQ=255)
                                                                                                                                 tGCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCa                        <  1:6336512/142‑1 (MQ=255)
                                                                                                                                 tGCCCAGTCATCTCGAAATGGATACCGCCGGCATTGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTCATGGTGTTGCCGTGCatcgggtcgg               >  2:100005/1‑151 (MQ=255)
                                                                                                                                  gCCCAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGCGCAccggatcgga              >  1:582051/1‑151 (MQ=255)
                                                                                                                                   cccagTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGaa             <  1:2152980/147‑1 (MQ=255)
                                                                                                                                   cccAGTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGaa             >  1:2286498/1‑151 (MQ=255)
                                                                                                                                     cagTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACt           <  2:2757056/149‑1 (MQ=255)
                                                                                                                                      agTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTc          <  2:5299190/150‑1 (MQ=255)
                                                                                                                                       gTCATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTcc         >  2:2330080/1‑151 (MQ=255)
                                                                                                                                         cATCTCGATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCcac       >  2:1228941/1‑151 (MQ=255)
                                                                                                                                            ctcgATATGGATACCGCCGGCGTGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGc    <  2:1836216/148‑1 (MQ=255)
                                                                                                                                            ctcgATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCTCAGc    <  2:5369109/148‑1 (MQ=255)
                                                                                                                                            ctcgATATGGATACCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGc    >  2:5977218/4‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGGCAATCAGGAATGCCAACTCCAGCGACTGATCGGCATTCATCCGCGGATCACAATGGGTGTGATACCGGTCAGACAACGCATCCTCGGTGATCGGCCGGGCACCGCCGATGCATTCGGTGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACTTGGAAGAACTGCTTCACCTCATCCAGAATCTGCGCAAAATCGCGGGTCTTGTAACCGCTGCTGGCCTTGATGGTGTTGCCGTGCATCGGATCGGAACTCCACAGCAC  >  CP009974/3756783‑3757075

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: