Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,001,812 G→T 6.4% intergenic (+1234/‑452) PP_4407 → / → PP_4409 gamma‑aminobutyrate transporter/site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,001,8120GT 6.4% 285.9 / 5.0 94intergenic (+1234/‑452)PP_4407/PP_4409gamma‑aminobutyrate transporter/site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (35/52);  new base T (3/3);  total (39/55)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GACAACGCAACATCATCGAGCGCATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATATACGAACAATTGA  >  NC_002947/5001663‑5001960
                                                                                                                                                     |                                                                                                                                                    
gACAACGCAACATCATCGAGCGCATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAg                                                                                                                                                      >  1:899563/1‑150 (MQ=255)
      gCAACATCATCGAGCGCATGTTTGGCTGGCTGCAGACGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTagc                                                                                                                                                <  2:1378623/150‑1 (MQ=255)
      gCAACATCATCGAGCGCATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTagc                                                                                                                                                >  1:2185981/1‑150 (MQ=255)
          catcatCGAGCGCATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTagccagc                                                                                                                                            <  2:1676435/150‑1 (MQ=255)
              atcGAGCGCATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGa                                                                                                                                        <  2:1637113/150‑1 (MQ=255)
              atcGAGCGCATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGa                                                                                                                                        <  2:1637054/150‑1 (MQ=255)
                  aGCGCATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGGGTTTGCGACATCTCCTTTCGTACGGAGCGTAGCCAGCCCGactgc                                                                                                                                    >  1:1690731/1‑148 (MQ=255)
                     gcATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACTTCTCCTTTCGTATTGAGCCTAGCCAGCCCGACTACTgt                                                                                                                                 <  2:874378/150‑2 (MQ=255)
                     gcATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTACCCAGCCCGACTACTgg                                                                                                                                 >  2:660375/1‑150 (MQ=255)
                      cATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGt                                                                                                                                <  1:845326/150‑1 (MQ=255)
                              gctggctgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGTCTTGTTCCATGCGGTGTTTGCGACATCTCCCTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAg                                                                                                                        <  1:672279/150‑1 (MQ=255)
                              gctggctgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAg                                                                                                                        <  2:1362248/150‑1 (MQ=255)
                              gctggctgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAg                                                                                                                        <  1:1801013/150‑1 (MQ=255)
                              gctggctgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAg                                                                                                                        <  1:2022293/150‑1 (MQ=255)
                               ctgcctgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGTCTAGCCAGCCCGACTACTGGTTTCACGAgg                                                                                                                       <  2:1404095/150‑1 (MQ=255)
                                tggctgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGgc                                                                                                                      >  1:2383867/1‑150 (MQ=255)
                                  gctgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGgccg                                                                                                                    >  2:1903854/1‑150 (MQ=255)
                                    tgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTa                                                                                                                                                  >  2:2160148/1‑118 (MQ=255)
                                    tgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTa                                                                                                                                                  <  1:2160148/118‑1 (MQ=255)
                                    tgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTa                                                                                                                                                  <  1:2160194/118‑1 (MQ=255)
                                    tgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTa                                                                                                                                                  >  2:2160194/1‑118 (MQ=255)
                                    tgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGgccgcc                                                                                                                  >  2:2000276/1‑150 (MQ=255)
                                    tgAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGgccgcc                                                                                                                  >  2:1668425/1‑150 (MQ=255)
                                    tcaatgagaACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTAATTTAGTACGGAGCCTAGCCAGCCCGACTACTGGGTTCACGAGgccgcc                                                                                                                  >  1:1264880/6‑150 (MQ=255)
                                      aaaGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCaa                                                                                                                <  2:1547500/150‑1 (MQ=255)
                                        agagAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTATCGTAAGAACCCGATCCAGCCGGTGTACTGGTTTCACGAGGCCGCCACGt                                                                                                              >  2:391673/1‑150 (MQ=255)
                                         gagaACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTg                                                                                                                                  >  1:1681450/1‑129 (MQ=255)
                                         gagaACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTg                                                                                                                                  <  2:1681450/129‑1 (MQ=255)
                                           gaACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTAcc                                                                                                           >  1:1434103/1‑150 (MQ=255)
                                             aCCGCCGGATCGTAGCGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTCTCGTACGCAGCCTAGCCAGCCCGCCTACTAGCTTCACTAGGCCGCCAAGTACCGc                                                                                                         <  2:628287/150‑1 (MQ=255)
                                               cgccgGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCNCTg                                                                                                       <  2:1500909/150‑1 (MQ=255)
                                                 ccgGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGt                                                                                                     <  2:1833686/150‑1 (MQ=255)
                                                  cgGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTc                                                                                                    <  1:657147/150‑1 (MQ=255)
                                                  cgGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTc                                                                                                    >  1:2049993/1‑150 (MQ=255)
                                                   ggATCGTAACGCGCTTCGACACATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCg                                                                                                   >  1:630509/1‑150 (MQ=255)
                                                   ggATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCg                                                                                                   <  2:899563/150‑1 (MQ=255)
                                                    gATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCgg                                                                                                  >  2:1439170/1‑150 (MQ=255)
                                                          aaCGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAatat                                                                                            <  2:2049993/150‑1 (MQ=255)
                                                              cgcTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATAtttga                                                                                        <  2:2446370/150‑3 (MQ=255)
                                                              cgcTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTaa                                                                                        <  2:2383867/150‑1 (MQ=255)
                                                              cgcTTCGACAAATGCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTaa                                                                                        <  1:1536862/150‑1 (MQ=255)
                                                                   cGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTa                                                                                   <  2:547446/150‑1 (MQ=255)
                                                                   cGACAAATTCGCGAAAAGCTACGCCGCCACGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGc                                                                                                         >  1:918143/1‑128 (MQ=255)
                                                                   cGACAAATTCGCGAAAAGCTACGCCGCCACGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGc                                                                                                         <  2:918143/128‑1 (MQ=255)
                                                                     aCAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACa                                                                                 >  1:2006921/1‑150 (MQ=255)
                                                                       aaaTTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTAc                                                                                  >  1:2161317/1‑147 (MQ=255)
                                                                       aaaTTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTAc                                                                                  <  2:2161317/147‑1 (MQ=255)
                                                                        aaTTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACatta                                                                              <  1:2000276/150‑1 (MQ=255)
                                                                           tCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACGttatta                                                                           <  1:1903854/150‑1 (MQ=255)
                                                                            cgcgAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCTGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTaa                                                                          <  1:139923/150‑1 (MQ=255)
                                                                             gcgAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAAt                                                                         <  2:1434103/150‑1 (MQ=255)
                                                                             gcgAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAAt                                                                         >  1:507614/1‑150 (MQ=255)
                                                                             gcgAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAAt                                                                         <  1:2196428/150‑1 (MQ=255)
                                                                               gAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACa                                                                                      <  1:999547/135‑1 (MQ=255)
                                                                               gAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACa                                                                                      >  2:999547/1‑135 (MQ=255)
                                                                                aaaaGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCa                                                                      <  1:420015/150‑1 (MQ=255)
                                                                                 aaaGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCAc                                                                     <  2:67981/150‑1 (MQ=255)
                                                                                              cATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAATTCGGCTGGTCGTATATATTTAACAGTACATTATTAATGGTGCATTAATAATGta                                                        <  2:1264880/150‑2 (MQ=255)
                                                                                               aTGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCAACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACAGTATTAATGCACCATTAATAAAGTga                                                       >  2:389011/1‑149 (MQ=255)
                                                                                                 ggTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTATGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTAAAAAGGACATTATTAATCCACCATTAATAAAGTgcgg                                                     >  2:1098783/1‑149 (MQ=255)
                                                                                                     tcACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGAAATCAAATTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGcacca                                                 >  1:147971/1‑150 (MQ=255)
                                                                                                       aCTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGc                                                                                                         <  2:682352/92‑1 (MQ=255)
                                                                                                       aCTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGc                                                                                                         >  1:682352/1‑92 (MQ=255)
                                                                                                       aCTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACcaca                                               >  1:2083365/1‑150 (MQ=255)
                                                                                                         tGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACaaa                                             >  1:1020858/1‑150 (MQ=255)
                                                                                                         tGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACaaa                                             <  1:518066/150‑1 (MQ=255)
                                                                                                           tcTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGc                                           <  1:1545312/149‑1 (MQ=255)
                                                                                                           gCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTGAACAGTACAGTATTAATGCACCATTAAGAAAGTGCGcaccaccacgc                                           >  2:1084869/1‑145 (MQ=255)
                                                                                                            cTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGcc                                          >  2:202874/1‑150 (MQ=255)
                                                                                                             ttGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGATCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGccc                                         >  1:871759/1‑150 (MQ=255)
                                                                                                             ttGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGATCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGccc                                         >  1:871748/1‑150 (MQ=255)
                                                                                                             ttGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGATCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGccc                                         >  1:871760/1‑150 (MQ=255)
                                                                                                             ttGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGccc                                         >  1:641156/1‑150 (MQ=255)
                                                                                                                     tGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACaaa                                             >  1:739272/1‑138 (MQ=255)
                                                                                                                     tGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACaaa                                             <  2:739272/138‑1 (MQ=255)
                                                                                                                     tGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGc                                           >  1:820505/1‑140 (MQ=255)
                                                                                                                     tGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGc                                           <  2:820505/140‑1 (MQ=255)
                                                                                                                        ggTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAg                              <  2:711467/150‑1 (MQ=255)
                                                                                                                         gtgtTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGc                             <  1:1439170/150‑1 (MQ=255)
                                                                                                                           gtTTGCGACATCTCCTTTCGTACGGATCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCACGATCCTAAGCCt                           <  2:871760/150‑1 (MQ=255)
                                                                                                                           gtTTGCGACATCTCCTTTCGTACGGATCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCt                           <  2:871748/150‑1 (MQ=255)
                                                                                                                           gtTTGCGACATCTCCTTTCGTACGGATCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCt                           <  2:871759/150‑1 (MQ=255)
                                                                                                                                 gACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGacaaa                     >  2:1210677/1‑150 (MQ=255)
                                                                                                                                     tctcCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGAcaaacaaa                 >  1:962085/1‑150 (MQ=255)
                                                                                                                                       tcCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGa                         <  1:927264/140‑1 (MQ=255)
                                                                                                                                       tcCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGa                         >  2:927264/1‑140 (MQ=255)
                                                                                                                                       tcCTTTCGTACGGAGCCTAGCCAAGAGGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAata               >  1:1520053/1‑150 (MQ=255)
                                                                                                                                         cTTTCGGACGGAGCCTAGCCAGCCCGACTACTGGTGTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGacaaacatgcgct             <  2:1690731/150‑7 (MQ=255)
                                                                                                                                          tttCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCATGTTGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATATAc            <  1:391673/150‑1 (MQ=255)
                                                                                                                                          tttCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATATAc            <  2:630509/150‑1 (MQ=255)
                                                                                                                                           ttCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGATCGTCTGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATANACg           <  2:1520053/150‑1 (MQ=255)
                                                                                                                                            tCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATATACGa          <  1:1084869/150‑1 (MQ=255)
                                                                                                                                             cGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATATACGaa         <  1:660375/150‑1 (MQ=255)
                                                                                                                                                    aGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATATACGAACAATTGa  <  2:2006921/150‑1 (MQ=255)
                                                                                                                                                    aGCCTAGCCAACCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATATACGAACAATTGa  <  2:1965108/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACAACGCAACATCATCGAGCGCATGTTTGGCTGGCTGAAAGAGAACCGCCGGATCGTAACGCGCTTCGACAAATTCGCGAAAAGCTACGCCGCCATGGTCTCACTGGCTTGTTCCATGCGGTGTTTGCGACATCTCCTTTCGTACGGAGCCTAGCCAGCCCGACTACTGGTTTCACGAGGCCGCCAAGTACCGCTGGTCGGAAATATTTAACAGTACATTATTAATGCACCATTAATAAAGTGCGCACCACAAAGCCCAAAATCCAAAGCCTGACAAACAAATATACGAACAATTGA  >  NC_002947/5001663‑5001960

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