Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,335	0	T	G	51.0%	50.8 / 72.3	51	G25G (GGT→GGG)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (8/17);  new base G (26/0);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 7.32e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-06

GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTC  >  NC_000913/523188‑523484                                                                                                                                                    |                                                                                                                                                        ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCACTTTCCCGGGGTTGGCAGGGGTa                                                                                                                                                      >  2:932681/1‑151 (MQ=255)    cTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCTGGGTCGGCCGGGGTAAG‑a                                                                                                                                                  >  2:230563/1‑151 (MQ=255)     ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑aa                                                                                                                                                 <  2:397271/151‑1 (MQ=255)     ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑aa                                                                                                                                                 <  1:1114962/151‑1 (MQ=255)      tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTGGGGAGGTGTTAG‑AAt                                                                                                                                                >  2:1097205/1‑151 (MQ=255)           gtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCAGGGGTTGGCCGGGGGAAG‑AAATgggg                                                                                                                                           >  1:270004/1‑149 (MQ=255)            tCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGgcg                                                                                                                                           <  2:579350/150‑1 (MQ=255)                            tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCggcgtggcg                                                                                                                          <  2:406833/151‑1 (MQ=255)                                gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCATTGTCCAGGGGTCGGCAGGGGTAAG‑AATTGGCGCGCCCACCGGCGCGGCGTGc                                                                                                                       >  1:1006658/1‑150 (MQ=255)                                gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAG‑AAATGGCGCGGCCCCCggcggtgcgggc                                                                                                                       >  2:319789/1‑142 (MQ=255)                                     cacaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAG‑AATTGGCGCGCCCACCGGCGGGGCGTGCTCgggg                                                                                                                 >  1:327722/1‑149 (MQ=255)                                         gggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGTGTAAG‑AAATTGCGCGCCCACCGGGCGGGCCTGCTCGgtgtgtt                                                                                                             >  2:914735/1‑150 (MQ=255)                                         gggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAG‑AATTTGCGCGCCCACCGGCGTGGCGTGCTCgggggtgt                                                                                                             >  1:1195459/1‑145 (MQ=255)                                         gggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGGGCTCGgtgtgt                                                                                                              >  1:637216/1‑150 (MQ=255)                                         gggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtg                                                                                                                 <  2:1194971/147‑1 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCATTTTCCAGGGGTCGGCAGGGGGAAG‑AATTGGGGGGCCCACCGGCGTGGGGTGCTCggggggggcg                                                                                                           >  2:493480/1‑143 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCAATTTCCAGGGGTTGGCCGGTGTAAG‑AATTGGCGCGCCCCACGGCGGGGCGTGCTCGGGGTGT‑Ccg                                                                                                          >  1:1023713/1‑151 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCCGGGGTAAA‑AATTGGCGCGCCCACCGGCGGGGCGGGGTCGGGGTGG‑Ccg                                                                                                          >  1:268312/1‑151 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGGAAG‑AATTTGCGCGCCCCCCGGCGTGGCGTGCTCGGTGGGT‑Tcg                                                                                                          >  1:147482/1‑151 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAG‑AATTTGCGCGCCCACCGGCGTGGCGGGCCCggggggttcg                                                                                                           >  1:826406/1‑143 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTACCCTTGTCCACGGTTCCGCAGGGGGAAG‑AATTTTGGGGCCCACCGGCGGGGCGGGCTCCGGGGGT‑Ccg                                                                                                          >  2:1054841/1‑151 (MQ=255)                                           gATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGTCCATTGTCCCGGGGTCGGCAGGGGTAAG‑AAATTGCGCCCCCCCCCGCGTGGCCCGCCCGGGGGGT‑Ccg                                                                                                          >  2:37498/1‑151 (MQ=255)                                                gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAG‑AAATTGCGCGCCCCCCGGCGGGGGGGGCTCGGGGTGG‑Ccgggggg                                                                                                     >  1:48961/1‑151 (MQ=255)                                                 aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAG‑AATTGGCGCGCCCACCCGCGTGGCGGGCCCGGGGGGT‑Ccgggcggg                                                                                                    >  1:228173/1‑151 (MQ=255)                                                     aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTGAG‑AAATGGCGCGCCCCCCCGCGTGGCGGGCTCGGTGGGTGCcgggggggtgg                                                                                                 >  2:1005430/1‑148 (MQ=255)                                                         cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAG‑AATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGGGTGT‑Gcgggcg                                                                                                      >  2:10909/1‑141 (MQ=255)                                                         cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑Ccgggcg                                                                                                      <  1:10912/141‑1 (MQ=255)                                                             aaGGTCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGc                                                                                        <  2:1023341/151‑1 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCCTCTTCACGGTTCGGGCGGGGGAAGCAATTGG‑GTGCCCACCCGCGTTGCCTGCCCGGGGTGT‑GCGGGGTGGGGGCCGTCGGGCAACACgg                                                                                 >  2:706273/1‑151 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGGAAG‑AATTGGCGCGCCCACCGGCGTGGGGGGCTCCGGGTGT‑CCGGGCGGGGGGACTTCGGGGCAccccg                                                                                 >  2:1008246/1‑149 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAG‑GAATGGGGCGCCCACCGGCGGGGCGTGCTCGGTTGGT‑CCGGGCGGGGTGGCTTCTGGCAAccccg                                                                                 >  2:187935/1‑149 (MQ=255)                                                                    cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGGTCGGGGTGT‑CCGGGCGGGGTGGCTTCGGGCAACCCgg                                                                                 >  1:144482/1‑151 (MQ=255)                                                                      ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGTAAG‑AATTGGGGCGCCCACCGGCGTGGCGGGGTCGGGGTGT‑CCGGGGGGGGTGACTTCGGGCAACCCggga                                                                               >  1:318514/1‑149 (MQ=255)                                                                       tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                              <  2:790232/151‑1 (MQ=255)                                                                       tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                              <  1:1126432/151‑1 (MQ=255)                                                                          atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATc                                                                           <  2:716752/151‑1 (MQ=255)                                                                                          cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTAGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCg                                                             <  1:983372/149‑1 (MQ=255)                                                                                             cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGt                                                        <  1:614285/151‑1 (MQ=255)                                                                                                 tCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGTTCCGCAGGGGTAAG‑AATTGGCGCGCCCCCCGGCGTGGGGTGGGCGGTTTGT‑TCCGGCGGGGTGGCTTCGGGCAACCCGGGAAATCCCCTGCTGGGGGGGGAgggg                                                       >  2:29077/1‑146 (MQ=255)                                                                                                 tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGAGGGGGAAG‑AATTGGCGCGCCCCCCGGCGTGGCGGGCTCGGGGTGT‑CCGGGCGGGGTGACTTCGGGCAACCCCGGAAATCCGCTGCCGGGGGCGGAgggg                                                       >  1:279667/1‑146 (MQ=255)                                                                                                        gaTATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCCGGTGTAAG‑AATTTGCGCGCCCCCCGGCGTGGCGGGCTCGGGGTGG‑TCGGGGGGGGTGACTTCGGGCCACCCGGGAAAACCCCTGCTGGGGGCGGAGGGGCGGCCCcgcg                                             >  2:537124/1‑151 (MQ=255)                                                                                                        gaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAA‑AATTGGCGGGCCCACCGGCGGGGCGTGCTCGGTGGGT‑CCGGGCGGGATGACTTCGGGCAACCCGGGAAATCCGCCGCTGGGGGGGAAGGGGCTGCCCGGCg                                             >  1:1000781/1‑151 (MQ=255)                                                                                                                tCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAc                                     <  2:332656/151‑1 (MQ=255)                                                                                                                        gCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcgc                              <  1:432584/150‑1 (MQ=255)                                                                                                                        gCGGTCCCATTGTCCAGGGCTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcgc                              <  1:1267588/150‑1 (MQ=255)                                                                                                                               cATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccgg                       >  1:1059418/1‑150 (MQ=255)                                                                                                                                 ttGTCCAGTGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc                    <  2:826080/151‑1 (MQ=255)                                                                                                                                   gTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGCGACGGACCTTGCGCTgcccggcccg                   >  1:1231212/1‑150 (MQ=255)                                                                                                                                        gggTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGGGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCg             >  2:486449/1‑151 (MQ=255)                                                                                                                                          gTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGtt           >  2:945503/1‑151 (MQ=255)                                                                                                                                                cAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtc      <  1:638245/150‑1 (MQ=255)                                                                                                                                                   gtgtAAG‑CATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGCAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGCCCTTGCGCTGCCCGGCCCGCTGCCGGTCCTtcccc   >  2:164265/1‑147 (MQ=255)                                                                                                                                                   gtgtAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct   >  1:957043/1‑150 (MQ=255)                                                                                                                                                   gtgtAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct   >  1:84558/1‑150 (MQ=255)                                                                                                                                                   gtgtAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctct   >  1:100590/1‑150 (MQ=255)                                                                                                                                                    tgtAAG‑AATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctctc  <  1:673391/150‑1 (MQ=255)                                                                                                                                                    tgtAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctctc  <  1:224202/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAG‑AATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGT‑CCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTC  >  NC_000913/523188‑523484

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