Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,048,596 G→A G32G (GGC→GGT gfcC ← capsule biosynthesis GfcC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,048,5960GA98.0% 176.3 / NA 49G32G (GGC→GGTgfcCcapsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/27);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGA  >  NC_000913/1048356‑1048842
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
gATATTGAGATTGAGTAACGGCTGGCGCACGGATTTAATCTTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGGCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGAt                                                                                                                                                                                                                                              <  2:264931/251‑1 (MQ=255)
 atatTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCTAACGGAAAGTGTCTGTTGTTCACCAGGCAGATa                                                                                                                                                                                                                                             <  2:376061/251‑1 (MQ=255)
        gatggagTAACTGCTGGCGCACGGATTTAATCGTGGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTc                                                                                                                                                                                                                                      <  1:322881/247‑1 (MQ=255)
         attgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTACGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATCAATAGTcc                                                                                                                                                                                                                                     >  2:40664/1‑250 (MQ=255)
           tgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCAc                                                                                                                                                                                                                                    <  1:19271/250‑1 (MQ=255)
                aaCTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTcc                                                                                                                                                                                                                              >  1:386758/1‑251 (MQ=255)
                    gctgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgc                                                                                                                                                                                                                          <  2:230348/251‑1 (MQ=255)
                          gcACGGATTTAATCGTTCCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCACGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCCGCAAGGCACCAGCCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCCACGGAACGTGTATGTTGTTGACCATGCAGCTAAATAGTCACCGTTCCTTgcgaaa                                                                                                                                                                                                                      >  2:437803/1‑246 (MQ=255)
                                aTTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATgg                                                                                                                                                                                                              >  1:187266/1‑251 (MQ=255)
                                       cGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCa                                                                                                                                                                                                       <  2:467267/251‑1 (MQ=255)
                                        gTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAACTAGTCACCGTTCCTTGCGCAAAGGCATGGGGGGTCAt                                                                                                                                                                                                      >  2:192611/1‑251 (MQ=255)
                                           gccgccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTAc                                                                                                                                                                                                   <  1:110359/251‑1 (MQ=255)
                                                cacaTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTaaa                                                                                                                                                                                              <  1:97862/251‑1 (MQ=255)
                                                  caTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGt                                                                                                                                                                                            <  1:452140/251‑1 (MQ=255)
                                                        cgtcAGCTTCTGCTGCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTAAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCt                                                                                                                                                                                       <  2:88507/250‑1 (MQ=255)
                                                                         cAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTg                                                                                                                                                                     >  2:15821/1‑251 (MQ=255)
                                                                         cAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTg                                                                                                                                                                     <  1:73343/251‑1 (MQ=255)
                                                                                caagcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTa                                                                                                                                                              <  1:499130/251‑1 (MQ=255)
                                                                                 aagcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTAt                                                                                                                                                             <  2:426027/251‑1 (MQ=255)
                                                                                      gggCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATg                                                                                                                                                        <  1:317044/249‑1 (MQ=255)
                                                                                            ggACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCgg                                                                                                                                                  <  1:143659/250‑1 (MQ=255)
                                                                                             gACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGAAAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGGTCCTTTCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCggg                                                                                                                                                 >  2:465454/1‑251 (MQ=255)
                                                                                                                   gcTTTAGCTTTCGCAGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGAGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAAAAGTCACCGTTCCTTGCGCAAAGGCATGGGGAGTCATTACGTAAAGTACGCTGGCAATGAAATACGAATGGAATTTATTCATGGCGGGGGCTTAAGAAACGGCATTTATACg                                                                                                                           >  2:199477/1‑251 (MQ=255)
                                                                                                                         gCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAAt                                                                                                                      >  1:481301/1‑250 (MQ=255)
                                                                                                                                ccgcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTg                                                                                                              <  1:511636/251‑1 (MQ=255)
                                                                                                                                ccgcAGAGGCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAACGTGTCTGTTGTTCACCAGGCAGATAAATAGGCACCGTTCCGTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTACTGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAAcccc                                                                                                                >  2:46272/1‑249 (MQ=255)
                                                                                                                                    agagTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGcc                                                                                                          <  1:47771/251‑1 (MQ=255)
                                                                                                                                      agTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCt                                                                                                         >  1:186944/1‑250 (MQ=255)
                                                                                                                                               gCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGcc                                                                                               <  1:36082/251‑1 (MQ=255)
                                                                                                                                               gCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGcc                                                                                               <  1:305009/251‑1 (MQ=255)
                                                                                                                                                  aGGCACCAGGCCAACAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGGTGTTCACCAGGCAGATAAATAGTCACCGTTACTTGCACAAAGGCATGGGGTGTCATTACGTAAAGAACGCTGGCAATGAAATACGACTGTAATTTATTCTTGGTGCGGGTTTACGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAAAATTTgtccac                                                                                            >  2:524262/1‑246 (MQ=255)
                                                                                                                                                     cACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTCACATTTGCCGACttt                                                                                         >  2:30890/1‑251 (MQ=255)
                                                                                                                                                             caccaGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGAACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAAcc                                                                                 >  2:350156/1‑251 (MQ=255)
                                                                                                                                                              accaGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCt                                                                                <  2:302436/251‑1 (MQ=255)
                                                                                                                                                              accaGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCt                                                                                <  1:483323/251‑1 (MQ=255)
                                                                                                                                                                          cgcgCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCg                                                                    >  1:443112/1‑251 (MQ=255)
                                                                                                                                                                           gcgcAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGa                                                                   <  1:156933/251‑1 (MQ=255)
                                                                                                                                                                                   gCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGAccc                                                           <  1:199477/251‑1 (MQ=255)
                                                                                                                                                                                       ttGTGTCACCAACTGGACAACATTTTCAACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGa                                                        >  1:404579/1‑250 (MQ=255)
                                                                                                                                                                                             cACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATgtgt                                                  <  1:70059/250‑1 (MQ=255)
                                                                                                                                                                                                     ggACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAAGAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCa                                           <  1:455456/249‑1 (MQ=255)
                                                                                                                                                                                                      gacaaaatgtagcACGGGGCCAACGGAAAGTGGCTGTTGTTCGCCAGGGAGATAAATAGTCACCGTTCCTTGCGCAAAGTCCTGGGGTGTCATTCCGTAAAGTACGCTGGCAATGAAATACGAGTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTTATTTCTACGGGAATAACCCCGGCGACTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCagtg                                        <  2:300858/239‑4 (MQ=255)
                                                                                                                                                                                                      gACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTg                                        >  1:268379/1‑251 (MQ=255)
                                                                                                                                                                                                                     ggggCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagt                          >  1:278325‑M1/1‑236 (MQ=255)
                                                                                                                                                                                                                      gggCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGGAAAGGCATGGGGTGTCATTACGTAAAGGACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaa                        <  2:103229‑M1/251‑17 (MQ=255)
                                                                                                                                                                                                                        gCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCTCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGATTTTGACGAACCGGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaacc                      <  2:187266‑M1/251‑19 (MQ=255)
                                                                                                                                                                                                                           aaCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaacc                      >  2:437881‑M1/1‑230 (MQ=255)
                                                                                                                                                                                                                                aaaGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagtt              >  1:190132‑M1/1‑225 (MQ=255)
                                                                                                                                                                                                                                            gttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGCACGCTCGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtgcacgat  >  2:451101‑M1/1‑213 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
GATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTGCCAGA  >  NC_000913/1048356‑1048842

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