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

GTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGC  >  NC_000913/357581‑358055
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
ttgtCAGCGGCTACCGTTATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    <  1:802019/238‑1 (MQ=255)
gtgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    <  2:342713/239‑1 (MQ=255)
gtgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    <  1:259171/239‑1 (MQ=255)
      gCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCCGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGAGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAtt                                                                                                                                                                                                                                              >  1:658086/1‑239 (MQ=255)
        ggCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                                                                                                                                                               <  1:630245/236‑1 (MQ=255)
         gCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCCTCTCCGGGCTGCGGAGAGCACTGAACTCGTAGGCCTGATAACCGTAGCGCATCAGGCAATTCCAGCCGCAAACCTGGGGCAGCGGGGCGCGGGGTTCATTTCCGCCAAACACCCCGCCTTTa                                                                                                                                                                                                                                            >  2:221285/1‑238 (MQ=25)
           tACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGGCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATcc                                                                                                                                                                                                                                         >  1:1164582/1‑239 (MQ=255)
                  attcattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTCTCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCACCCGCAGACCTGTGTCCGCGGCGACCGTGGTTCCTTTCCGCCCACCACCGCCCATTTTTCCACCGcc                                                                                                                                                                                                                                   >  2:1328848/1‑238 (MQ=21)
                                            tcctccAACTCCATGTGCAGAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATCCCCGATACGGCGTTCAACCCTTTTCCCACCGACGTAGAACACCAACCTCGTAGGCGTGAGAAGCGGAGTCGCCCAGGCGGTTTGAGGCGCAGAACTGCGGCGGCGGCTGCCGGGATTTATTTCCTCCAACAACCCCGAATTTATACTACTCACTGTGCAAAAATCCCCTCCCCGCGGCt                                                                                                                                                                                                     >  2:1102823/1‑239 (MQ=255)
                                            tcctccAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACAACCGCGCAATTATCCAACGCCATGTTCAAAAATGCCTTCgcggcg                                                                                                                                                                                                         >  1:729728/1‑236 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGCACGCCTTATCCGGCCTACGTAGAGCACTGAAATCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCCGCCGCAGACCTGTGTCAGCGGCTAACGTGATTCATTTCCGCCAAAAACCCCGCATTTATCCCACGCCATTTGCAAAAATGCCCTCGCGGCGGCTGGCTCCCCGCTa                                                                                                                                                                                         >  2:952476/1‑239 (MQ=18)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCCTTTATCCAACCCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                                                                                                                                         >  1:84942/1‑239 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGcca                                                                                                                                                                                         >  1:1083380/1‑237 (MQ=255)
                                                        aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                                                                                                                                         >  1:1305601/1‑239 (MQ=255)
                                                         tgtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCCTGATTCATTTCCGCCAACAACCCCCCATTTATCCAACGCCATGCAAAAATGCCTTCGCGGCGGCGGTCTGCCACCCAttt                                                                                                                                                                                      >  1:1204456/1‑239 (MQ=18)
                                                         tgtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGAGGCTACCGTGATTCATTTCCGCCAACAACCCCGCATTTATCCCACGCCCTGCGCAAAAATGCCTTCGCGGCGGCTGGCTGCCAGCttg                                                                                                                                                                                        >  2:243771/1‑237 (MQ=255)
                                                          gtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAAAGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTGCCGCCCACAACCGCCCCTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGACGGGCGTCTGCCAGCTAtt                                                                                                                                                                                       >  1:1059123/1‑239 (MQ=17)
                                                                     ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCCCCCATTTATCCAACGCCATGTGCAAAAACGCCTTCGCGGCGGGCGCCTGCCCGCTCTTTTTtcgccccc                                                                                                                                                                            >  2:1240645/1‑236 (MQ=17)
                                                                     ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGTAGAGCACTGGACCCGTAGGCCTGATAAACGGAGCGCATCAGGCAATTCCAGCCGCAGACCTTGGTCAGCGGCTACCGTGATTCATTTCCGCCAACCACCGCGCCTTTAGCCAACGCCATGTGCAACAATGCCTTCGCCGCGGCTGTCTGCCAGCCATGTTTccgcctc                                                                                                                                                                             >  2:501507/1‑236 (MQ=17)
                                                                              cGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    <  1:1006826/161‑1 (MQ=255)
                                                                              cGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                                                                                                                                                    >  2:1006826/1‑161 (MQ=255)
                                                                                           gCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCCGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACCCCCTTGGCAAAAATGCCTTCGCGGCGGCCGTCTGCCAGCCATTTTTTCCCCCCAACAAAAACGCCCGTCCCTccc                                                                                                                                                       >  2:696177/1‑237 (MQ=14)
                                                                                                aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGcc                                                                                                                                                                                                                                   >  1:1300732/1‑160 (MQ=40)
                                                                                                aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGcc                                                                                                                                                                                                                                   <  2:1300732/160‑1 (MQ=40)
                                                                                                  tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                                               >  1:175765/1‑239 (MQ=255)
                                                                                                  tttATGCCGGATGCGGCGTGAACCCCTTATCCGGGCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGAGGAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAAAGCCTTCGCGGCGGCTGTCTCGCAGCTATTTTTCCGCCGCAAAAAAACCCCCGGtctctctcgaagttgg                                                                                                                                               >  2:1248628/1‑229 (MQ=17)
                                                                                                      tGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAAGGCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAAACCTGTGTCAGCGGCTCCCGGGGTTCATTTTCGCCCTCAAACGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGACGCGGGTGTCTGCCACCTCTTTTTTCGCCGCAACCAAAACCCCGTTGTCTCCCGTCGTGCTgggg                                                                                                                                          >  2:1264451/1‑239 (MQ=14)
                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGGACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCCCGCATTTATCCAACGCCATGTGCCAAAATGCCTTCGCGGCGGCCGTCTTCCCGCTCTTTTTTCCCCGCCACAAAACCGCCGCTCTTTCCCGGAAGGGCGGGGc                                                                                                                                         >  2:1227977/1‑239 (MQ=14)
                                                                                                       gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCCGCTATTTTTTCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGCGGGGc                                                                                                                                         >  2:8340/1‑239 (MQ=255)
                                                                                                           gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGCgggggctag                                                                                                                                     >  1:128078/1‑234 (MQ=255)
                                                                                                            aTGCGGCGTGAACCCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCCTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCCGTCTCTCCAGTAGTGGCgggcggcggg                                                                                                                                    >  1:743392/1‑232 (MQ=255)
                                                                                                                ggcgagcacctcTTACCCGGCCTCCGAAGGGCCCTGAACTCGTAGCCCTGAAAGGCGTGGCGCTCCAGCATACCCCCGCCCCGGCCCTGTGTCAGGGGCTACGGTGATTCTTCCCCCCCACCACCCCCGCTTTAACCCACCCCCATGTCAAAAAATCCCTTCGGGGCGGTTGTCGGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                                                                <  2:87675/227‑1 (MQ=255)
                                                                                                                              aTCCGGCCTACGTAGACCTCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCATCCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAAcccc                                                                                                                                                                                                                                   <  2:1275339/130‑1 (MQ=14)
                                                                                                                                   gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGCGGGGCAACGAGAAAAAGTTTAATccccgctttg                                                                                                           >  1:525869/1‑232 (MQ=18)
                                                                                                                                        cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCaaaa                                                                                                                                                                                                                        >  1:1275339/1‑130 (MQ=255)
                                                                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTg                                                                                                                                                                                                >  2:893224/1‑142 (MQ=255)
                                                                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTg                                                                                                                                                                                                <  1:893224/142‑1 (MQ=255)
                                                                                                                                                  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGGAAAAAGATTTAACCCCCCATGGTGTTTGGCGATCcct                                                                                            >  1:87052/1‑239 (MQ=18)
                                                                                                                                                     ccTCGTAGGCCTGATAAGCGTACCGCATCAGGCAATCCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                                          <  2:627455/237‑1 (MQ=255)
                                                                                                                                                                        gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGGAAAAACTTTTAAACACTCTTGTTGTGTTGGGAAGGCTGCTGGGAACCACGTgggtatgg                                                                      >  1:1192581/1‑233 (MQ=18)
                                                                                                                                                                              gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAATAAGATTTAACCCCTCACTTTGTTTGGGAATAGCTGGTGGGTACAAAGTGGGAAGGGGAGTg                                                                 >  1:470435/1‑239 (MQ=21)
                                                                                                                                                                                                                            tccatccccGCCAACACCCGCACTTTACCCCACCCCCATTTGAAAAAATGCCTTGGGGGCGGTTGTCGGCCGTCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa                   <  1:286968/232‑1 (MQ=255)
                                                                                                                                                                                                                             tctttccccccAACACCCGCGCTTTTACCCAACGCCATGTCAAAAAATGCCTTCGCGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                  <  1:1298912/230‑1 (MQ=255)
                                                                                                                                                                                                                             tctttcccGCCAACCCCCGCGCATTTCCCAACCCCATGGTCAAAAAATCCCTTCGCGCGGGCTGTCGGCCAGTTATTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                  <  2:1305601/233‑1 (MQ=255)
                                                                                                                                                                                                                             tctttcccGCAAACACCCGCGCATTTACCCAACCCCATGTAAAAAAATGCCTGCGGGGCGGCGGTCCGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                  <  1:8340/233‑1 (MQ=255)
                                                                                                                                                                                                                                ttcccGCCAACACCCCCACTTTTACCCACCGCCATGTAAAAAAATGCTTCCGCGGCGGCTGCCTGCCATTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGaccac               <  1:134787/236‑1 (MQ=255)
                                                                                                                                                                                                                                    cccaaacaccCGCGCTTTTACCCACCGCCTGGTGCAAAAAGGCCTCCGGGCGGTTTTTCGCCCAGTTTTTTTCCCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCtgt           <  1:508315/235‑1 (MQ=255)
                                                                                                                                                                                                                                             ccGCGAATATACCCAAGCCATGGTCAAAAACGCCTTCCGGGGCGGTTTTCTGCTGGTTTTTTCCCCGCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGc  <  2:1149998/239‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
GTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGC  >  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: