Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,571,096	A→G	26.4%	L160P (CTC→CCC)	yjiR ←	putative DNA‑binding transcriptional regulator/putative aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,571,096	0	A	G	26.4%	84.3 / 4.5	45	L160P (CTC→CCC)	yjiR	putative DNA‑binding transcriptional regulator/putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (8/24);  new base G (10/2);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCG  >  NC_000913/4570874‑4571321                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        aGTTCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGAACACCCAG‑CCCCCGGTCAGGATTAGGCG‑GGcga                                                                                                                                                                                                                         >  1:667984/1‑236 (MQ=255)    tCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCa                                                                                                                                                                                                                          <  1:517607/234‑1 (MQ=255)                    tGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCAG‑CCGCCGTCGGGCCTCAGGCG‑GGGAATCCGCTCTCCCAATACTCg                                                                                                                                                                                                      >  1:1010473/1‑238 (MQ=255)                      aTGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGGGACCACCCAGTCCCCCGT‑GTGCCTCAGGCGCGGCAA‑CTGCTCTTCCAATACTctacg                                                                                                                                                                                                   >  1:923598/1‑235 (MQ=255)                        gccgtttgcTGGATCGGTTGGAATTCCAATCCCTTTCACGCCCATGCCGCGCAGCATCGGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACCCACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGAc                                                                                                                                                                                                <  2:413244/231‑1 (MQ=255)                                 tGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCAG‑CCGCCGTCGAGCATCAGGCG‑GGGAAACTGTGCTCGCCAAAATCGC‑TGCCCGGCTaaaa                                                                                                                                                                                       >  1:552895/1‑237 (MQ=255)                                 aggATGGTTTGGAATTTCAATCCCTTTCCCCCCCAGCCCGCCCAGCATCTGCATCGAACCTTAATAACAGGGGTATTCGCCCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTTATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACa                                                                                                                                                                                       <  2:650137/238‑1 (MQ=255)                                             attttCAATCTCTTTCACCCCAATGCCGCGCACCATCTGCACCGAACCGTATTACCAGGGGGATTCGCCCGCG‑ACAATATCCCCCGGTTTACACACGGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTaa                                                                                                                                                                                         <  2:1032461/223‑1 (MQ=255)                                               tttCAATCACTTTCACGCCCAGGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTACCGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCgaga                                                                                                                                                                         <  1:615755/239‑1 (MQ=255)                                                       aCTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGGGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCCC‑TGACCGGGTCACAAGGCATaa                                                                                                                                                                               >  1:1032461/1‑225 (MQ=255)                                                           tCACGCCCATGCCGCGCAGCATCTGCATCGAACCTGAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTACCGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑a                                                                                                                                                            <  1:837146/239‑1 (MQ=255)                                                                                    cATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTaa                                                                                                                                                                                         <  1:348798/186‑1 (MQ=255)                                                                                              taataaCGGGGGGTTTCGCCCCGA‑AAAATTCCCCCCGGTTCACCCCCTGCCATTACCGCCAGGAAACGCGATTTAGGGCGCCCGCTGGTGATGAGGATGTCATCGGGGGTGCCCCCCGGG‑CCGCTGTCGAGCATCAGGGG‑GGCAATCTGCTCTCGCAATACTCGC‑TGCCCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAg                                                                                                                         <  2:856444/239‑1 (MQ=255)                                                                                              taataaCAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGGGGACACCCAG‑CCGCCGGCGGGGAACAGGGG‑GGGCAACTGCGCCCCCAATACTCGC‑TGACCGGGGTACAAGACATaa                                                                                                                                                                               >  2:348798/1‑186 (MQ=255)                                                                                               aataaCAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTAtg                                                                                                                                                       <  1:652594/208‑1 (MQ=255)                                                                                                      gggggATTCGACCGCG‑ACAATTTCGCCCGTTTTACACCCTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGct                                                                                                                                                                                                                   <  2:718574/143‑1 (MQ=255)                                                                                                        gtgtattgACCGCGCACAATACCCCCTGTTTTAACCCCGGCCATTACCGCCGCGGAGATCGATTTATGCCAGCCGCTGGTGATTGTGATGCCACCGGCGGTGCCCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCGCGCAATACTCGC‑TGCCCGGCTAACAAGTCATAACCGAGAACGGTTGGCAG‑ATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGC‑G‑Gtt                                                                                                                  <  2:1214074/228‑1 (MQ=255)                                                                                                         ggATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCCAG‑CCCCCCTCCGGCCTTCCGCG‑GGCGACCTGCTCTCCCCATACCCCCGCTACCGGC‑AACAAATCATAACCCGGAAAAGGTTGCAG‑AGTATATAGCACCACCCCGCGCATCTCCCGCCACAAC‑GCGTTTcggc                                                                                                             >  2:578844/1‑236 (MQ=255)                                                                                                                aCCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGGGGCCCCCCAG‑CCGCCGTTGGGCATCCGGGG‑GGGAATCTCCTCTCCCacata                                                                                                                                                                                                         >  1:718574/1‑139 (MQ=255)                                                                                                                        cAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCg                                                                                                                                                                                                               <  1:77533/130‑1 (MQ=255)                                                                                                                                  cGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCAG‑CCGCCGTCGGGCATCAAGCG‑GGCAATCCGCTCTCGCAATAATCGc                                                                                                                                                                                                     >  2:77533/1‑130 (MQ=255)                                                                                                                                     tttACACCTGGCAATTAACGCCAGCGACTTCGAGTTATGGCAGCCGCTGGGGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGgt                                                                                        <  1:1017262/232‑1 (MQ=255)                                                                                                                                     tttACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAAc                                                                                                                                                                       <  1:545923/156‑1 (MQ=255)                                                                                                                                     tttACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCt                                                                                                                                        <  1:484493/186‑1 (MQ=255)                                                                                                                                               gCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCCAG‑CCGCTGTCGGGCATCAGGCG‑GGCAATCTTCTCTCGCAATACTCGC‑TGTCCGGCGAACAAAGCATAACCCAGAACGGGTTGCAG‑AGTATTCTGCACCCCCCCGCGTAGCTCCCGc                                                                                                                              >  2:484493/1‑186 (MQ=255)                                                                                                                                               gCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCCG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTTCTCTCCCCATAACCCC‑TTACCGGCTCACAAGGCATAACCCGGAACGGGTTGTAG‑GTTTTGCTTCCCCCCCCGGCTTTGCTCCCGCCCCACC‑G‑GGGTTACGCCGTGCGCGTTGTCATCCCGCGGGCtttttttctttgt                                                                       >  2:1017262/1‑230 (MQ=255)                                                                                                                                               gCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCAG‑CCCCCGTCGGGCATCAAGCG‑GGCAAATTGCCCTCCCAAAACATCC‑TGACCCGCTTACAAATCCTAACCGCGAAAGGgttttaga                                                                                                                                                             >  2:545923/1‑151 (MQ=255)                                                                                                                                               gCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGGGGCCACCCAG‑CCGCCGTCGAGCa                                                                                                                                                                                                                                      >  2:220331/1‑85 (MQ=255)                                                                                                                                                ccATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑a                                                                                                                                                            <  1:599266/155‑1 (MQ=255)                                                                                                                                                         ccAGCGACATGGGTTTAGGGCAGCCGGTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGGGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGaa                                                            <  1:578844/239‑1 (MQ=255)                                                                                                                                                              gACATCGATTTATGGCGGCCGCGGGTGATGATGATGTCATCGGCGGTGCCCACCGAG‑CCGCTGGCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATCCTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAAc                                                           <  1:1163859/235‑1 (MQ=255)                                                                                                                                                                  tCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCAG‑CCCCCGTCGAGCATCAGGGG‑GGCAATCTGCTCTCGCCATACTCCC‑TGACCCGCCAACAAATCATAAACCAGAAACGGTTGCAGAATTAGTCT‑CACCCCCCCGCCTTACTCCCGCCCCCGCGG‑GTTCCAGGTGGGGGCCTCCGCCACCCGGGGGGCTTTT‑CTTTACCGAGCAATAGGAAC‑GTCGcca                                                   >  2:1155879/1‑237 (MQ=255)                                                                                                                                                                      gTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCa                                                                                                                                                                                  <  1:103294/112‑1 (MQ=255)                                                                                                                                                                      gTTATGCCAGCAGTTGGTAAGGATCATTCCATGGGCGGTCACCCCCGAC‑CCGCTGCCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATCCTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATgc                                                <  2:456570/238‑1 (MQ=255)                                                                                                                                                                         agggCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATGCGCCt                                            <  1:486176/237‑1 (MQ=255)                                                                                                                                                                                ccGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCCGTCGGGCAACAGGGG‑GGGAATCTGCTCTCGCAAAACTCGC‑TGACCGGCTAACAAGGCATAAACGAGaa                                                                                                                                                                        >  2:103294/1‑112 (MQ=255)                                                                                                                                                                                 cGCTGGTGATGATGATGTCATCGGCGGTGCCCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATGCGCCTCCa                                         <  1:607022/234‑1 (MQ=255)                                                                                                                                                                                    tGGTGATGATGATGTCATCGGCGGTGCCCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCCCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTg                                                                               <  1:1133774/194‑1 (MQ=255)                                                                                                                                                                                       tgattatgaTGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATGCGCCTCCAGCATATCCAg                               <  1:1130556/238‑1 (MQ=255)                                                                                                                                                                                              gatgTCATCGGCGGTGACCACCGAG‑CCGCTGTCGGGCATCAGGCG‑GGCAATCTTCTCTCGCAATAATCGC‑TGACCGGCTAACAAATCATAAACGAGAAAGGGTTGCAG‑ATTATGCTGCACCACCCCGCTTAACTCCCCCCCCAGC‑G‑GGTTCAGGCTGGGCGCTTCGACATCCGGGGTGCTTTTGCTTTACggga                                                                     >  2:1133774/1‑192 (MQ=255)                                                                                                                                                                                                                  cgtgg‑CCGCTGCCGGGCATCAGGCG‑GCCAATCTGCCCACGCAATCCTCGT‑TGCCCGGCTAACAAGTCATAACCGAAACGGGTTGGCAG‑ATTAGGCTGCCCCCCCCGGCTTAGCCCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCAcc   <  2:477706/235‑1 (MQ=255)                                                                                                                                                                                                                   cGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATGCGCCTCCAGCCAATCCAGAACCTGATCCCCCTGGGGAATTTCCACCg  >  2:573304/1‑239 (MQ=255)                                                                                                                                                                                                                      g‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCa     >  1:1201448/1‑233 (MQ=255)                                                                                                                                                                                                                           cTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑a                                                                                                                                                            <  2:601689/81‑1 (MQ=255)                                                                                                                                                                                                                           cTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑a                                                                                                                                                            >  1:601689/1‑81 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        AGTTCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCG‑ACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAG‑CCGCTGTCGAGCATCAGGCG‑GGCAATCTGCTCTCGCAATACTCGC‑TGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAG‑ATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGC‑G‑GTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAAT‑GGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCG  >  NC_000913/4570874‑4571321

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