Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,897,594 G→C intergenic (‑5372/‑3225) RPPX_08470 ← / → RPPX_08480 tryptophan synthase subunit beta/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,897,5940GC81.1% 6656.3 / inf 4472intergenic (‑5372/‑3225)RPPX_08470/RPPX_08480tryptophan synthase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (375/353);  new base C (1563/1553);  total (2255/2217)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 4480 total aligned reads displayed.

GCAACCGCCCCGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTC  >  CP009974/1897445‑1897741
                                                                                                                                                     |                                                                                                                                                   
gtcactgccccggacaacgtctacgtcggtagcaacgacccggttgtgatgtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGc                                                                                                                                                    >  1:5022302‑M2/129‑151 (MQ=255)
      gCCCCGGACAACGTCTATGTCGGCACTAACGCGCCGGTGGTCAATGCCATCGATGCCGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCa                                                                                                                                              >  2:1560833/1‑151 (MQ=255)
       caccagACAACTTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCAGCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCaa                                                                                                                                             <  2:2329355/146‑1 (MQ=255)
       caccagACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCaa                                                                                                                                             <  2:3153978/146‑1 (MQ=37)
         ccagACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACg                                                                                                                                           >  2:60526/4‑151 (MQ=255)
          cGGACAACGTCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGa                                                                                                                                          >  2:1712959/1‑151 (MQ=255)
           agACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCGGACAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGCGCCAGGCACCCCGGGCAACGaa                                                                                                                                         >  2:3789534/2‑151 (MQ=25)
                  gtctacgtcggtagcaacgacccggttgtgatgtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAggcgctc                                                                                                                                  >  2:3515977‑M2/111‑148 (MQ=255)
                  gTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGata                                                                                                                                  <  1:143237/151‑1 (MQ=255)
                   tCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGacat                                                                                                                                 >  1:6408128/1‑148 (MQ=255)
                   tCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGacac                                                                                                                                 >  2:5048607/1‑148 (MQ=255)
                   tCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGatat                                                                                                                                 >  1:1729183/1‑151 (MQ=25)
                    cTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACatc                                                                                                                                <  2:658394/151‑1 (MQ=255)
                            gtagcaacgacccggttgtgatgtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGt                                                                                                                        >  2:3389510‑M2/101‑151 (MQ=255)
                            gCACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAAACCCCCGTCAGCGCCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATACCATCAAGGt                                                                                                                        >  1:4261480/1‑151 (MQ=255)
                             cACTAACGCGCCGGTGGTCAATGCCATCGATGCCGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTc                                                                                                                       <  2:6150339/151‑1 (MQ=37)
                             cACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTc                                                                                                                       <  1:4609764/151‑1 (MQ=255)
                               gtAACGAGCCGGTGGTCAATGCCAGCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCAc                                                                                                                     <  2:3922525/150‑1 (MQ=255)
                                caatgcgCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGCCCAGGTCAcc                                                                                                                    >  2:2475545/7‑151 (MQ=37)
                                caacgcgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCAcc                                                                                                                    >  1:6733926/7‑151 (MQ=255)
                                 aaCCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCActa                                                                                                                   <  1:3977024/151‑3 (MQ=255)
                                  acgcgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCAt                                                                                                                  <  1:1279353/147‑1 (MQ=255)
                                  acgcgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCACGCACCCCGGGCAACGAAGGCGCTATCATCAACGTCACCac                                                                                                                  >  1:5800129/5‑150 (MQ=255)
                                    cAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCActatta                                                                                                                >  1:6841210/1‑146 (MQ=255)
                                      cccggttgtgatgtcgatcgcaaccgttgagggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACg                                                                                                              <  1:6260802‑M2/61‑1 (MQ=255)
                                        cggttgtgatgtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGc                                                                                                            >  1:1882640‑M2/89‑151 (MQ=255)
                                          gtggtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGctg                                                                                                          >  2:6472183/1‑149 (MQ=37)
                                           tggtCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGa                                                                                                         >  1:5913451/1‑151 (MQ=255)
                                            ggtCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGAc                                                                                                        <  2:3007610/151‑1 (MQ=37)
                                              tgatgtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCa                                                                                                      <  1:6033542‑M2/69‑1 (MQ=255)
                                              tCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCa                                                                                                      >  1:412217/1‑151 (MQ=38)
                                               ccaTGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAg                                                                                                     <  2:4185675/149‑1 (MQ=255)
                                                aaTGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGa                                                                                                    <  2:5919933/151‑1 (MQ=25)
                                                 aTGCCATCGATGCCGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGAc                                                                                                   <  2:77232/151‑1 (MQ=21)
                                                  gtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGact                                                                                                  <  1:2608469‑M2/73‑2 (MQ=255)
                                                     gatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCg                                                                                               <  2:1269097‑M2/76‑1 (MQ=255)
                                                      aTCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCgg                                                                                              <  1:2177997/151‑1 (MQ=37)
                                                      aTCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCgg                                                                                              <  1:5367457/151‑1 (MQ=21)
                                                      aTCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCgg                                                                                              >  2:2751753/1‑151 (MQ=21)
                                                       tcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGt                                                                                             <  1:1594306‑M2/78‑1 (MQ=255)
                                                        ggcaaccggtgaaggcgcgggtggcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTg                                                                                            <  2:4041893‑M2/79‑1 (MQ=255)
                                                         gACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGtag                                                                                           >  2:138947/1‑149 (MQ=21)
                                                          caaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGc                                                                                          >  2:1169971‑M2/71‑151 (MQ=255)
                                                           cgcgGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGcc                                                                                         >  1:2442129/1‑150 (MQ=37)
                                                              ggTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCgag                                                                                      >  1:5241755/1‑151 (MQ=37)
                                                              ggTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGtagccgag                                                                                      >  1:2926954/1‑144 (MQ=18)
                                                              cgTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGCCACCATCACGGCCGACCAGACGTCGGTGGCCgag                                                                                      >  2:5302405/2‑151 (MQ=255)
                                                               gTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCgaga                                                                                     >  2:1923757/1‑151 (MQ=18)
                                                               gTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCgaga                                                                                     <  2:5792397/151‑1 (MQ=18)
                                                                tCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGaa                                                                                    <  1:2338754/151‑1 (MQ=18)
                                                                  gaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACg                                                                                  <  2:5649499‑M2/89‑1 (MQ=255)
                                                                  aGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACg                                                                                  >  1:6614054/1‑151 (MQ=18)
                                                                   gcggcgCAGATGCATGGAAGTTCGAAAACCTGAACCCGCCCACGACTCCAGTCAGCACCGAGGTGACCGCCGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGt                                                                                 <  1:475426/151‑1 (MQ=18)
                                                                   aaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGAAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGt                                                                                 <  1:7084196‑M2/90‑1 (MQ=255)
                                                                    cggcgCTGATGCGTGGAAGGTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGACACCATCCCGGCCGACCAGACGTCGGTGCCCGAGAACgcc                                                                                >  2:5191723/1‑149 (MQ=255)
                                                                     ggcgCGGATGCATGGAAGTTCGAAAACCGGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCa                                                                               <  2:3283982/151‑1 (MQ=21)
                                                                     ggcgCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCCCGGCCGACCAGACCTCGGTGGCCGAGAACGTCa                                                                               >  2:888163/1‑151 (MQ=25)
                                                                     ggcgCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCa                                                                               >  1:3051442/1‑151 (MQ=18)
                                                                       cgcAGATGCGTGGAAGATCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGATGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCaaa                                                                             <  2:2069565/151‑2 (MQ=21)
                                                                        gcAGAGGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCaaac                                                                            <  2:4776185/151‑3 (MQ=21)
                                                                          tgATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAAccg                                                                          <  1:342883/150‑2 (MQ=17)
                                                                           gATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGccga                                                                         <  2:6881371/151‑3 (MQ=255)
                                                                               cgtGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCtt                                                                     >  2:3645544/3‑151 (MQ=18)
                                                                               cATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCtt                                                                     <  1:6118922/151‑1 (MQ=9)
                                                                                  ggAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCCTCACGGCCGACCAGACCTCGGCGGCCGCGAACGTCCAACCGACCTTCAc                                                                  >  2:4623234/1‑151 (MQ=37)
                                                                                  ggAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACATTCAc                                                                  >  2:5424076/1‑151 (MQ=18)
                                                                                  ggAAGTTCGAAAACCTGAACCTGGACAAAACCCCCGTCAGCACCGAGGTAACCAACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCAc                                                                  <  1:2031947/151‑1 (MQ=18)
                                                                                   gAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCAcc                                                                 <  2:6396397/151‑1 (MQ=255)
                                                                                     aGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGt                                                               >  1:5125031/1‑151 (MQ=18)
                                                                                        tCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCa                                                            >  1:5723724/1‑151 (MQ=2)
                                                                                              acctgaacctgGACAAAACCCCCGTCAGCACCGAGGTAACCGACTAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCaa                                                      >  1:5468089/1‑151 (MQ=18)
                                                                                                ctgaacctgGACAAGCCCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAAcc                                                    <  2:7044119/151‑1 (MQ=14)
                                                                                                ctgaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAAcc                                                    <  1:2764624/151‑1 (MQ=2)
                                                                                                 tgaacctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCa                                                   >  1:7203111/1‑151 (MQ=14)
                                                                                                 tgaacctgGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCACCCa                                                   >  2:6291112/1‑151 (MQ=17)
                                                                                                    acctgGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGcc                                                <  1:5577170/151‑1 (MQ=17)
                                                                                                        ggACAAGCCCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTg                                            <  2:6737954/151‑1 (MQ=11)
                                                                                                          aCAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTAAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTAGc                                          <  2:2885154/151‑1 (MQ=14)
                                                                                                            aaGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGctc                                        <  1:1002582/151‑3 (MQ=11)
                                                                                                              gaccccggTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAc                                      >  1:6497067/8‑151 (MQ=17)
                                                                                                               acccccgTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACg                                     >  1:783188/7‑151 (MQ=14)
                                                                                                                ccccggTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGa                                    >  2:5548460/6‑151 (MQ=11)
                                                                                                                   cggTCAGCACCGGGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCt                                 <  2:1019024/149‑1 (MQ=17)
                                                                                                                    ggTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTg                                >  1:3583396/2‑151 (MQ=17)
                                                                                                                      tCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgct                              >  2:4914711/1‑149 (MQ=17)
                                                                                                                         gCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGt                           <  2:4292092/151‑1 (MQ=11)
                                                                                                                         gCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCAACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGt                           >  2:2590823/1‑151 (MQ=17)
                                                                                                                           aCCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGCGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                         >  1:4895470/1‑151 (MQ=11)
                                                                                                                           aCCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGa                         >  1:3297974/1‑151 (MQ=17)
                                                                                                                            ctaccgTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAAGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAc                        >  2:7179904/6‑151 (MQ=11)
                                                                                                                                 gTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGa                   <  1:318892/151‑1 (MQ=14)
                                                                                                                                   gACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGc                 >  1:931771/1‑151 (MQ=11)
                                                                                                                                   cactgacgaGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGc                 <  1:1139657/147‑1 (MQ=17)
                                                                                                                                       gacgaGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGCGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTAAGcaaca             >  2:4445949/1‑151 (MQ=11)
                                                                                                                                          gaGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACg          >  1:3843782/1‑151 (MQ=17)
                                                                                                                                          gaGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaacaatg          >  1:4496262/1‑149 (MQ=14)
                                                                                                                                            gCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGct        <  1:6940720/151‑2 (MQ=14)
                                                                                                                                             ccAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGctc       <  1:3729985/151‑3 (MQ=14)
                                                                                                                                              cAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCa      >  1:4369249/1‑151 (MQ=17)
                                                                                                                                               aGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg     <  2:2201369/151‑1 (MQ=17)
                                                                                                                                                     |                                                                                                                                                   
GCAACCGCCCCGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTC  >  CP009974/1897445‑1897741

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: