Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 2,738,465 T→C A226A (GCT→GCC)  RPPX_12165 → benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099742,738,4650TC98.9% 123.9 / ‑5.5 89A226A (GCT→GCC) RPPX_12165benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (41/47);  total (41/48)
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.88e-01

ATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAT  >  CP009974/2738319‑2738612
                                                                                                                                                  |                                                                                                                                                   
agcaccGACATGATCGTCGACCAGTCGCCGGACGGTCTGGAAGTGCTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCa                                                                                                                                                 <  2:6350865/146‑1 (MQ=1)
                          gcgtGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCGACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCGGCGCAAGCAGCGCGAg                                                                                                                       <  2:3066841/148‑1 (MQ=0)
                                                  tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                                                               <  1:5635968/144‑1 (MQ=0)
                                                   ggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGa                                                                                              <  2:1230223/145‑1 (MQ=1)
                                                   ggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGAACGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGa                                                                                              <  2:5015331/145‑1 (MQ=0)
                                                    gagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAg                                                                                             >  1:799612/6‑151 (MQ=1)
                                                    gagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAg                                                                                             <  2:300526/146‑1 (MQ=1)
                                                     agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                                                                                            >  1:889692/5‑150 (MQ=1)
                                                      gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                           >  2:4785075/4‑149 (MQ=1)
                                                      gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                           <  2:3059201/148‑3 (MQ=1)
                                                      gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                           >  1:989278/4‑149 (MQ=1)
                                                      gcgaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                           <  2:4745553/148‑3 (MQ=1)
                                                                          gggtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAg                                                                       <  2:1749783/147‑1 (MQ=0)
                                                                           ggtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTCAGCAgg                                                                      >  2:5591227/4‑151 (MQ=0)
                                                                                       ctgaCCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCt                                                          <  2:227142/148‑1 (MQ=0)
                                                                                            ggcggagaaCGGCGCCGCCGGCTATCACGGTAGCTCCGTGCGCTGGAACTACGGCGCCAGTCGGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTAGTcc                                                     <  2:5441757/144‑1 (MQ=0)
                                                                                                gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                 >  2:444382/4‑151 (MQ=0)
                                                                                                gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                 <  2:2985751/148‑1 (MQ=1)
                                                                                                gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGCGGTTTCTATTCCtccc                                                 >  2:6034326/4‑148 (MQ=1)
                                                                                                gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGGGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGCGGTGGTTTCTATTCCTTCg                                                 >  2:1579574/4‑151 (MQ=1)
                                                                                                      ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                           >  2:3189141/1‑151 (MQ=2)
                                                                                                      ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                           <  1:167707/151‑1 (MQ=2)
                                                                                                      ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                           >  2:5233361/1‑151 (MQ=2)
                                                                                                       gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg                                          >  2:3847568/1‑151 (MQ=11)
                                                                                                        cgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGa                                         >  2:2269449/1‑151 (MQ=11)
                                                                                                         gcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGac                                        <  1:6186046/151‑1 (MQ=2)
                                                                                                          ccGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGaca                                       >  1:2309547/1‑151 (MQ=2)
                                                                                                           cGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTGTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGacac                                      >  2:6330214/1‑151 (MQ=11)
                                                                                                           cGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGacac                                      >  2:5934668/1‑151 (MQ=2)
                                                                                                           cGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGacac                                      >  2:3111905/1‑151 (MQ=2)
                                                                                                             aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                    <  2:5273855/151‑1 (MQ=2)
                                                                                                                 cTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTGCTATGCCTTCGACCATGGACACctgctg                                >  1:5431143/1‑151 (MQ=11)
                                                                                                                    tcatgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctc                             <  1:2180544/144‑1 (MQ=11)
                                                                                                                     catgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctct                            <  1:4840473/145‑1 (MQ=11)
                                                                                                                     catgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctct                            >  1:1359000/7‑151 (MQ=11)
                                                                                                                     catgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctct                            <  1:723442/145‑1 (MQ=11)
                                                                                                                      atgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTg                           <  2:1541052/146‑1 (MQ=11)
                                                                                                                      atgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTg                           >  2:5920953/6‑151 (MQ=11)
                                                                                                                       tgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTgg                          >  1:984601/5‑151 (MQ=11)
                                                                                                                       tgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTgg                          >  1:5035408/5‑151 (MQ=11)
                                                                                                                        gttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGa                         <  1:5155402/148‑1 (MQ=11)
                                                                                                                         ttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGAc                        >  1:5394370/3‑151 (MQ=11)
                                                                                                                          taGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACt                       <  1:1701059/150‑1 (MQ=14)
                                                                                                                           aGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTc                      >  2:5275618/1‑151 (MQ=11)
                                                                                                                           aGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTc                      <  2:3138730/151‑1 (MQ=14)
                                                                                                                           aGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTc                      >  1:1404316/1‑151 (MQ=14)
                                                                                                                            gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                     <  2:1831970/151‑1 (MQ=14)
                                                                                                                            gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                     >  1:1333077/1‑151 (MQ=14)
                                                                                                                            gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                     <  2:901153/151‑1 (MQ=14)
                                                                                                                            gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                     <  1:4982551/151‑1 (MQ=14)
                                                                                                                            gCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCg                     <  1:1611668/151‑1 (MQ=14)
                                                                                                                             gtacGTGCACTGGAACTACGACGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCgg                    <  2:1282283/148‑1 (MQ=14)
                                                                                                                             cTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCgg                    <  2:496252/151‑1 (MQ=14)
                                                                                                                              tCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGt                   >  1:5950414/1‑151 (MQ=14)
                                                                                                                               ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                  >  1:1366002/1‑151 (MQ=14)
                                                                                                                               ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                  <  2:2268092/151‑1 (MQ=14)
                                                                                                                               ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                  >  1:6446310/1‑151 (MQ=14)
                                                                                                                                 gggCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTggg                <  1:3233661/149‑1 (MQ=14)
                                                                                                                                 gTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTggg                >  2:1007726/1‑151 (MQ=14)
                                                                                                                                 gTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTggg                <  1:2269449/151‑1 (MQ=14)
                                                                                                                                 gTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTgcg                >  2:250729/1‑149 (MQ=14)
                                                                                                                                  tGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGc               <  1:3415349/151‑1 (MQ=14)
                                                                                                                                   gCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCg              <  2:112146/151‑1 (MQ=14)
                                                                                                                                   gCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCg              >  1:229187/1‑151 (MQ=14)
                                                                                                                                   gCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCg              <  2:5866107/151‑1 (MQ=14)
                                                                                                                                   gCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCg              <  2:1027647/151‑1 (MQ=14)
                                                                                                                                   gCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCCCTGGACTCGGTGGGCg              >  1:5507227/1‑151 (MQ=14)
                                                                                                                                    cACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGa             <  2:1313550/151‑1 (MQ=14)
                                                                                                                                    cACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGa             >  1:141045/1‑151 (MQ=14)
                                                                                                                                    cACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGATATCTATTCCTTCGACCACGGACACCTGCTGCTCTGGGCTCGGTGGGCgg             >  2:4132451/1‑150 (MQ=11)
                                                                                                                                     aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa            <  1:4535958/151‑1 (MQ=14)
                                                                                                                                     aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa            <  1:1022231/151‑1 (MQ=14)
                                                                                                                                       tGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGCCTCGGTGGGCGAAcc          >  2:921242/1‑151 (MQ=14)
                                                                                                                                       tGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcc          >  2:921220/1‑151 (MQ=14)
                                                                                                                                       gggAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcc          <  2:2762877/150‑1 (MQ=14)
                                                                                                                                        ggAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAccc         <  2:1590873/151‑1 (MQ=17)
                                                                                                                                         gagcgaCGCCGCCACGCAGAACCAGCGCAAGCAGCGCGAGCCGGGCTACGAAATCAAAACCATGAGCGCCGGCACTTGGGCTAGGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc        <  2:4672185/146‑1 (MQ=14)
                                                                                                                                         gAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc        <  1:4841094/151‑1 (MQ=17)
                                                                                                                                         gAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc        <  1:4840984/151‑1 (MQ=17)
                                                                                                                                         gAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc        >  1:5209068/1‑151 (MQ=17)
                                                                                                                                           aCTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGa      <  1:1537785/151‑1 (MQ=17)
                                                                                                                                            cTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAg     <  1:6473766/151‑1 (MQ=17)
                                                                                                                                            cTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAg     >  2:2610078/1‑151 (MQ=14)
                                                                                                                                             tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg    <  1:5231904/151‑1 (MQ=17)
                                                                                                                                             tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg    >  2:4604187/1‑151 (MQ=14)
                                                                                                                                             tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg    <  1:3184615/151‑1 (MQ=17)
                                                                                                                                             tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg    >  1:2235362/1‑151 (MQ=17)
                                                                                                                                              aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGa   <  2:3147537/151‑1 (MQ=17)
                                                                                                                                              aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTAGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAACCCGAGGa   >  2:2632952/1‑151 (MQ=17)
                                                                                                                                               cgccgcCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCGGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt  >  2:5786471/1‑151 (MQ=17)
                                                                                                                                               cgccgcCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt  >  1:4207250/1‑151 (MQ=17)
                                                                                                                                               cgccgcCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt  >  2:2618806/1‑151 (MQ=17)
                                                                                                                                                  |                                                                                                                                                   
ATCATCGACATGATCGTCGATCAGTCGCCTGAAGGCCTGGAAGTGCTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAT  >  CP009974/2738319‑2738612

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