Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 2,738,522 T→C S245S (AGT→AGC RPPX_12165 → benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099742,738,5220TC99.5% 523.3 / NA 214S245S (AGT→AGCRPPX_12165benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (112/101);  minor base A (0/1);  total (112/102)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

TCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGCGCA  >  CP009974/2738379‑2738669
                                                                                                                                               |                                                                                                                                                   
tcctACATCTACGAAGGTAACTGGACGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCCCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGgct                                                                                                                                              <  2:1706468/148‑2 (MQ=0)
                                    gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                                                                          >  1:3622923/4‑151 (MQ=0)
                                    gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                                                                          <  2:1672564/148‑1 (MQ=1)
                                         cGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAt                                                                                                     >  2:3785492/1‑151 (MQ=2)
                                          ggCGCTGACGGCTATCTTGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACTAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGTGCTAAGCAGGGTGGTGGTTTCTATTTCTTCTACCAtt                                                                                                    >  2:1718050/1‑150 (MQ=2)
                                          ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                                                                                    >  2:1419084/1‑151 (MQ=2)
                                          ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                                                                                    >  2:4213504/1‑151 (MQ=11)
                                          ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGGTTCTATTCCTTCGACCATg                                                                                                    >  1:2965556/1‑151 (MQ=11)
                                              ccGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGaca                                                                                                >  1:2410558/1‑151 (MQ=2)
                                               cGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGCTTCTATTCCTTCGACCATGGacac                                                                                               >  2:503022/1‑151 (MQ=2)
                                                gACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACAcc                                                                                              >  2:2642154/1‑151 (MQ=11)
                                                 aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                             <  1:5771364/151‑1 (MQ=2)
                                                 aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                             <  1:1718050/151‑1 (MQ=11)
                                                 aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                             <  2:1532646/151‑1 (MQ=2)
                                                  cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                                            <  2:580292/151‑1 (MQ=2)
                                                        tcatgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctc                                                                                      <  2:7240359/144‑1 (MQ=2)
                                                         catgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctct                                                                                     <  2:2207156/145‑1 (MQ=11)
                                                          atgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTg                                                                                    >  2:2936361/6‑151 (MQ=11)
                                                          atgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCCCTg                                                                                    >  2:5110652/6‑151 (MQ=11)
                                                           tgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTgg                                                                                   <  1:3399385/147‑1 (MQ=11)
                                                              taGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACt                                                                                >  2:4928561/2‑151 (MQ=14)
                                                              taGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTAGACt                                                                                <  2:1624006/150‑1 (MQ=11)
                                                                gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                                                                              >  1:7056799/1‑151 (MQ=14)
                                                                gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                                                                              >  1:1813370/1‑151 (MQ=14)
                                                                gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCAGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                                                                              <  1:3357047/151‑1 (MQ=14)
                                                                   ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                                                                           >  2:3283742/1‑151 (MQ=14)
                                                                   ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                                                                           >  1:5644927/1‑151 (MQ=14)
                                                                     gTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTggg                                                                         <  2:2694895/151‑1 (MQ=14)
                                                                     gTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTggg                                                                         >  2:5405172/1‑151 (MQ=14)
                                                                     gTGCACAGGCACTACGCCGCCACGCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGCACACCTGCTGCTCTGGACTCGGTggg                                                                         <  1:983503/151‑1 (MQ=14)
                                                                       gCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCg                                                                       <  2:911594/151‑1 (MQ=14)
                                                                         aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCGAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa                                                                     >  2:735129/1‑151 (MQ=14)
                                                                            ggAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGGCGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAccc                                                                  <  2:4752381/151‑1 (MQ=14)
                                                                            ggAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCACTGGACTCGGTGGCCGCAccc                                                                  >  2:6276987/1‑151 (MQ=14)
                                                                             gAACGACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc                                                                 <  1:5971270/151‑1 (MQ=14)
                                                                               aCTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTCAGCAGGGTGGAGGTTTCGATTCCTTAGACCATGGACACCTGCTGCTCTGGGCGCGGTGGGCGAAccccaa                                                               >  2:4338100/1‑149 (MQ=14)
                                                                                cTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAg                                                              <  2:2439463/151‑1 (MQ=17)
                                                                                 tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                             <  2:471/151‑1 (MQ=17)
                                                                                 tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                             >  2:5456609/1‑151 (MQ=17)
                                                                                 tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                             <  2:6376205/151‑1 (MQ=17)
                                                                                 tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                             <  2:1073/151‑1 (MQ=17)
                                                                                 tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                             <  1:7007762/151‑1 (MQ=17)
                                                                                 tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                             >  1:6574006/1‑151 (MQ=17)
                                                                                 tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGCGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                             >  1:6705862/1‑151 (MQ=17)
                                                                                  aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGa                                                            >  1:1521046/1‑151 (MQ=17)
                                                                                     ccgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCg                                                         >  2:1820455/3‑151 (MQ=17)
                                                                                       gCCACTCAGAGCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGcc                                                       <  2:653442/151‑1 (MQ=17)
                                                                                       gCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGcc                                                       >  1:4443579/1‑151 (MQ=17)
                                                                                         cactcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCa                                                     <  1:1360913/147‑1 (MQ=18)
                                                                                         cACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCa                                                     >  2:1299645/1‑151 (MQ=17)
                                                                                           ctcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGc                                                   >  1:7190507/3‑151 (MQ=17)
                                                                                            tcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGTACTCGGTGGGCGAACCCCGAGGATCGCGCAGcc                                                  >  1:3304194/2‑150 (MQ=17)
                                                                                            gcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCt                                                  <  2:5199194/150‑1 (MQ=18)
                                                                                             cAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAACCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCtt                                                 <  2:72308/151‑1 (MQ=18)
                                                                                               gAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTAt                                               >  2:5591999/1‑151 (MQ=18)
                                                                                                 aCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGa                                             >  1:4056280/1‑151 (MQ=18)
                                                                                                  ccAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAg                                            <  1:4667323/151‑1 (MQ=255)
                                                                                                   cAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgc                                           >  1:5073036/1‑151 (MQ=18)
                                                                                                    aGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgcg                                          <  1:4759933/151‑1 (MQ=18)
                                                                                                      cgcAAGCAGCGCGAGGCAGGCGACGAACTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCGCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                                        <  1:2912464/151‑1 (MQ=255)
                                                                                                         aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                                     >  1:144763/1‑151 (MQ=255)
                                                                                                         aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                                     <  2:2080555/151‑1 (MQ=255)
                                                                                                         aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                                     >  2:185084/1‑151 (MQ=255)
                                                                                                         aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                                     >  1:259376/1‑151 (MQ=255)
                                                                                                           gcagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAc                                   <  1:1397041/151‑1 (MQ=255)
                                                                                                            cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAcc                                  <  2:6636869/151‑1 (MQ=18)
                                                                                                            cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCGCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAcc                                  <  2:315575/151‑1 (MQ=18)
                                                                                                             agcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCa                                 >  1:800493/1‑151 (MQ=255)
                                                                                                               cgcgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAAt                               <  1:10098/151‑1 (MQ=255)
                                                                                                                gcgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAAtt                              <  2:680383/151‑1 (MQ=21)
                                                                                                                  gAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTgg                            >  2:5981005/1‑151 (MQ=255)
                                                                                                                  gAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTgg                            >  1:1624724/1‑151 (MQ=255)
                                                                                                                  gAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTgg                            >  1:5434409/1‑151 (MQ=255)
                                                                                                                   aggcaggcGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGc                           >  1:5099731/1‑151 (MQ=21)
                                                                                                                    ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                          <  1:6973138/144‑1 (MQ=21)
                                                                                                                    ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGCTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                          >  1:2233752/8‑151 (MQ=255)
                                                                                                                         gcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGCGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGa                     >  1:7169408/3‑151 (MQ=255)
                                                                                                                             cGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTc                 >  1:3759308/1‑151 (MQ=255)
                                                                                                                             cGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTc                 >  1:4965360/1‑151 (MQ=255)
                                                                                                                             cGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTc                 <  2:5920433/151‑1 (MQ=255)
                                                                                                                               aaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCgg               <  1:426523/149‑1 (MQ=255)
                                                                                                                               aaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCgg               <  1:6031759/149‑1 (MQ=255)
                                                                                                                                aaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGc              <  2:118376/150‑1 (MQ=25)
                                                                                                                                aaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGc              >  1:1960718/2‑151 (MQ=255)
                                                                                                                                 agcAAAACCATGAGCGCCGGCTCTTGGGCGAAGCAGGGTGGTGGTTTCTATTCCTGCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGcc             <  2:7068178/149‑1 (MQ=255)
                                                                                                                                 aTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGcc             >  1:3223122/1‑151 (MQ=255)
                                                                                                                                  tCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCa            >  1:1133660/1‑151 (MQ=255)
                                                                                                                                  tCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCa            >  1:3780260/1‑151 (MQ=255)
                                                                                                                                   cAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAg           >  1:3287692/1‑151 (MQ=21)
                                                                                                                                   cAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAg           <  2:3186072/151‑1 (MQ=255)
                                                                                                                                   cAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAg           >  1:6848161/1‑151 (MQ=25)
                                                                                                                                    aaaaCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAgg          <  1:5732789/148‑1 (MQ=21)
                                                                                                                                    aaaaCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGACCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCacg          >  2:3889098/4‑149 (MQ=255)
                                                                                                                                     aaaCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGc         <  1:6731112/149‑1 (MQ=255)
                                                                                                                                       aacATGAGAGCAGGCTTGTGGGCTAAGCAAGGTTGTAGTTGCTATACCTCGGACAATGGAAAACTGCTAGTGTGGAATCGGTGGGAGAACCCCGGGGATAGCCCAGCTTATGAGGGCCGTGACAGATTGGCTCAAGATTTCGGCCAGGccc       <  2:2345754/149‑1 (MQ=255)
                                                                                                                                       aCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGccc       <  1:18649/151‑1 (MQ=255)
                                                                                                                                        ccATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAGCCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCcg      <  2:4985520/151‑1 (MQ=25)
                                                                                                                                        ccATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCcg      >  2:6576731/1‑151 (MQ=255)
                                                                                                                                          aTGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCcgcg    >  1:6348901/1‑151 (MQ=255)
                                                                                                                                            gAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGCGCa  <  2:3695646/151‑1 (MQ=37)
                                                                                                                                               |                                                                                                                                                   
TCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGCGCA  >  CP009974/2738379‑2738669

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