Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 5,105,366 A→G N162N (AAT→AAC)  RPPX_22650 ← benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099745,105,3660AG85.0% 194.5 / 49.6 133N162N (AAT→AAC) RPPX_22650benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (12/8);  new base G (52/61);  total (64/69)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 133 total aligned reads displayed.

TTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTC  >  CP009974/5105220‑5105513
                                                                                                                                                  |                                                                                                                                                   
ttctCGGCGATCAGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCa                                                                                                                                                 >  2:2585471/4‑151 (MQ=1)
         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                                                                                                                                        >  2:1220759/4‑151 (MQ=1)
         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                                                                                                                                        >  1:150966/4‑151 (MQ=1)
         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                                                                                                                                        >  1:3043802/4‑151 (MQ=1)
         gtcaGCTTCCAGTTACCTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCaa                                                                                                                                        >  1:3043435/4‑151 (MQ=1)
                    gTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGccc                                                                                                                             <  1:884125/151‑1 (MQ=1)
                     ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc                                                                                                                            <  1:1110405/148‑1 (MQ=1)
                     ttaccTTCGTAGATGTAGGAGCTAGCTCCACGCAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGcccc                                                                                                                            <  1:3706044/148‑1 (MQ=1)
                         ctgcGTAGGTGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCAGAACAAAAAGCCCCGGTa                                                                                                                        <  2:2022719/148‑1 (MQ=1)
                                            agctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCAcc                                                                                                     <  2:3887493/144‑1 (MQ=1)
                                             gctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGACGCCCCGGTACGACTCGAATCGTGCCACCt                                                                                                    <  1:1196736/145‑1 (MQ=1)
                                              ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCtt                                                                                                   >  1:3533153/6‑151 (MQ=1)
                                              ctccacgcAGCACTTCCAGACCGTCCGGCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCtt                                                                                                   <  1:3582024/146‑1 (MQ=1)
                                                                  ccgtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGTTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGGGCCGt                                                                               >  1:2209783/7‑151 (MQ=1)
                                                                   cgtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTc                                                                              >  1:4543423/6‑151 (MQ=0)
                                                                   cgtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTc                                                                              >  2:3315559/6‑151 (MQ=0)
                                                                    gtccggCGACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGCGCTCGACCAGTGGTTTCACATCGGCAGTCAGACTGCCCAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTCGGTCAGGCCCTGGGAGCCCTCg                                                                             >  2:4275632/5‑151 (MQ=1)
                                                                       gggggACTGGTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAg                                                                          <  2:3613312/147‑1 (MQ=0)
                                                                              tGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGc                                                                   >  2:445610/1‑151 (MQ=11)
                                                                               gATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt                                                                  >  2:3963779/1‑151 (MQ=11)
                                                                                aTCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTg                                                                 >  1:1745460/1‑151 (MQ=11)
                                                                                  cggcgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTc                                                               <  1:2794826/148‑1 (MQ=11)
                                                                                     cgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt                                                                  <  2:610701/145‑1 (MQ=11)
                                                                                     cgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCt                                                                  >  1:610701/1‑145 (MQ=11)
                                                                                     cgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCggg                                                            >  1:2116026/1‑151 (MQ=11)
                                                                                     cgaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCggg                                                            <  2:3062949/151‑1 (MQ=11)
                                                                                      gaTCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCgggg                                                           >  1:4409847/1‑151 (MQ=11)
                                                                                        tCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTa                                                         <  1:1779292/151‑1 (MQ=11)
                                                                                        tCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTa                                                         <  1:1057282/151‑1 (MQ=11)
                                                                                        tCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTa                                                         >  2:3771655/1‑151 (MQ=11)
                                                                                        tCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTa                                                         <  1:2460073/151‑1 (MQ=11)
                                                                                         cATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAg                                                        >  1:1465485/1‑151 (MQ=11)
                                                                                         cATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAg                                                        >  2:1802813/1‑151 (MQ=11)
                                                                                           tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                                      >  1:3363511/1‑151 (MQ=11)
                                                                                           tGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGcc                                                      <  2:1504569/151‑1 (MQ=11)
                                                                                            gTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCg                                                     <  1:4467167/151‑1 (MQ=11)
                                                                                            gTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCg                                                     >  2:4378176/1‑151 (MQ=11)
                                                                                             tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                    <  2:2343252/151‑1 (MQ=11)
                                                                                             tCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCCTGGTCAGGTCGCGGGAGCCGTCGCAGTTCAAGCTGTCGGGGTAGCCgg                                                    >  2:4227764/1‑151 (MQ=11)
                                                                                             tCGATGATCCTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCgg                                                    <  2:4231015/151‑1 (MQ=11)
                                                                                                atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                 >  2:3041955/1‑151 (MQ=11)
                                                                                                atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                 >  2:922196/1‑151 (MQ=11)
                                                                                                atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                 >  1:345807/1‑151 (MQ=11)
                                                                                                atgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGt                                                 >  1:2468429/1‑151 (MQ=11)
                                                                                                 tgatCTTCGCCGCCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                <  1:4716245/151‑1 (MQ=11)
                                                                                                 tgatCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                <  1:162535/151‑1 (MQ=11)
                                                                                                 tgatCTTCGCCGACTCGCCCAGCTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGtt                                                >  2:4654113/1‑151 (MQ=11)
                                                                                                   atCTTCGCCGGCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGc                                              <  2:1824649/151‑1 (MQ=11)
                                                                                                   atCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGc                                              >  2:2289882/1‑151 (MQ=11)
                                                                                                    tCTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCt                                             >  2:4641988/1‑151 (MQ=11)
                                                                                                     cTTCGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTg                                            <  1:2941509/151‑1 (MQ=11)
                                                                                                        ggCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTgggg                                         <  2:1005725/150‑1 (MQ=11)
                                                                                                        cGCCGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTgggg                                         >  2:1180023/1‑151 (MQ=11)
                                                                                                         gCCGGCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                        <  2:2213650/151‑1 (MQ=11)
                                                                                                         gCCGAGTGGCGCGGGTGCGCGACCGGCGGCTTGACATCGGCGTTCAGGCGGCCGAAGAAAAAGCCCCGGTACGACTCAAAACGGGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                        <  2:4385034/151‑1 (MQ=2)
                                                                                                         gCCGACTCGCCCAGGTGCGCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGt                                        <  2:1596594/151‑1 (MQ=11)
                                                                                                          ccGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTc                                       >  2:1379919/1‑151 (MQ=11)
                                                                                                           cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGCCCCGTCGCCGTTGAACCTGTCCGCGTCGCCGGCGTTGCTGGGGGCt                                      >  2:2743638/1‑151 (MQ=11)
                                                                                                           cGACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCt                                      >  2:5089971/1‑151 (MQ=11)
                                                                                                            gACTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCtt                                     <  1:313509/151‑1 (MQ=11)
                                                                                                             ccTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCttt                                    <  2:5070714/150‑1 (MQ=11)
                                                                                                             aCTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCttt                                    >  1:1775858/1‑151 (MQ=11)
                                                                                                              cTCGCCCAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTc                                   >  1:4064556/1‑151 (MQ=11)
                                                                                                                  cccAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCt                               <  1:3104213/151‑1 (MQ=11)
                                                                                                                   ccAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCtt                              >  2:1604955/1‑151 (MQ=14)
                                                                                                                   ccAGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCtt                              <  1:4976132/151‑1 (MQ=14)
                                                                                                                     gggTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTca                            <  2:194868/150‑1 (MQ=14)
                                                                                                                     aGGTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTca                            <  1:1275797/151‑1 (MQ=14)
                                                                                                                      ggTGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcag                           <  2:1578873/151‑1 (MQ=14)
                                                                                                                        tGCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcagca                         >  2:3927121/1‑151 (MQ=14)
                                                                                                                         gCTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTcagcag                        >  1:4013846/1‑151 (MQ=14)
                                                                                                                          cTCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGt                       >  1:2374371/1‑151 (MQ=14)
                                                                                                                           tCGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGtt                      >  2:1380660/1‑151 (MQ=14)
                                                                                                                            cGACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGGGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCCGttt                     >  1:4644516/1‑151 (MQ=14)
                                                                                                                             gACCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTg                    >  1:3102984/1‑151 (MQ=14)
                                                                                                                              aCCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGc                   <  1:4897466/151‑1 (MQ=14)
                                                                                                                              aCCAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGc                   <  2:1922833/151‑1 (MQ=14)
                                                                                                                               ccAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGcc                  <  2:5075057/151‑1 (MQ=14)
                                                                                                                               ccAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGcc                  <  2:4544586/151‑1 (MQ=14)
                                                                                                                                cAGCGGCTTCGCATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCg                 <  2:3483636/151‑1 (MQ=14)
                                                                                                                                cAGCGGCTTCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCg                 <  1:3401887/151‑1 (MQ=14)
                                                                                                                                 agcggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTCGGTCAGGTCGTGGGAGCCGTCGCACTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGc                >  2:661255/7‑151 (MQ=14)
                                                                                                                                  gcggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCt               <  2:1795694/146‑1 (MQ=14)
                                                                                                                                   cggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtg              <  1:3024940/147‑1 (MQ=14)
                                                                                                                                    ggcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgt             >  1:4850394/4‑151 (MQ=14)
                                                                                                                                     gcttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgtt            <  1:4842708/149‑1 (MQ=14)
                                                                                                                                      cttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgttg           <  2:2367358/150‑1 (MQ=14)
                                                                                                                                       ttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgttgt          <  2:1234176/151‑1 (MQ=14)
                                                                                                                                       ttCACATCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCtgttgt          >  1:3629493/1‑151 (MQ=14)
                                                                                                                                         cacaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGttgttg        >  1:5202230/1‑151 (MQ=14)
                                                                                                                                          acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa       >  2:379036/1‑151 (MQ=14)
                                                                                                                                          acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa       <  1:2686384/151‑1 (MQ=14)
                                                                                                                                          acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa       >  1:4608444/1‑151 (MQ=14)
                                                                                                                                          acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGa       >  1:4830295/1‑151 (MQ=14)
                                                                                                                                          acaTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTGGTTGa       >  1:4579950/1‑151 (MQ=14)
                                                                                                                                           caTCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGaa      <  1:665036/151‑1 (MQ=14)
                                                                                                                                            aTCGGCGTTCAGGGTGCCGAACGGACGGTCCGGGTACGACTCCACACGTGCCACTTCCGTCAGGTCGTGGGAGTCGTCGCAGTTGGAGCTGGCGGGGGAGGCGGCGTTGCTGGAGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGagc     <  2:3469906/151‑3 (MQ=255)
                                                                                                                                             tCGGCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACg    <  1:840488/151‑1 (MQ=14)
                                                                                                                                              gggCGTTCAGGCTGCCGAGCAGAAAGCCCCGGGACGACTCAGAACGTGCCACGTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAGGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTGCGCCTTCAGCAGTTTGCCGCTGTTGTTGAACGt   <  2:4163975/150‑1 (MQ=14)
                                                                                                                                               ggCGTTCAGGCTGCCGAACAAAAAGCCCCGGTACGACTCAAAACGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTc  <  1:2287929/151‑1 (MQ=14)
                                                                                                                                                  |                                                                                                                                                   
TTTTCGGCGGTGAGCTTCCAGTTGCCTTCGTAGATGTACGAGCTGGCACCGCGCAGCACTTCCAGGCCTTCAGGCGACTGATCGACGATCATGTCGATGATCTTCGCCGACTCGCCCAGGTGCTCGACCAGTGGTTTCACATCGGCATTCAGACTGCCAAACAGGAAGCCCCGGTACGACTCGAATCGTGCCACCTTGGTCAGGTCGTGGGAGCCGTCGCAGTTGAAGCTGTCGGGGTAGCCGGCGTTGCTGGGGTCTTTCACCTTCAGCAGTTTGCCGCTGTTGTTGAACGTC  >  CP009974/5105220‑5105513

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