Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,105,405 T→A R149R (CGA→CGT RPPX_22650 ← benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,105,4050TA94.9% 474.8 / 23.9 241R149R (CGA→CGTRPPX_22650benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (6/6);  new base A (111/112);  total (117/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 241 total aligned reads displayed.

TACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTT  >  CP009974/5105256‑5105548
                                                                                                                                                     |                                                                                                                                               
taggAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATc                                                                                                                                                <  1:4872992/148‑1 (MQ=1)
          ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCtt                                                                                                                                      <  2:2173939/146‑1 (MQ=1)
                             accgtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCg                                                                                                                   >  2:1885131/8‑151 (MQ=0)
                                gtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCg                                                                                                                <  2:1303925/147‑1 (MQ=1)
                                         cTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAg                                                                                                       <  1:5874292/151‑1 (MQ=11)
                                         cTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAg                                                                                                       <  2:5293205/151‑1 (MQ=11)
                                          tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                                                      <  1:4520260/151‑1 (MQ=11)
                                                acgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgg                                                                                                <  1:1770899/151‑1 (MQ=11)
                                                 cgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCAACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCggg                                                                                               >  2:6656563/1‑151 (MQ=11)
                                                  gaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgggg                                                                                              >  2:5189395/1‑151 (MQ=11)
                                                   aTCATGTCGATGGTCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGt                                                                                             <  2:785388/151‑1 (MQ=11)
                                                   aTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGt                                                                                             <  2:906415/151‑1 (MQ=11)
                                                    tCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTa                                                                                            <  1:6003320/151‑1 (MQ=11)
                                                     cATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAg                                                                                           <  1:6995/151‑1 (MQ=11)
                                                      aTGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGc                                                                                          >  1:6128117/1‑151 (MQ=11)
                                                       tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                                                                         >  1:5788437/1‑151 (MQ=11)
                                                        gTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCg                                                                                        >  1:3939670/1‑151 (MQ=11)
                                                        gTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCg                                                                                        <  2:1923770/151‑1 (MQ=11)
                                                        gTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCg                                                                                        <  1:416015/151‑1 (MQ=11)
                                                         tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                                                       <  1:6337498/151‑1 (MQ=11)
                                                          cGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGCCTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCTGc                                                                                      >  2:5445050/1‑151 (MQ=11)
                                                          cGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGc                                                                                      >  1:1010144/1‑151 (MQ=11)
                                                          cGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGc                                                                                      >  1:5387322/1‑151 (MQ=11)
                                                          cGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGc                                                                                      <  1:2609441/151‑1 (MQ=11)
                                                            atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                                                    >  1:6692084/1‑151 (MQ=11)
                                                            atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                                                    >  1:4046096/1‑151 (MQ=11)
                                                             tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                                                   <  2:5900619/151‑1 (MQ=11)
                                                             tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                                                   >  2:6596830/1‑151 (MQ=11)
                                                               atCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGc                                                                                 <  1:6758249/151‑1 (MQ=11)
                                                                tCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCt                                                                                <  1:4905646/151‑1 (MQ=11)
                                                                 cTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTg                                                                               <  1:7047679/151‑1 (MQ=11)
                                                                     gCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTggggg                                                                           >  2:1707433/1‑150 (MQ=11)
                                                                     gCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                                                           <  2:1912434/151‑1 (MQ=11)
                                                                     gCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                                                           >  2:3133233/1‑151 (MQ=11)
                                                                      ccGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTc                                                                          >  1:973589/1‑151 (MQ=11)
                                                                      ccGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTc                                                                          >  2:1535677/1‑151 (MQ=11)
                                                                      ccGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTc                                                                          >  2:2226181/1‑151 (MQ=11)
                                                                       cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                                                                         >  2:236829/1‑151 (MQ=11)
                                                                       cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                                                                         <  2:1148721/151‑1 (MQ=11)
                                                                       cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                                                                         <  2:1586528/151‑1 (MQ=11)
                                                                        gACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCtt                                                                        <  1:6218094/151‑1 (MQ=11)
                                                                           tCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCa                                                                     >  2:2940468/1‑151 (MQ=11)
                                                                           tCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCa                                                                     >  2:1230745/1‑151 (MQ=11)
                                                                            cGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAc                                                                    <  1:6578372/151‑1 (MQ=11)
                                                                             gCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAcc                                                                   <  2:6965084/151‑1 (MQ=11)
                                                                              cccAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCt                                                                  >  1:2185702/1‑151 (MQ=11)
                                                                                cAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTc                                                                >  2:3869785/1‑151 (MQ=14)
                                                                                cAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTc                                                                <  1:369595/151‑1 (MQ=14)
                                                                                cAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTc                                                                >  2:6358574/1‑151 (MQ=11)
                                                                                cAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTc                                                                >  1:3232881/1‑151 (MQ=14)
                                                                                  ggTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGCTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCACTTCAAGCTGTCGGGGTAGCCGGCGTTGCCGGGGTCTTTCACCTTcac                                                              >  2:1079564/1‑150 (MQ=11)
                                                                                     gCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGACGGCGTTGCTGGGGTCTTTCACCTTcagcag                                                           <  1:1163881/151‑1 (MQ=14)
                                                                                       tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCCGCAGtt                                                         >  2:3676521/1‑151 (MQ=14)
                                                                                       tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt                                                         >  2:4568949/1‑151 (MQ=14)
                                                                                       tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt                                                         <  2:3781455/151‑1 (MQ=14)
                                                                                       tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt                                                         <  2:3488041/151‑1 (MQ=14)
                                                                                       tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGACCCGTCGCAGTTGAAGCTGTCCGGGTAGCCGGCGTTCCTGCGGTCTTTCACCTCCAGCACtt                                                         >  2:2742327/1‑151 (MQ=14)
                                                                                        cGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGttt                                                        <  2:7076796/151‑1 (MQ=14)
                                                                                        cGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCCGCGTTGCTGGGGTCTTTCACCTTCAGCAGttt                                                        <  2:2548841/151‑1 (MQ=14)
                                                                                         gACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTg                                                       >  1:2438276/1‑151 (MQ=14)
                                                                                         gACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGGCCTTCACCTTCTGCAGGTTg                                                       >  2:7061621/1‑151 (MQ=14)
                                                                                          aCCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGc                                                      <  1:157218/151‑1 (MQ=14)
                                                                                           ccAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGcc                                                     <  2:1011087/151‑1 (MQ=14)
                                                                                            cagcggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCg                                                    >  2:6934943/8‑151 (MQ=14)
                                                                                              gcggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCt                                                  >  1:3774529/6‑151 (MQ=14)
                                                                                              gcggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGAGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCt                                                  <  1:995709/146‑1 (MQ=14)
                                                                                               cggcttCGCATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtg                                                 <  2:4460417/147‑1 (MQ=14)
                                                                                               cggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtg                                                 >  1:3617762/5‑151 (MQ=14)
                                                                                                 gcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgtt                                               >  1:5336798/3‑151 (MQ=14)
                                                                                                    tCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTgttgtt                                            >  1:3056098/1‑151 (MQ=14)
                                                                                                      acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa                                          >  1:3680399/1‑151 (MQ=14)
                                                                                                      acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa                                          <  2:3170272/151‑1 (MQ=14)
                                                                                                       caTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGaa                                         >  1:2453829/1‑151 (MQ=14)
                                                                                                        aTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAAc                                        <  1:6096147/151‑1 (MQ=14)
                                                                                                         tCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACg                                       <  2:5214996/151‑1 (MQ=14)
                                                                                                          cGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGt                                      >  2:1408895/1‑151 (MQ=14)
                                                                                                             cgttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCa                                   <  1:1503754/149‑1 (MQ=14)
                                                                                                             cgttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCa                                   <  2:3088293/149‑1 (MQ=14)
                                                                                                              gttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAg                                  >  1:4178297/2‑151 (MQ=14)
                                                                                                                tgaggcTGCCAAACAAAAAGCCCCGGTACGACTCAAGCCGGGCCACCTGGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGACGCTGTTGTTGAACGTCCAGcc                                <  2:5333612/146‑1 (MQ=14)
                                                                                                                 caggcTGCCGAACAAAGAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCg                               <  2:6813483/147‑1 (MQ=14)
                                                                                                                    gcTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTgg                            <  2:4169947/150‑1 (MQ=14)
                                                                                                                     cTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGa                           <  1:5141018/151‑1 (MQ=14)
                                                                                                                          gaacaaaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACggg                      <  1:1988425/144‑1 (MQ=255)
                                                                                                                           aacaaaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGc                     >  2:2405631/7‑151 (MQ=255)
                                                                                                                           aacaaaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGCCCAGCCGTGGAACGGGc                     >  2:5794531/7‑151 (MQ=14)
                                                                                                                            acaaaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCa                    <  1:3500387/146‑1 (MQ=255)
                                                                                                                            acaaaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCa                    >  2:5871760/6‑151 (MQ=255)
                                                                                                                              aaaacgCCCCGGTACGACTCAACACGCGCCACCTTGGTCACGTCGTGGGGGCCGTCGCACTTGAAGCTGCCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTCTTGCTGAGCGTCCAGCCGTGGAACGGGCAgg                  <  2:5556246/146‑1 (MQ=255)
                                                                                                                               aaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGt                 <  1:4621241/149‑1 (MQ=255)
                                                                                                                               aaaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGt                 <  1:7105614/149‑1 (MQ=255)
                                                                                                                                aaaGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGta                >  1:4242050/2‑151 (MQ=255)
                                                                                                                                    ccccGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGc            <  1:6859951/151‑1 (MQ=255)
                                                                                                                                    ccccGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGc            <  2:2531536/151‑1 (MQ=255)
                                                                                                                                    ccccGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGc            <  1:4090435/151‑1 (MQ=255)
                                                                                                                                     cccGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCt           >  2:2234149/1‑151 (MQ=255)
                                                                                                                                      ccGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTg          >  1:6512210/1‑151 (MQ=255)
                                                                                                                                           aCGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGt                                      <  2:2441544/118‑1 (MQ=255)
                                                                                                                                             gACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGt   >  1:2135488/1‑151 (MQ=255)
                                                                                                                                              aCTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGtt  >  1:3327825/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGAACGGGCAGGTATAGCTGGAACGGTT  >  CP009974/5105256‑5105548

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: