Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 478,339 G→T 12.5% intergenic (+84/+26) PP_0393 → / ← cca 2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine pyrophosphokinase/multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947478,3390GT12.5% 368.2 / 3.6 139intergenic (+84/+26)PP_0393/cca2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine pyrophosphokinase/multifunctional tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (60/55);  new base T (4/19);  total (64/75)
Fisher's exact test for biased strand distribution p-value = 2.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

CAGATCGATCAGGTACTGGCCCCCGTGGCCTTTGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTCACATCCACCGCTCGCG  >  NC_002947/478190‑478477
                                                                                                                                                      |                                                                                                                                          
gtcaggcactgggcgaagcactcaaaggtgaacggctgaaagcgctcgaggcttaccagcgcggttagacctacCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCggg                                                                                                                                            >  1:3018838‑M1/75‑150 (MQ=255)
                tGGCCCCCGTGGCCTTTGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGGAAAACCCCCCCCACGACGGGCTGAGCCCCACCCCTGGTGTGTGGCCCCGGGAGGGGCGGGGTTAACCGCCCAGAgg                                                                                                                            >  1:4113082/1‑150 (MQ=255)
                    ccccGTGGCCTTTGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGGTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                        >  2:1508994/1‑150 (MQ=255)
                      ccGTGGCCTTTGAGTGGCGCGGTTTGCTGTTGACGCCCGCTGGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTa                                                                                                                      <  2:6731733/150‑1 (MQ=255)
                       cGTGGCCTTTGTGGGGCGGGGGTTGCTGTTGACGCCCGGTTGAGCCTGTATGGGCCTCGTGGGGGTTAAACCCGCTCCTACGACGGGCTGCGCCAAACCTCTTATGTTGCCCCCCTTTCGGCGGGGTTTTCCCCGCGAAGAGGCCGGTAg                                                                                                                     <  2:2612122/150‑1 (MQ=255)
                         tGGCCTTTTTGGGGCGGGGGTTGCATTTCAGGCCCGTTGGCGCCTTTATCGGCCTCTTCGGGGTTAAACCCGCACCTACGCGGGTCTGACCCCAACCTCTGTTGTGTGCCCCGTTTTGGTGGGGTTTTACCCGCGAAGAGGCCGGTAGGt                                                                                                                   <  2:5207383/150‑1 (MQ=255)
                         tGGCCTTTTTGGGGCGCGGTTTGCCGTTGACGCCCGCTTGCGCCTTTAGGGGCCTCTGCGGGGTTAAACCCGCTCCTACGCGGGGCTGAGCCAAACCTCTGTTGTGGGCCCCCTGTCGGCGGGGTTTTACCCGCGAAGAGGCCGGTAGGt                                                                                                                   <  2:2599859/150‑1 (MQ=255)
                          ggCCTTTGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAAACCGCACCTACGACGGGCTGAGCCCAACCTCTGATTTGTGCCCCCGGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGTAGGTc                                                                                                                  >  1:6760186/1‑150 (MQ=39)
                             cTTTGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTaa                                                                                                               >  2:5578458/1‑150 (MQ=255)
                              tttGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAAACCGCACCAAAGACGGGCTGAGCCCAACCACTGATATTTTCCCCCAGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGTAGGTCTAAc                                                                                                              >  2:1951205/1‑150 (MQ=255)
                              tttGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAAACCCCTACAACGACCGGCTGAGCCCAAACTCTGATGTGTGCCCCCAGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAAc                                                                                                              >  2:3021254/1‑150 (MQ=255)
                              tttGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAAACCCCGCCGACGACGGGCTGAGCCCAAACACAGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAAc                                                                                                              >  2:3021259/1‑150 (MQ=38)
                                     ggCGCGGGTTGCTGTGGCCGCCCGCTTGAGCCTGTATCGGCCGCTCCGCGGGTAACCCCGCTCCCACGAGGGGCTGCGCCCACCCTCTGAGGTGGGCCCCCTGTCGGAGCGGGTTTACCGGCGAAGAGGCCGGTAGGTCTAACCGCGCTg                                                                                                       <  2:6573710/150‑1 (MQ=255)
                                             tcgaggcttaccagcgcggttagacctacCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCggg                                                                                                                                            >  1:1921818‑M1/30‑105 (MQ=255)
                                             tcgaggcttaccagcgcggttagacctacCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCggg                                                                                                                                            <  2:1921818‑M1/76‑1 (MQ=255)
                                                cttttGCCGCCCTCTTGTGCGTTTATGGGCCTCTGGGGGGTTAACCCCGCTCCTCCGCGGGGCTGCGCCAACCCTGTTTTGTGGGCCCCCTGTCGGTGGGGGTTAACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGa                                                                                            <  2:2449324/147‑1 (MQ=255)
                                                ctgTTGCCCCCCGTTTGAGCCTTTATCGGCCTCTCCGGGGGTAAACCCGCCCCTCCGCCGGGCTGCGCCCAACCTCTGTTGTGGGCCCCGTTTTGGCGGGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGa                                                                                            <  2:781298/148‑1 (MQ=255)
                                                cAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGCGGGCAGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGCGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGa                                                                                            >  2:3502395/1‑150 (MQ=255)
                                                     gccGCCCGCTTGCGCCTGTATCGGCCTCTTCGGGGGTAACCCCGCTCCTACGCCGGGCTGCGCCACACCTCTGATGTGGGCCCCCGGTAGGAGCGGGTTTACCCGCGACGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCt                                                                                       <  2:6559838/148‑1 (MQ=255)
                                                     gACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCGACGACGGGCGGAGCCCAACCGCTGATGGGTGCCCCCGGTAGGAGCGGGTTTACCCCCGAAGAGGCCGGTAGGTCTAACCCCGCTGGTAAGCCGCGaccccg                                                                                       >  2:2203114/1‑145 (MQ=38)
                                                         cccGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCa                                                                                   >  2:5705798/1‑150 (MQ=255)
                                                            gCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGcc                                                                                <  2:4465640/150‑1 (MQ=255)
                                                              ttGCGCGTGTATCGGCCGCTCGGCGGGTAAACCCGCTCCTACGAGGGGCTGCGCCAAACCTCTGTGGTGTCCCCCCTTTTGGAGCGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGt                                                                              <  2:3156318/150‑1 (MQ=255)
                                                              ttGCGCCTTTCTCGGCCTCTTCGGGGTAAACCCCGCCCCTCCGACGGGCTGCGCCCACCGTCTGTTGTGGGCCCCCTGTCGGTGGGGTTTAACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGt                                                                              <  1:6245689/150‑1 (MQ=255)
                                                              ttGAGCCTGTATGGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGTTTTACCCGCGATGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGt                                                                              <  2:3155822/150‑1 (MQ=255)
                                                                  gCCTGTATCGGCCTCTTCGGGGGTAAACCCGCTCCTACGACGGGCTGAGCCAAACCTCTGATGTGTGCCCCCTGTTGGAGGGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCAc                                                                          <  2:5256591/150‑1 (MQ=39)
                                                                    cTGTATCGGCCTCTTCGCGGGTAAACCCCCTCCTACGACGGGCTGAGCCCAACATCTGATGTGTGCCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCt                                                                        >  2:6768829/1‑150 (MQ=255)
                                                                      gTATCGGCCTCTTCGCGGTTAAACCCGCTCCTACGACGGGCTGCGCCCACCCTGTGATGTGTGCCCCGTGTTGGAGCGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCttt                                                                      <  2:6383209/150‑1 (MQ=38)
                                                                      gTATCGGCCTCTGCCCGGTAAACCCCCCTCATACGCCGGGCTGCCCCAAACCTCTGATGTGTCCCCCCTGCAGGAGCGGTTTAACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCttt                                                                      <  1:6800757/150‑1 (MQ=255)
                                                                      gTATCGGCCTCGTCGCGGGTACCCGGGCTCCTCCGTCGGGCGGTGCCCACCCGCTGATGGGGGCCCGCTGTTGGTGCGGGTTTACCCGCGATGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCttt                                                                      <  2:4839574/150‑1 (MQ=255)
                                                                      gTATCGGCCTCGTCGCGGGAAAACCCGCTCCTACGAGGGGCTGAGCCCAACCTCTGTTGTGTGCCCCCTCTAGGAGCGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCttt                                                                      <  2:6383099/150‑1 (MQ=255)
                                                                        ttCGGCCTCTTGGGGGGAAACCCCGCTCCTACGAGGGGCTGAGCCAAACCTGTGTTGTGTGCCCCGTTTAGGCGGGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGa                                                                    <  1:5584613/149‑1 (MQ=255)
                                                                        acCGGCCTCTCCGCGGTAAACCCCGCTCATAGCAGGGGCTGACCCCAACCTCTGTGGTGGGCCCCCTGTAGGCGGGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGa                                                                    >  1:6771437‑M1/3‑150 (MQ=255)
                                                                         tCGGCCTCTTGGGGGGAAAACCCGCTCCTCGGACGGGCTGAGCCCACCCTCTGTTGTGGGCCCCCTGTAGGCGGGGGTTCACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAg                                                                   <  1:6285693/150‑1 (MQ=255)
                                                                          CGGCCTCTTCGGGGGTAAACCCGCTCCTACGCGGGGCTGAGCCCACCCTCTGAGGTGTGCCCCCTGTAGTAGCGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGt                                                                  >  2:4885744‑M1/1‑150 (MQ=255)
                                                                             CCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCAGGAGGAGCGGGGTTACCCGCGAAGAGGCCGGTAGGTCGAACCGCGCTGGTAAGCCTCGAGCGCTGTCAGCCGTTCACCCTTGCGTGCt                                                               <  2:3862058‑M1/150‑1 (MQ=255)
                                                                               tctTCGCGGGTAAACCCGCTACTAAGACGGGCTTAGCCCAACCTCTGTTGTGTTCCCCCTGTAGGATCGGGTTTAACCGCGAAGAGGCCGGTAGGTCTAAACGCCCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCGTTGTGTGCTTc                                                             >  2:1718658/1‑150 (MQ=255)
                                                                               TCTTCGCGGTAAAACCCGCTCCAACGACGGGCTCACCCCAACCTCTGATGTGTCCCCCCTGTAGGACCGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTc                                                             >  2:4353880‑M1/1‑150 (MQ=255)
                                                                               TCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCGCTGTTGTGTGCCCCCGGTAGGAGCGGGTTTACCCGCGAAGGGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCGTTGAGTGCTTc                                                             <  2:6076421‑M1/150‑1 (MQ=255)
                                                                                 TTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAAACTCTGATGTGTGCCCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTGAGCCTTTGAGTGCTTCGc                                                           <  2:5974125‑M1/150‑1 (MQ=255)
                                                                                  TCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGTTGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACATTTGAGTGCTTCGcc                                                          <  2:6130531‑M1/150‑1 (MQ=255)
                                                                                    GCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCa                                                        <  2:2528536‑M1/150‑1 (MQ=255)
                                                                                    GCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCa                                                        <  2:3498500‑M1/150‑1 (MQ=255)
                                                                                    GCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCa                                                        <  2:3499240‑M1/150‑1 (MQ=255)
                                                                                      gggAAACCCCGCCCCTCCGCGGGGCGGCGCCCACCCTGGGTGGGGGGCCCCCTGTGGGCGGGGGTTTCCCCGCGACGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGt                                                      <  1:3862058/150‑1 (MQ=255)
                                                                                      GGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGt                                                      <  2:3558561‑M1/150‑1 (MQ=255)
                                                                                      GGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGt                                                      >  2:5712958‑M1/1‑150 (MQ=255)
                                                                                       ggttcaccccGCTCCAACGCCGGGCTGCCCCCACCCTCTGTTGTGGGCCCCCTGCCGGCGGGGTTTTCCCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGt                                                      <  1:2783532/143‑1 (MQ=255)
                                                                                        GTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGc                                                    <  2:2167896‑M1/150‑1 (MQ=255)
                                                                                            ACCCGCTCCTACGAAGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCATCGAGCGCTTTAAGCCGTTCACCTTAGAGTGCTTCGCCCAGTGCCTGa                                                <  2:5854096‑M1/150‑1 (MQ=255)
                                                                                              CCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAc                                               <  2:449879‑M1/150‑1 (MQ=255)
                                                                                              CCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAc                                               <  2:3250542‑M1/150‑1 (MQ=255)
                                                                                               CCGCTCCTACGCGGGGCTGCGCCCACCCTCTTTTGTGTCCCCCGTGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAcc                                              >  2:1893476‑M1/1‑150 (MQ=255)
                                                                                               CCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAcc                                              <  1:1204764‑M1/150‑1 (MQ=255)
                                                                                               CCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAcc                                              <  2:4495198‑M1/150‑1 (MQ=255)
                                                                                               CCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAcc                                              <  2:2512605‑M1/150‑1 (MQ=255)
                                                                                               CCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAcc                                              <  1:300883‑M1/150‑1 (MQ=255)
                                                                                               CCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAcc                                              <  2:5843468‑M1/150‑1 (MQ=255)
                                                                                               CCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAcc                                              <  1:1407453‑M1/150‑1 (MQ=255)
                                                                                               CCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGAcc                                              <  1:1204905‑M1/150‑1 (MQ=255)
                                                                                                     CTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAg                                        <  1:2563662‑M1/150‑1 (MQ=255)
                                                                                                     CTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAg                                        <  2:4861854‑M1/150‑1 (MQ=255)
                                                                                                     CTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAg                                        <  1:3103602‑M1/150‑1 (MQ=255)
                                                                                                     CTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAg                                        <  1:3524221‑M1/150‑1 (MQ=255)
                                                                                                     CTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAg                                        <  1:3742480‑M1/150‑1 (MQ=255)
                                                                                                      TACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGGAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGc                                       <  1:1667217‑M1/150‑1 (MQ=255)
                                                                                                      TACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGc                                       <  2:388684‑M1/150‑1 (MQ=255)
                                                                                                      TACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGc                                       <  1:2462821‑M1/150‑1 (MQ=255)
                                                                                                       ACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGcc                                      >  1:5843468‑M1/1‑150 (MQ=255)
                                                                                                        CGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGccc                                     <  1:3604053‑M1/150‑1 (MQ=255)
                                                                                                         gacgGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCg                                    >  2:2067872/1‑150 (MQ=255)
                                                                                                         gacgGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCg                                    >  1:1668225/1‑150 (MQ=255)
                                                                                                          acgGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGCTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCgg                                   >  1:3862339/1‑150 (MQ=255)
                                                                                                             ggCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCtt                                >  1:4711869/1‑150 (MQ=255)
                                                                                                              gCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTg                               >  1:5640490/1‑150 (MQ=255)
                                                                                                                tGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCa                             >  2:2699487/1‑150 (MQ=255)
                                                                                                                 gAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCAc                            >  2:5054622/1‑150 (MQ=255)
                                                                                                                    cccAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCaa                         >  2:6763264/1‑150 (MQ=255)
                                                                                                                    cccAACCTCGGATGTGTGCCCCCTGTTGGAGGGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCaa                         <  2:287424/150‑1 (MQ=255)
                                                                                                                     ccAACCTCTTATGTGGGCCCCCTGTCGGTGCGGGTTTCCCCGCGACGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAAc                        <  1:5716029/150‑1 (MQ=255)
                                                                                                                     ccAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAAc                        >  2:1297286/1‑150 (MQ=255)
                                                                                                                      cAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACg                       >  2:3169900/1‑150 (MQ=255)
                                                                                                                      cAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACg                       <  2:2102437/150‑1 (MQ=255)
                                                                                                                      cAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACg                       >  2:238020/1‑150 (MQ=255)
                                                                                                                       aaCCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACgg                      >  1:3658075/1‑150 (MQ=255)
                                                                                                                        aCCTCTTATTGGTGCCCCATTTAGCAGCGGCTTTACCCGCTAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGc                     <  2:2944454/150‑1 (MQ=255)
                                                                                                                           tctGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTc                  >  1:4131613/1‑150 (MQ=255)
                                                                                                                           tctGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTc                  >  1:4838531/1‑150 (MQ=255)
                                                                                                                             tGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTcac                >  1:872050/1‑150 (MQ=255)
                                                                                                                             tGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTcac                >  1:5824449/1‑150 (MQ=255)
                                                                                                                             tGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTcac                >  2:6253260/1‑150 (MQ=255)
                                                                                                                              gATTTGTGCCCCCTGTAGCTGGGGTTTTCCCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTcaca               <  2:1140941/150‑1 (MQ=255)
                                                                                                                              gATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTcaca               >  1:2093282/1‑150 (MQ=255)
                                                                                                                                  tgtgCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTCACATCCa           >  2:2331027/1‑150 (MQ=255)
                                                                                                                                   gtgCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTCACATCCAc          >  1:5701040/1‑150 (MQ=255)
                                                                                                                                        cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTCACATCCACCGCTc     >  1:5288127/1‑150 (MQ=255)
                                                                                                                                        cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTCACATCCACCGCTc     >  1:6814150/1‑150 (MQ=255)
                                                                                                                                         ccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGt                                           <  2:3724081/111‑1 (MQ=255)
                                                                                                                                          cTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTCACATCCACCGCTcgc   <  2:516749/150‑1 (MQ=255)
                                                                                                                                           tGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTCACATCCACCGCTcgcg  >  2:5927002/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                          
CAGATCGATCAGGTACTGGCCCCCGTGGCCTTTGAGTGGCGCGGGTTGCAGTTGACGCCCGCTTGAGCCTGTATCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGGCTGAGCCCAACCTCTGATGTGTGCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAGGTCTAACCGCGCTGGTAAGCCTCGAGCGCTTTCAGCCGTTCACCTTTGAGTGCTTCGCCCAGTGCCTGACCGGTCAGCCCGGCTTGCACCAACGGCTTCACATCCACCGCTCGCG  >  NC_002947/478190‑478477

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