Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,532,063 C→T intergenic (‑412/+374) sgcX ← / ← yjhP putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,0630CT98.2% 179.4 / ‑6.6 57intergenic (‑412/+374)sgcX/yjhPputative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (29/27);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GACCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCACTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTG  >  NC_000913/4531825‑4532306
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
gaCCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTATCTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATCCCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCATGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTAt                                                                                                                                                                                                                                         <  1:298372/251‑1 (MQ=255)
 aCCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATg                                                                                                                                                                                                                                        <  1:392173/251‑1 (MQ=255)
   cAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTACGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGATTAGATAATCAGCTCAAATCATTCTCCGCTTTATg                                                                                                                                                                                                                                        >  2:153321/1‑249 (MQ=255)
      tGATTATTGTTGAAATAATGAGGATGACTATAACGCGGTAACTGTTATAGCACCTCTTTGTTTCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCCTCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAACCATCGGATCTCCTGACTAGATTCAGAGGGAGTATATATAATCAGCTCAAATCATTCTCCGCTTTATGAGcc                                                                                                                                                                                                                                    <  2:369339/250‑1 (MQ=255)
        aTTATGGTTGAAATAATGAGGAGGAATAAAACGCGGTAACTGTTATAGCACCACTTTTTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAg                                                                                                                                                                                                                                 <  2:384438/251‑1 (MQ=255)
                   aataatGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGc                                                                                                                                                                                                                      <  1:259288/251‑1 (MQ=255)
                               gAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCACTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACg                                                                                                                                                                                                          >  1:212778/1‑251 (MQ=255)
                                   aaaaCGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGcc                                                                                                                                                                                                                                    >  2:17741/1‑221 (MQ=255)
                                   aaaaCGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGcc                                                                                                                                                                                                                                    <  1:17741/221‑1 (MQ=255)
                                    aaaCGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTg                                                                                                                                                                                                     <  1:174528/251‑1 (MQ=255)
                                        gcgGTAACTGTTATAGCACCTCTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTTATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGTAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACg                                                                                                                                                                                                 <  1:467131/251‑1 (MQ=255)
                                           gTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCCCGCtt                                                                                                                                                                                              >  2:90296/1‑251 (MQ=255)
                                            tAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCtt                                                                                                                                                                                              >  2:160051/1‑250 (MQ=255)
                                                 gTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGc                                                                                                                                                                                        >  1:219371/1‑251 (MQ=255)
                                                        caccacTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTtcctcc                                                                                                                                                                                 >  1:143857/1‑251 (MQ=255)
                                                         accacTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCg                                                                                                                                                                                >  2:73514/1‑251 (MQ=255)
                                                          ccacTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCgg                                                                                                                                                                               >  1:23727/1‑251 (MQ=255)
                                                                 gttgCAGGACGCTCTCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACt                                                                                                                                                                        <  1:221277/251‑1 (MQ=255)
                                                                    gCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGGCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCGCAAATCATGCTCCGCTTTATGAGCCCAGGCTGGGAGTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACtttt                                                                                                                                                                     <  1:43520/251‑1 (MQ=255)
                                                                             tCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGGATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGa                                                                                                                                                             >  1:39000/1‑250 (MQ=255)
                                                                                    tCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAgg                                                                                                                                                      >  1:108261/1‑250 (MQ=255)
                                                                                     cATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGa                                                                                                                                                     >  1:371673/1‑250 (MQ=255)
                                                                                         aCGTCCTGCACGGGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCTTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGTATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATTAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGaacct                                                                                                                                                 <  2:322779/250‑4 (MQ=255)
                                                                                               tGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAAAGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGGTTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCCCGCTTATGTCATTCCTCCGGGGCGACTTTTTGCCAGGTCTCAAGGAAGCTGGcaaac                                                                                                                                          >  2:284087/1‑247 (MQ=255)
                                                                                                   cGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTg                                                                                                                                      <  2:404596/251‑1 (MQ=255)
                                                                                                               aGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGcttctt                                                                                                                          <  1:73514/251‑1 (MQ=255)
                                                                                                                gccgctcCATCAGCGGATCTTCGAATGGTATGGCGCGTGGCAGGGCGCTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATGCAGAGGGATTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCTGATATCGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGTAGCTGGCTAACCCTGATACCGCttcttc                                                                                                                         <  2:57115/247‑1 (MQ=255)
                                                                                                                        aTCAGCGGATCGTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGTTCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGCCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGtt                                                                                                                 >  1:36085/1‑251 (MQ=255)
                                                                                                                          cAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGttt                                                                                                                <  1:62361/250‑1 (MQ=255)
                                                                                                                          cAGCGGATCTTAGAATGGAAGGGCGCGTGGCAGGTCGGTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTg                                                                                                               <  2:187900/251‑1 (MQ=255)
                                                                                                                             cGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTCTGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTg                                                                                                               >  2:122630/1‑248 (MQ=255)
                                                                                                                                               ggCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCt                                                                                           <  1:90296/250‑1 (MQ=255)
                                                                                                                                                     tGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTAc                                                                                    >  1:36740/1‑251 (MQ=255)
                                                                                                                                                         aGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACtttt                                                                                                                                                                     >  1:132008/1‑166 (MQ=255)
                                                                                                                                                         aGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACtttt                                                                                                                                                                     <  2:132008/166‑1 (MQ=255)
                                                                                                                                                                      ccAGCCTAAGGTAAGAAAACATCGTAGCTCCTGATTAGATTCAGAGGGAGTATAGAGAATCATCTCAAATCATTCTCCGGTTTATGATCCCAGGTTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCGTCCTCCGGGACGACGTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGAGATCGCATCTTCTCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACt                                                                   <  2:169252/251‑1 (MQ=255)
                                                                                                                                                                       cAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTg                                                                  >  1:201842/1‑251 (MQ=255)
                                                                                                                                                                       cAGCATAATGTAAGAAAACAACGGATCTCCTGAATAGATTCAGAGGGAGTAAAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTTGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGTCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCGGATACCGCTTCTTCCCATGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTg                                                                  <  2:208496/251‑1 (MQ=255)
                                                                                                                                                                              aTGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGcc                                                                                                                                                                                                                     <  2:358197/97‑1 (MQ=255)
                                                                                                                                                                              aTGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGcc                                                                                                                                                                                                                     >  1:358197/1‑97 (MQ=255)
                                                                                                                                                                              aTGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACACTATTGATCACTGGCACtgt                                                           >  2:321239/1‑251 (MQ=255)
                                                                                                                                                                              aTGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACtgt                                                           >  2:154421/1‑251 (MQ=255)
                                                                                                                                                                                  aaGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGGCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCACAGGGTTTGAGATTGCCGTTGCTGCAGCTGCATTACGCACACTATTGATCACTGGAACTGTGTGc                                                       >  2:150605/1‑251 (MQ=255)
                                                                                                                                                                                   agcaaaCATCGGATCTCCTGAATCGCTTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGCCTTTTTGCCAGGAATCACGGCAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGCGATTGCCGTTGCTGTCGCTGCATTACGCCCCATATTGATCACTGGCACTGTGTGc                                                       >  2:465704/4‑250 (MQ=255)
                                                                                                                                                                                     aaaaCATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGCTACCGCTTCTTCCCATGGTTTGAGATGGCCGTAGCTGACACTGCATTACGCCCAATATTGATCACTGGCACTGTGTGCATc                                                    >  2:271963/1‑251 (MQ=255)
                                                                                                                                                                                       aaCATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTg                                                  >  1:90566/1‑251 (MQ=255)
                                                                                                                                                                                         aaacggatctcccttatagtgtCAGAGGGAGTATAGATAACCAGATTAAATCATTCTCCGCTTTATGCGCCAATGCTGGCATTGCCGTTACATAACGGTTTGCACGCTTAAGTGCTTCCTCCGGGCTCTCTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATTCCGCTTCTTCCCAGGGTTTTAGATTGCCGTTGCTGACGGTGCATGACGTACAATATTGAGCTCCGGCACTGTGTGCAGCTggc                                                <  2:78932/230‑2 (MQ=255)
                                                                                                                                                                                                 ctcCTGAATAGATTCAGCGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACAt                                        >  2:100674/1‑251 (MQ=255)
                                                                                                                                                                                                   ccTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAATCCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACa                                         <  1:15244/248‑1 (MQ=255)
                                                                                                                                                                                                   ccTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGAc                                          <  1:335828/247‑1 (MQ=255)
                                                                                                                                                                                                                agagGGAGTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGc                         >  1:149717/1‑251 (MQ=255)
                                                                                                                                                                                                                       gTATAGATAATCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAg                  <  1:271555/251‑1 (MQ=255)
                                                                                                                                                                                                                               aaTCAGCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTc          <  1:321239/251‑1 (MQ=255)
                                                                                                                                                                                                                               aaTCAGCGCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGtt           <  1:66792/250‑1 (MQ=255)
                                                                                                                                                                                                                                    gCTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCAtgtt     >  2:176267/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cTCAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCAtgttg    >  1:232984/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cAAATCATTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCACAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTAGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTtgtg  >  1:440974/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
GACCAGTGATTATGGTTGAAATAATGAGGATGAATAAAACGCGGTAACTGTTATAGCACCACTTTGTTGCAGGACGCTCCCTCTTCATTACGTCCTGCACGTGGATACCCCAGCATCTCCATCAGCGGATCTTAGAATGTAATGGCGCGTGGCAGGTCGTTGGTTCCCAGCATAATGTAAGAAAACATCGGATCTCCTGAATAGATTCAGAGGGAGTATAGATAATCAGCTCAAATCACTCTCCGCTTTATGAGCCCAGGCTGGCATTGCCGTCAGATAACGTTTTGCACGCTTATGTGCTTCCTCCGGGCCGACTTTTTGCCAGGAATCAAGGAAGCTGGCTAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTG  >  NC_000913/4531825‑4532306

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: