Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,819 G→C 21.1% A291G (GCG→GGG)  cynR ← transcriptional activator of cyn operon, autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8190GC21.1% 62.5 / 3.4 39A291G (GCG→GGG) cynRtranscriptional activator of cyn operon, autorepressor
Reads supporting (aligned to +/- strand):  ref base G (19/11);  new base C (4/4);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGT  >  NC_000913/357581‑358046
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      
gtgtCAGCGGCTACCGTGCTTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTGCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACGGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                        <  1:1395457/239‑1 (MQ=255)
 tcgcacggtcTACCGTAATTCAACCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                        <  2:342382/229‑1 (MQ=255)
 ggccAGCGGCTACCCTGCTGCATCCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                        <  2:1319199/235‑1 (MQ=255)
   tCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACcgcg                                                                                                                                                                                                                                      >  1:868183/1‑238 (MQ=255)
      gCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACTCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAACGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGGGATTCATTTCCGCCAACAACCGCGCAtt                                                                                                                                                                                                                                  >  1:760715/1‑239 (MQ=255)
      gCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTACAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCCGACCTGTGGCAGCGG‑CCACCGGGATTCATTTTCGCCCACAACcgcacctt                                                                                                                                                                                                                                  >  2:337930/1‑234 (MQ=21)
        ggCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                                                                                                                                                   <  1:871918/236‑1 (MQ=255)
        ggCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                                                                                                                                                   <  1:234157/236‑1 (MQ=255)
                       ttcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCCGACCTGTGGCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAAGTGc                                                                                                                                                                                                                 >  1:559815/1‑239 (MQ=255)
                           cGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATCAGCGTAGCGCATCAGGGAATTCCAGCCGCAGACCTGTGAAAGCGG‑CTACCGTGATTCAATTCCGCCAAACACCCCGCATTTATCCAACGCCATCTGCaaaa                                                                                                                                                                                                             >  1:1343726/1‑239 (MQ=21)
                            gCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGGGATTCATTTCCCCCAACAACCGCGCATTTTTCCCACGCCCTGTGCaaaaa                                                                                                                                                                                                            >  1:330663/1‑239 (MQ=255)
                                    ccGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCCTATCCCGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGCTTCCTTTCCGCCAACACCCCCGCATTTTTCCAACGCCATGTGCAAAAATGCCCTc                                                                                                                                                                                                     >  2:144315/1‑238 (MQ=18)
                                    ccGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCCGCGTGCACCCCTTATCCGGCCTACGTAGAGCCCTGAACTCGTAGGCCTGGTAAGCGTAGCGCCTCAGGCAATTCCCGCCGCAGACCTGTGTCAGCGG‑CTACCGGGATTCATTTCCGCCAACAACCCCCCATTTATCCAACGCCATGTGCAAAAATGCCTTc                                                                                                                                                                                                     >  1:180845/1‑238 (MQ=255)
                                        gcATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCgc‑ggc                                                                                                                                                                                               <  1:570322/239‑1 (MQ=255)
                                                   aCGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGcc                                                                                                                                                                                    <  2:847998/239‑1 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCCGCCGCAGACCTGGGTCAGCGG‑CTACCGTGATTCATTTCCGCCAAACACCCCGCATTTTTCCCACGCCATCTGCAAACAAGCCTTCGC‑CGCGGCCTTCTGCCCG‑CTa                                                                                                                                                                              >  2:490813/1‑239 (MQ=17)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTa                                                                                                                                                                              >  1:1224910/1‑239 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTa                                                                                                                                                                              >  1:1086005/1‑239 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGGTATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGGGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTa                                                                                                                                                                              >  1:957857/1‑239 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCCGCCGTCTGCCAga                                                                                                                                                                                 >  2:1148276/1‑236 (MQ=255)
                                                              aaaaaTGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCCTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCCCCAAACACCACGCATTTATCCCAACCCATGTGCCAAAATGACTTCCC‑GGCGGCTGGCCTCCCT‑CTAtttttt                                                                                                                                                                        >  2:858753/1‑238 (MQ=17)
                                                                     ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCCG‑CTATTTTTCCGCCGca                                                                                                                                                                 >  1:339700/1‑239 (MQ=255)
                                                                     ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAAGCCTTCGC‑GGCGGCTGTCTGCCAG‑CCATTTTTCctccgcc                                                                                                                                                                 >  2:1008839/1‑238 (MQ=255)
                                                                          gCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTTCGTAGAGCACTGCACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCCTGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CCATTTTTCCGCCGCCACaaa                                                                                                                                                            >  1:927026/1‑239 (MQ=255)
                                                                          gCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGCGCTTCATTTCCGCCAAAAACCGCGCATTTCTCCAACGCCATGTGACAAAATGCCTTCCC‑CCGGGCTGGCTGCCCG‑CTCTTTTTCCCGCGCAACcaa                                                                                                                                                            >  2:1356996/1‑239 (MQ=14)
                                                                                                      tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTATTTTTCCGCCGCAACAAAACCGCCGGTCTCTCCAGTAGTGGCgg‑gg                                                                                                                               >  1:790784/1‑239 (MQ=255)
                                                                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGCgg‑ggggacg                                                                                                                          >  1:595315/1‑234 (MQ=255)
                                                                                                            aTGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGGTCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTCTCCCACGCCCTTTGCAAAAATGCCTTCGC‑CGCGGCCGTCTGCCAG‑CTCTTTTTTCTCCCCACCACCACCGCAGCtccctctccacgccgcgacgataac                                                                                                                           >  2:778285/1‑219 (MQ=1)
                                                                                                                     gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTTCCCG‑CTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGAGGCgg‑ggggcgcgtcaacgatt                                                                                                                >  2:973789/1‑224 (MQ=18)
                                                                                                                                        cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACCACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGCGG‑CGCTAGGGACAAAGAttttaaccaccacttgtttgt                                                                                             >  1:1243284/1‑221 (MQ=18)
                                                                                                                                        cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGGAGTGGCGG‑GTCACGCGAGATaacgatgaatacctccagttgttt                                                                                             >  2:745385/1‑216 (MQ=18)
                                                                                                                                             agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGcc                                                                                                                                                                                    <  1:207971/149‑1 (MQ=255)
                                                                                                                                             agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGcc                                                                                                                                                                                    >  2:207971/1‑149 (MQ=255)
                                                                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGCGGCAAGGGAAAAAGCATTAAACCCGCCTcgttgttgtgacacaa                                                                                     >  1:335984/1‑228 (MQ=18)
                                                                                                                                                              aaaaaaaggcgcgccccaccAGACAATCCCCCCCGCACCCCGGTGCCGGCGG‑CCCCCGGTATTTATTCCCCCAAACCCCCGCGTATTTACCAAACCCCATGGAAAAAAACCCCTGCGG‑GGGGTCTGTCT‑CCAGTCTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                                       <  2:279784/220‑1 (MQ=255)
                                                                                                                                                                       ggcagcaattgggcaaaccacAGCCGCAGACCTGTGGCAGCGG‑CTACCGTGTTCCATCCCCGCCAACACCCCCCCATTCCCCCAACGCCATGTAAAAAAATGCTTGCGG‑GGCGGTTGTCTGCCAG‑TTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                                                <  2:868183/217‑1 (MQ=255)
                                                                                                                                                                                                   agccttgtgcAGCGGCCTACGGTTATTCTTTCCCCCCAACACCCCCCCTTTTACCCACCGCCATGTCAAAAATTGCCTCCGC‑GGCGGTTGTCGCCAAG‑TTATTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                                    <  1:1168698/229‑1 (MQ=255)
                                                                                                                                                                                                                            ttCATCCCCGCCACCACCCCCGCTTTTACCCACCGCCTGGTGAAAAAATGCTTCCGG‑GGCGGCTGTCTGCCAC‑CTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa          <  1:778285/239‑1 (MQ=255)
                                                                                                                                                                                                                             tCATTTCCCCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTTCGC‑GGCGCCTTTCTGCCAT‑CTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac         <  1:1008839/239‑1 (MQ=255)
                                                                                                                                                                                                                             cactccccGCCAACACCCGCTCTTTTCCCCAACGACATGTACAAATAGGCCTTCGC‑GGGGCCTCTCTGCCAG‑TTATTTTCCCCCCACAGCAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac         <  1:745385/233‑1 (MQ=255)
                                                                                                                                                                                                                                 cccccccAACAACCGCGCATTTACCCACCGCAGTGTGCAAAAATGCTTCCGG‑GGGGGTTGTCGCCCAG‑TTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGaccacc     <  1:1196506/234‑1 (MQ=255)
                                                                                                                                                                                                                                  ccccccAACACCCCCTCA‑TTATCCAACCCCATGGGAAAAAAACCCCCCGCGGGGGGCTGTCTGCCGG‑TTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCt    <  1:41922/235‑1 (MQ=255)
                                                                                                                                                                                                                                    cGCCAACCACCGCTCATTTACCCACGCCCGTGTGCAAAAATGCCTTGGC‑GGCGTCTGCCTGCCAG‑TTTTTTTCCCGCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCtgt  <  1:858753/239‑1 (MQ=255)
                                                                                                                                                                                                                                    cGCCAACACCCGCCCATTTCCCCACCCCAATGTGCAAAAATGCCTCCGG‑GGCGGCTGTCTGCCAT‑ATTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCtgt  <  1:1148276/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      
GTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGG‑CTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAG‑CTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGT  >  NC_000913/357581‑358046

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