Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 357,819 G→C 19.8% 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,8190GC19.8% 68.0 / 4.0 43A291G (GCG→GGG) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (23/11);  new base C (7/2);  total (30/13)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

GTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGC  >  NC_000913/357581‑358055
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
ttgtCAGCGGCTACCGTTATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    <  1:802019/238‑1 (MQ=255)
gtgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    <  2:342713/239‑1 (MQ=255)
gtgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    <  1:259171/239‑1 (MQ=255)
      gCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCCGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGAGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAtt                                                                                                                                                                                                                                              >  1:658086/1‑239 (MQ=255)
        ggCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                                                                                                                                                               <  1:630245/236‑1 (MQ=255)
         gCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCCTCTCCGGGCTGCGGAGAGCACTGAACTCGTAGGCCTGATAACCGTAG‑CGCATCAGGCAATTCCAGCCGCAAACCTGGGGCAGCGGGGCGCGGGGTTCATTTCCGCCAAACACCCCGCCTTTa                                                                                                                                                                                                                                            >  2:221285/1‑238 (MQ=25)
           tACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGGCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATcc                                                                                                                                                                                                                                         >  1:1164582/1‑239 (MQ=255)
                  attcattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTCTCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCACCCGCAGACCTGTGTCCGCGGCGACCGTGGTTCCTTTCCGCCCACCACCGCCCATTTTTCCACCGcc                                                                                                                                                                                                                                   >  2:1328848/1‑238 (MQ=21)
                                            tcctccAACTCCATGTGCAGAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATCCCCGATACGGCGTTCAACCCTTTTCCCACCGACGTAGAACACCAACCTCGTAGGCGTGAGAAGCGGAGTCGC‑CCAGGCGGTTTGAGGCGCAGAACTGCGGCGGCGGCTGCCGGGATTTATTTCCTCCAACAACCCCGAATTTATACTACTCACT‑GTGCAAAAATCCCCT‑CCCCGCGG‑Ct                                                                                                                                                                                                     >  2:1102823/1‑239 (MQ=255)
                                            tcctccAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACAACCGCGCAATTATCCAACGCCAT‑GTTCAAAAATGCCTT‑Cgcggcg                                                                                                                                                                                                         >  1:729728/1‑236 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGCACGCCTTATCCGGCCTACGTAGAGCACTGAAATCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCCGCCGCAGACCTGTGTCAGCGGCTAACGTGATTCATTTCCGCCAAAAACCCCGCATTTATCCCACGCCAT‑TTGCAAAAATGCCCT‑CGCGGCGG‑CTGGCTCCCCGCTa                                                                                                                                                                                         >  2:952476/1‑239 (MQ=18)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCCTTTATCCAACCCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTa                                                                                                                                                                                         >  1:84942/1‑239 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGcca                                                                                                                                                                                         >  1:1083380/1‑237 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTa                                                                                                                                                                                         >  1:1305601/1‑239 (MQ=255)
                                                         tgtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCCTGATTCATTTCCGCCAACAACCCCCCATTTATCCAACGCCA‑‑‑TGCAAAAATGCCTT‑CGCGGCGG‑CGGTCTGCCACCCAttt                                                                                                                                                                                      >  1:1204456/1‑239 (MQ=18)
                                                         tgtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGAGGCTACCGTGATTCATTTCCGCCAACAACCCCGCATTTATCCCACGCCCT‑GCGCAAAAATGCCTT‑CGCGGCGG‑CTGGCTGCCAGCttg                                                                                                                                                                                        >  2:243771/1‑237 (MQ=255)
                                                          gtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAAAGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTGCCGCCCACAACCGCCCCTTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGACGG‑GCGTCTGCCAGCTAtt                                                                                                                                                                                       >  1:1059123/1‑239 (MQ=17)
                                                                     ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCCCCCATTTATCCAACGCCAT‑GTGCAAAAACGCCTT‑CGCGGCGG‑GCGCCTGCCCGCTCTTTTTtcgccccc                                                                                                                                                                            >  2:1240645/1‑236 (MQ=17)
                                                                     ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGTAGAGCACTGGACCCGTAGGCCTGATAAACGGAG‑CGCATCAGGCAATTCCAGCCGCAGACCTTGGTCAGCGGCTACCGTGATTCATTTCCGCCAACCACCGCGCCTTTAGCCAACGCCAT‑GTGCAACAATGCCTT‑CGCCGCGG‑CTGTCTGCCAGCCATGTTTccgcctc                                                                                                                                                                             >  2:501507/1‑236 (MQ=17)
                                                                              cGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    <  1:1006826/161‑1 (MQ=255)
                                                                              cGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    >  2:1006826/1‑161 (MQ=255)
                                                                                           gCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCCGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACCCCCT‑TGGCAAAAATGCCTT‑CGCGGCGG‑CCGTCTGCCAGCCATTTTTTCCCCCCAACAAAAACGCCCGTCCCTccc                                                                                                                                                       >  2:696177/1‑237 (MQ=14)
                                                                                                aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGcc                                                                                                                                                                                                                                   >  1:1300732/1‑160 (MQ=40)
                                                                                                aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGcc                                                                                                                                                                                                                                   <  2:1300732/160‑1 (MQ=40)
                                                                                                  tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                                               >  1:175765/1‑239 (MQ=255)
                                                                                                  tttATGCCGGATGCGGCGTGAACCCCTTATCCGGGCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGAGGAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAAAGCCTT‑CGCGGCGG‑CTGTCTCGCAGCTATTTTTCCGCCGCAAAAAAACCCCCGGtctctctcgaagttgg                                                                                                                                               >  2:1248628/1‑229 (MQ=17)
                                                                                                      tGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAAGGCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAAACCTGTGTCAGCGGCTCCCGGGGTTCATTTTCGCCCTCAAACGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGACGCGG‑GTGTCTGCCACCTCTTTTTTCGCCGCAACCAAAACCCCGTTGTCTCCCGTCGTGCT‑gggg                                                                                                                                          >  2:1264451/1‑239 (MQ=14)
                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGGACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCCCGCATTTATCCAACGCCAT‑GTGCCAAAATGCCTT‑CGCGGCGG‑CCGTCTTCCCGCTCTTTTTTCCCCGCCACAAAACCGCCGCTCTTTCCCGGAAGGGC‑GGGGc                                                                                                                                         >  2:1227977/1‑239 (MQ=14)
                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCCGCTATTTTTTCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGC‑GGGGc                                                                                                                                         >  2:8340/1‑239 (MQ=255)
                                                                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGC‑gggggctag                                                                                                                                     >  1:128078/1‑234 (MQ=255)
                                                                                                            aTGCGGCGTGAACCCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCCT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCCGTCTCTCCAGTAGTGGC‑gggcggcggg                                                                                                                                    >  1:743392/1‑232 (MQ=255)
                                                                                                                ggcgagcacctcTTACCCGGCCTCCGAAGGGCCCTGAACTCGTAGCCCTGAAAGGCGTGG‑CGCTCCAGCATACCCCCGCCCCGGCCCTGTGTCAGGGGCTACGGTGATTCTTCCCCCCCACCACCCCCGCTTTAACCCACCCCCAT‑GTCAAAAAATCCCTT‑CGGGGCGG‑TTGTCGGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGaaa                                                                                                                                <  2:87675/227‑1 (MQ=255)
                                                                                                                              aTCCGGCCTACGTAGACCTCTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCATCCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAAcccc                                                                                                                                                                                                                                   <  2:1275339/130‑1 (MQ=14)
                                                                                                                                   gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGC‑GGGGCAACGAGAAAA‑AG‑TTTAATccccgctttg                                                                                                           >  1:525869/1‑232 (MQ=18)
                                                                                                                                        cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCaaaa                                                                                                                                                                                                                        >  1:1275339/1‑130 (MQ=255)
                                                                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTg                                                                                                                                                                                                >  2:893224/1‑142 (MQ=255)
                                                                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTg                                                                                                                                                                                                <  1:893224/142‑1 (MQ=255)
                                                                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAA‑AGGAAAA‑AGATTTAA‑CCCCCCATGGTGTTTGGCGATCcct                                                                                            >  1:87052/1‑239 (MQ=18)
                                                                                                                                                     ccTCGTAGGCCTGATAAGCGTAC‑CGCATCAGGCAATCCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                                          <  2:627455/237‑1 (MQ=255)
                                                                                                                                                                        gTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCA‑GAGGAAA‑AACTTTTAAACACTCTTGTTGTGTTGGGAAGGCTGCTGGGAACCACGTgggtatgg                                                                      >  1:1192581/1‑233 (MQ=18)
                                                                                                                                                                              gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCA‑GAGAAATAAGATTTAA‑CCCCTCACTTTGTTTGGGAATAGCTGGTGGGTACAAAGTGGGAAGGGGAGTg                                                                 >  1:470435/1‑239 (MQ=21)
                                                                                                                                                                                                                            tccatccccGCCAACACCCGCACTTTACCCCACCCCCAT‑TTGAAAAAATGCCTT‑GGGGGCGG‑TTGTCGGCCGTCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa                   <  1:286968/232‑1 (MQ=255)
                                                                                                                                                                                                                             tctttccccccAACACCCGCGCTTTTACCCAACGCCAT‑GTCAAAAAATGCCTT‑CGCGGCGG‑TTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                  <  1:1298912/230‑1 (MQ=255)
                                                                                                                                                                                                                             tctttcccGCCAACCCCCGCGCATTT‑CCCAACCCCATGGTCAAAAAATCCCTT‑CGCGCGGG‑CTGTCGGCCAGTTATTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                  <  2:1305601/233‑1 (MQ=255)
                                                                                                                                                                                                                             tctttcccGCAAACACCCGCGCATTTACCCAACCCCAT‑GTAAAAAAATGCCTG‑CGGGGCGG‑CGGTCCGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                  <  1:8340/233‑1 (MQ=255)
                                                                                                                                                                                                                                ttcccGCCAACACCCCCACTTTTACCCACCGCCAT‑GTAAAAAAATGCTTC‑CGCGGCGG‑CTGCCTGCCATTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGaccac               <  1:134787/236‑1 (MQ=255)
                                                                                                                                                                                                                                    cccaaacaccCGCGCTTTTACCCACCGCCTG‑GTGCAAAAAGGCCTC‑CG‑GGCGGTTTTTCGCCCAGTTTTTTTCCCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCtgt           <  1:508315/235‑1 (MQ=255)
                                                                                                                                                                                                                                             ccGCGAATATACCCAA‑GCCAT‑GGTCAAAAACGCCTTCCGGGGCGG‑TTTTCTGCTGGTTTTTTCCCCGCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGc  <  2:1149998/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
GTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGC  >  NC_000913/357581‑358055

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