Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,105,354 T→C K166K (AAA→AAG)  RPPX_22650 ← benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,105,3540TC87.0% 142.2 / 28.2 101K166K (AAA→AAG) RPPX_22650benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (6/7);  new base C (39/48);  total (46/55)
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.73e-01

Only 100 of 101 total aligned reads displayed.

TTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGC  >  CP009974/5105220‑5105496
                                                                                                                                      |                                                                                                                                              
ttctCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                                                                                                                                >  2:668889/4‑151 (MQ=1)
ttctCGGCGGTCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                                                                                                                                <  2:3818708/148‑1 (MQ=1)
ttctCGGCGGTCAGCTTCCAGTCACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTTGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                                                                                                                                <  2:5234826/148‑1 (MQ=1)
         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                                                                                                                       >  2:2566899/4‑151 (MQ=1)
         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                                                                                                                       <  1:2997544/148‑1 (MQ=1)
                    gTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc                                                                                                            <  2:2447120/151‑1 (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)
                                               tccacgcAGCACTTCCAGGCGGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTg                                                                                 <  2:3053862/147‑1 (MQ=1)
                                               tccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTg                                                                                 >  1:3908769/5‑151 (MQ=1)
                                                ccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTgg                                                                                >  1:3929407/4‑151 (MQ=1)
                                                       ggacgtCCGGGCCTCCAGGCGACTGATCGACGATCATGGCGATGACCTTCGCCGCCCCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTAAGGATGAAGAAAAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTc                                                                         <  2:2219530/146‑1 (MQ=0)
                                                       gcgctgccAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTc                                                                         <  2:3610142/145‑1 (MQ=1)
                                                                    gtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCg                                                            <  2:3314535/147‑1 (MQ=1)
                                                                     tccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGGTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGACTCGTGCCACCTTGGTCAGGGCGTGGGAGCCGTCGc                                                           >  2:5211982/4‑151 (MQ=1)
                                                                              tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                  <  2:1751245/151‑1 (MQ=11)
                                                                              tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                  <  1:3954260/151‑1 (MQ=11)
                                                                              tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                  <  1:4700014/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)
                                                                                  cgacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTc                                              >  2:3141291/1‑151 (MQ=11)
                                                                                   gacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCg                                             <  1:641181/151‑1 (MQ=11)
                                                                                   gacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCg                                             <  1:5247457/151‑1 (MQ=11)
                                                                                   gacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCg                                             <  1:3126539/151‑1 (MQ=11)
                                                                                   gacgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCg                                             >  1:3550411/1‑151 (MQ=11)
                                                                                    acgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgg                                            >  1:2776416/1‑151 (MQ=11)
                                                                                    acgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCACGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgg                                            >  2:5190524/1‑151 (MQ=11)
                                                                                     cgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCggg                                           <  1:128520/151‑1 (MQ=11)
                                                                                      gaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgtgg                                          >  1:410258/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                                         <  2:1869471/151‑1 (MQ=11)
                                                                                       aTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGt                                         >  1:4082978/1‑151 (MQ=11)
                                                                                           tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                     <  1:4528427/151‑1 (MQ=11)
                                                                                           tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                     <  1:4623695/151‑1 (MQ=11)
                                                                                           tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                     <  1:2923602/151‑1 (MQ=11)
                                                                                            gTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCg                                    <  1:4637503/151‑1 (MQ=11)
                                                                                             tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                   <  1:2537538/151‑1 (MQ=11)
                                                                                             tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                   >  2:1990908/1‑151 (MQ=11)
                                                                                             tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                   <  1:495812/151‑1 (MQ=11)
                                                                                             tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                   >  1:289304/1‑151 (MQ=11)
                                                                                             tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCAgg                                   >  1:287132/1‑151 (MQ=11)
                                                                                               gatgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCg                                 >  1:4295127/1‑151 (MQ=11)
                                                                                                 tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                               <  1:2795052/151‑1 (MQ=11)
                                                                                                 tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                               <  1:3400664/151‑1 (MQ=11)
                                                                                                 tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                               >  1:3592006/1‑151 (MQ=11)
                                                                                                   atCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGc                             <  1:4534713/151‑1 (MQ=11)
                                                                                                    tCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCt                            <  2:3101252/151‑1 (MQ=11)
                                                                                                     cTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTg                           >  2:2594966/1‑151 (MQ=11)
                                                                                                     cTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTg                           <  1:2678590/151‑1 (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:511626/1‑151 (MQ=11)
                                                                                                           cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                     <  1:1682256/151‑1 (MQ=11)
                                                                                                           cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                     >  1:2684962/1‑151 (MQ=11)
                                                                                                           cGACTCGCCCAGCTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                     <  1:3817212/151‑1 (MQ=11)
                                                                                                            gACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCTAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCtt                    <  2:5223292/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)
                                                                                                             aCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCACGTCGTGGGAGCCGACGCAGTTGAAGCAGACGGGGTAGCCGACGTTGCTAGGGACtat                   >  2:610324/1‑149 (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               <  1:2829674/151‑1 (MQ=11)
                                                                                                                 gCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAcc               >  1:5245481/1‑151 (MQ=14)
                                                                                                                 gCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCAcc               <  2:2013951/151‑1 (MQ=11)
                                                                                                                  cccAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCt              >  1:1023682/1‑151 (MQ=14)
                                                                                                                   ccAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCtt             >  1:2304386/1‑151 (MQ=11)
                                                                                                                    cAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTc            <  2:4049353/151‑1 (MQ=14)
                                                                                                                    cAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTc            >  2:2628357/1‑151 (MQ=11)
                                                                                                                     aGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTca           <  2:3955464/151‑1 (MQ=14)
                                                                                                                      ggTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcag          >  2:2115500/1‑151 (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     >  2:4129790/1‑151 (MQ=14)
                                                                                                                           tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt     >  2:530034/1‑151 (MQ=14)
                                                                                                                           tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt     <  1:2821682/151‑1 (MQ=14)
                                                                                                                            cGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGttt    <  2:4516520/151‑1 (MQ=14)
                                                                                                                             gACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTg   >  2:3328693/1‑151 (MQ=14)
                                                                                                                             gACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTg   <  1:588844/151‑1 (MQ=14)
                                                                                                                              aCCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGc  <  1:3215567/151‑1 (MQ=14)
                                                                                                                                      |                                                                                                                                              
TTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGC  >  CP009974/5105220‑5105496

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: ‑