Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,105,378 A→G F158F (TTT→TTC)  RPPX_22650 ← benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,105,3780AG93.2% 220.0 / 13.8 134F158F (TTT→TTC) RPPX_22650benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base G (54/70);  total (59/75)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

GTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGA  >  CP009974/5105229‑5105523
                                                                                                                                                     |                                                                                                                                                 
gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                                                                                                                                                  >  2:2566899/4‑151 (MQ=1)
           gTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc                                                                                                                                       <  2:2449157/151‑1 (MQ=1)
                           taggAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATc                                                                                                                       <  2:1263687/148‑1 (MQ=1)
                                     ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTTCCACCtt                                                                                                             >  2:1941651/6‑151 (MQ=1)
                                     ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCtt                                                                                                             >  1:2652791/6‑151 (MQ=1)
                                      tccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTg                                                                                                            >  1:3908769/5‑151 (MQ=1)
                                              ggacgtCCGGGCCTCCAGGCGACTGATCGACGATCATGGCGATGACCTTCGCCGCCCCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTAAGGATGAAGAAAAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTc                                                                                                    <  2:2219530/146‑1 (MQ=0)
                                                            tccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGGTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGACTCGTGCCACCTTGGTCAGGGCGTGGGAGCCGTCGc                                                                                      >  2:5211982/4‑151 (MQ=1)
                                                                     tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                             <  1:4700014/151‑1 (MQ=11)
                                                                     tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                             <  2:1751245/151‑1 (MQ=11)
                                                                     tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                             <  1:3954260/151‑1 (MQ=11)
                                                                     tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                             >  2:3165189/1‑151 (MQ=11)
                                                                      gATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt                                                                            >  1:2442686/1‑151 (MQ=11)
                                                                       aTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTg                                                                           <  1:812696/151‑1 (MQ=11)
                                                                        tCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGt                                                                          <  1:817078/151‑1 (MQ=11)
                                                                          gacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCg                                                                        <  1:5247457/151‑1 (MQ=11)
                                                                          gacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCg                                                                        >  1:3550411/1‑151 (MQ=11)
                                                                          gacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCg                                                                        <  1:3126539/151‑1 (MQ=11)
                                                                           acgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgg                                                                       >  1:2776416/1‑151 (MQ=11)
                                                                           acgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCACGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgg                                                                       >  2:5190524/1‑151 (MQ=11)
                                                                             gaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgggg                                                                     >  1:5192008/1‑151 (MQ=11)
                                                                             gaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGGACGACTCACAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgggg                                                                     >  2:3823158/1‑151 (MQ=11)
                                                                              aTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGt                                                                    <  1:1099162/151‑1 (MQ=11)
                                                                              aTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGt                                                                    >  1:4082978/1‑151 (MQ=11)
                                                                                  tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                                                <  1:4528427/151‑1 (MQ=11)
                                                                                  tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                                                <  1:2923602/151‑1 (MQ=11)
                                                                                   gTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCg                                                               <  1:4637503/151‑1 (MQ=11)
                                                                                    tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                              >  1:289304/1‑151 (MQ=11)
                                                                                    tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                              >  2:1990908/1‑151 (MQ=11)
                                                                                    tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                              <  1:2537538/151‑1 (MQ=11)
                                                                                    tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                              <  1:495812/151‑1 (MQ=11)
                                                                                    tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCAgg                                                              >  1:287132/1‑151 (MQ=11)
                                                                                      gatgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCg                                                            >  1:4295127/1‑151 (MQ=11)
                                                                                        tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                          <  1:3400664/151‑1 (MQ=11)
                                                                                        tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                          <  1:2795052/151‑1 (MQ=11)
                                                                                          atCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGc                                                        <  1:4534713/151‑1 (MQ=11)
                                                                                           tCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCt                                                       <  2:3101252/151‑1 (MQ=11)
                                                                                            cTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTg                                                      <  1:2678590/151‑1 (MQ=11)
                                                                                            cTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTg                                                      >  2:2594966/1‑151 (MQ=11)
                                                                                             ttCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTgg                                                     <  1:1785677/151‑1 (MQ=11)
                                                                                              tCGCCGACTCGCCCAGGTGCTCGACCAGTGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTggg                                                    >  1:5424821/1‑151 (MQ=11)
                                                                                              tCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTggg                                                    >  1:1078796/1‑151 (MQ=11)
                                                                                              tCGCCGACGCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTggg                                                    <  1:240061/151‑1 (MQ=11)
                                                                                               cGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAGAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTgggg                                                   <  2:3295461/151‑1 (MQ=11)
                                                                                               cGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGTGCCGTCGCAGTTGAAGCTGTCCGGGTAGCCGGCGTTGCTgggg                                                   >  2:1514115/1‑151 (MQ=11)
                                                                                                gCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                                  >  1:4866563/1‑151 (MQ=11)
                                                                                                  cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                                                >  1:2684962/1‑151 (MQ=11)
                                                                                                  cGACTCGCCCAGCTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                                                <  1:3817212/151‑1 (MQ=11)
                                                                                                    aCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                             >  2:3117428/1‑120 (MQ=11)
                                                                                                    aCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                             <  1:3117428/120‑1 (MQ=11)
                                                                                                    aCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCttt                                              >  2:610384/1‑151 (MQ=11)
                                                                                                    aCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCttt                                              >  1:75877/1‑151 (MQ=11)
                                                                                                     cTCGCCCAGGTGCTCGCCCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTc                                             <  2:483519/151‑1 (MQ=11)
                                                                                                     cTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTc                                             >  2:3616139/1‑151 (MQ=11)
                                                                                                      ccGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCGCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCa                                            <  2:5177899/150‑1 (MQ=11)
                                                                                                       cGCCCGGGTGCTCGCCCAGCGGCTTCACATCGGCGTTCCGGCTGCCGAACAAAAAGCCTCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAc                                           <  1:4809361/151‑1 (MQ=11)
                                                                                                       cGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAc                                           >  2:356838/1‑151 (MQ=11)
                                                                                                       cGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAc                                           <  2:367888/151‑1 (MQ=11)
                                                                                                        gCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAcc                                          <  1:2846409/151‑1 (MQ=11)
                                                                                                        gCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAcc                                          <  2:2013951/151‑1 (MQ=11)
                                                                                                        gCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAcc                                          >  1:5245481/1‑151 (MQ=14)
                                                                                                         cccAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCt                                         >  1:1023682/1‑151 (MQ=14)
                                                                                                           cAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTc                                       <  2:4049353/151‑1 (MQ=14)
                                                                                                             ggTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcag                                     <  2:2157034/151‑1 (MQ=14)
                                                                                                             ggTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcag                                     >  1:4574808/1‑151 (MQ=14)
                                                                                                             ggTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGCTGCTGGGGTCCTTCACCTCcag                                     >  2:5295321/1‑151 (MQ=11)
                                                                                                                gCTCGACCAGCAGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcagcag                                  <  1:1671327/151‑1 (MQ=14)
                                                                                                                 cTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGt                                 <  2:1078796/151‑1 (MQ=14)
                                                                                                                  tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt                                <  1:2821682/151‑1 (MQ=14)
                                                                                                                  tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt                                >  2:4129790/1‑151 (MQ=14)
                                                                                                                  tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt                                >  2:530034/1‑151 (MQ=14)
                                                                                                                    gACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTg                              <  1:588844/151‑1 (MQ=14)
                                                                                                                    gACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTg                              >  2:3328693/1‑151 (MQ=14)
                                                                                                                     aCCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGc                             <  1:3215567/151‑1 (MQ=14)
                                                                                                                           ggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgt                       <  1:3040126/148‑1 (MQ=14)
                                                                                                                           ggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgt                       >  2:3036942/4‑151 (MQ=14)
                                                                                                                           ggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGAtgt                       >  2:5071489/4‑151 (MQ=14)
                                                                                                                             cttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgttg                     <  1:4861306/150‑1 (MQ=14)
                                                                                                                              ttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgttgt                    <  2:847305/151‑1 (MQ=14)
                                                                                                                               tCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTgttgtt                   <  2:4653999/151‑1 (MQ=14)
                                                                                                                                cacaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGttgttg                  <  2:2782115/151‑1 (MQ=14)
                                                                                                                                cacaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGttgttg                  <  1:1384509/151‑1 (MQ=14)
                                                                                                                                 ccaTGGGCGTTCAGGCTGCCGAACAAAGAGCCCGGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa                 <  2:5388972/150‑1 (MQ=14)
                                                                                                                                   aTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCCCCTTCAGCAGTTTGCCGCTGTTGTTGAAc               >  2:244633/1‑151 (MQ=14)
                                                                                                                                   aTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAAc               <  1:2946515/151‑1 (MQ=14)
                                                                                                                                       gcgttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGtgc           >  1:5047285/4‑149 (MQ=14)
                                                                                                                                       gcgttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTcc           >  1:1972256/4‑151 (MQ=14)
                                                                                                                                        cgttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCa          >  1:132092/3‑151 (MQ=14)
                                                                                                                                         gttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGATCGTcgtg         >  2:3337609/2‑148 (MQ=14)
                                                                                                                                         gttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAg         <  1:1926972/150‑1 (MQ=14)
                                                                                                                                          ttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGc        <  2:324767/151‑1 (MQ=14)
                                                                                                                                          ttCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGc        >  1:2577512/1‑151 (MQ=14)
                                                                                                                                           tCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGcc       <  1:4170435/151‑1 (MQ=14)
                                                                                                                                            caggcTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCg      <  1:3961918/147‑1 (MQ=14)
                                                                                                                                             aggcTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGt     <  1:5041301/148‑1 (MQ=14)
                                                                                                                                              ggcTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTg    <  2:1288994/149‑1 (MQ=14)
                                                                                                                                               gcTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTgg   <  1:3605870/150‑1 (MQ=14)
                                                                                                                                               gcTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTgg   <  1:2422155/150‑1 (MQ=14)
                                                                                                                                                cTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGACTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGa  >  1:5394675/1‑151 (MQ=14)
                                                                                                                                                cTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCCTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGCTCTTTCACCTCCCGCAGTTTGCCACTGTTGTTGAACGTCCAGCCGtcgc  >  2:3116007/1‑148 (MQ=14)
                                                                                                                                                     |                                                                                                                                                 
GTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTCCAGCCGTGGA  >  CP009974/5105229‑5105523

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