BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,068,572	G→C	16.0%	A423G (GCG→GGG)	rutG ←	pyrimidine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,068,572	0	G	C	16.0%	134.6 / 2.7	61	A423G (GCG→GGG)	rutG	pyrimidine permease
Reads supporting (aligned to +/- strand): ref base G (22/28); new base C (8/3); total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

TGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAG  >  NC_000913/1068338‑1068810

gaaagaaTGAAAGCTTCCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATA‑TATTGAATTTGCCCTGTTTTGTAGGCCGGATAAGGCGTTT‑CCGCCGCAC‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCa                                                                                                                                                                                                                                                <  1:668970/232‑1 (MQ=255)

           aGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGCTTTCTGTTTATGGCTCCTGATGAACAACTTCTTGCGGGGGAAAGCCAACCCACTTTCTGCTTAAAACCCCATTGAGTAAAATc                                                                                                                                                                                                                                    >  1:656113/1‑238 (MQ=37)

               tttaataggatacaggCTTACCGGGCCTACCAGAAGATTG‑AATATTATTGATTTGGCCTTGTTTTGTGGGCCGGAAAAGGTTTTG‑CCGCCCCAC‑CCGGCATAAAAAATGCTTCCGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGcc                                                                                                                                                                                                                               <  1:1035949/224‑1 (MQ=255)

                gCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCCCCTGTTGAACAACTTTTGGCGGTGGGACGTCAACCCAACTTCTGCTTTACAAACCAttttgaaaaaactttc                                                                                                                                                                                                                                >  2:668970/1‑223 (MQ=25)

                 ccTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGTTGAACAACTTCTTGTGGTGGAACGTCAACCAACTTTTTGCTTACCAACGCCTTGGGTAAAATcagcccaa                                                                                                                                                                                                                             >  1:845844/1‑232 (MQ=37)

                 ccTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGGTTTTTGTTTATGGCTCCTGGTGAACAACCTCTGGCCGGGGAACCTCAACCAACTTTTTCCTTAAAACAGCATTTGGTAAAAACGCCCCaa                                                                                                                                                                                                                             >  2:1064744/1‑239 (MQ=25)

                  ctgtggggtgttttttATCGGGCCAACAAGAAGTTTGAAATT‑TTTTGATTTTCCCTTGTTTGTGAGGCGGAAAAAGGTTTTT‑CCGCCGCAG‑CGGGCATGAAAAATGC‑TGTGGTGTCCAACAACCGGCCCCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCaaa                                                                                                                                                                                                                            <  2:8476/225‑1 (MQ=255)

                                  tcaggtctcccaagagtaTTGCATTT‑TATTGAGTTGGTCTTGTTTTGTGGGGCGGATAGGTGGTTT‑CCCCCGGCT‑CTGGAATAAACAGTGTGTCCGTTTA‑CAAAATTCGGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AAtc                                                                                                                                                                                                           <  2:656113/222‑1 (MQ=255)

                                  aTCAGGCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGCCAACCAACTTTCTGCTTAACAACCCCTTTGGTTAAATctgtttcttttaaccatttcccgg                                                                                                                                                                                                             >  1:838045/1‑215 (MQ=37)

                                       gCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTTATGAACAACTTCTCGCGGTGGAACGTCAACCAAATTTCTGCTTAACAACGCATTGTTGAAAATCGCGCCGAAGGGTCCTGTccctttctcttt                                                                                                                                                                                                       >  1:377030/1‑230 (MQ=37)

                                       gCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCCCCTGATGAACAACTTCTGGCGGTGGAACGGCAACCAACTTTCTTCTTAACAACCCATTTAGTAAAAACGCGCCaaaaggttcgtttacaaaaccccc                                                                                                                                                                                                       >  2:1004357/1‑218 (MQ=37)

                                         cTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAAAAACTTCTGGCGGGGGAACGTCAAACAACTTTTCTCTTAAAAAACCATTTAGGAAAAAACCGCCCAAggggtcttctgtttcttaaaaac                                                                                                                                                                                                     >  2:1035949/1‑218 (MQ=25)

                                                    aggttgttaaaggggacgggctcTGTTTGGTAGGCCGAATAAGGCTTTT‑ACGCCGCAT‑CCGGAATGAACAATGCTTACGTTGC‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAAc                                                                                                                                                                                           <  1:186514/216‑1 (MQ=255)

                                                             tatTGAATTGGCACTGTTTGTGAGGCCGGAAAGGGCGTTT‑ACGCCGCAT‑CGGGCAAGAACAATGCGTACGTTGA‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGt                                                                                                                                                                                 <  1:749807/239‑1 (MQ=255)

                                                                  aaTTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTTGCGGTGGAACGTCAACCAACTTTCTGCTTAAAAACCCATTGAGTAAAAACGCGCCAAAAGTTGCTGTACC‑AACCCCCCCCCACCAaaaaaaccccacccggaaag                                                                                                                                                                            >  2:1052171/1‑223 (MQ=38)

                                                                      tgccttGTTTGGTAGGCCGGATAGGGCGTTT‑ACGCCGCAT‑CCGGAATGAAGAAGGCGTCCGTTGC‑CAAAAACCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCaaa                                                                                                                                                                        <  1:1064744/234‑1 (MQ=255)

                                                                       gCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTTGCGGTGGAACGTCAACCAACTTTCTGCTTTACAACGCATTAAGTAAAATAGCGCCAAAGGGTTCTTTCCC‑AATCCCCCCCAACGAAAAAACCCCCCACGTgaaggagaaa                                                                                                                                                                       >  2:842600/1‑231 (MQ=38)

                                                                       gCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAAAAAATTCTGGCGGTGGAAAGTAAACCAACTTTCTTCTTTACAAAGCATTGAGTAAAAACCGGCCAAAGTGTGCTTTACC‑AAGCCCACCCAACGCAAAAAAACCCACAGAGAGAGAaaaa                                                                                                                                                                       >  2:434562/1‑239 (MQ=37)

                                                                       gCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCCGATGAAAAACTTCTGGCGGGGGGACGTTAAACAAATTTCTGCTTAACAAAGAATTTGGTAAAAGAGCGCCCAAGGTTTCCTTACC‑AAACCCCCCCAAACGAAAAcccaccacgcgggaggcaaaa                                                                                                                                                                       >  2:1027641/1‑220 (MQ=25)

                                                                       gCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCCCCCGATGAACAACTTCTTGCGGTGGAACCGCCACCCACTTTCTGCTTTACAAACCATTGGGTAAAAACCGCCCAAAAGGTGCTGTTCCAAACCCCCCCCAACAAAAAACccccaccggggagggaaaa                                                                                                                                                                        >  2:881638/1‑224 (MQ=25)

                                                                       gCACTGTCTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGGAAAAACGCCCAACGTGAGGACaaaa                                                                                                                                                                       >  1:33509/1‑239 (MQ=255)

                                                                            gTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTTTGCTTAACAACCCATTGAGTAAAATAGGGCCAAAGGGTGCTGTACC‑AATCCCTCCCAACGTAAAACCCCCCAGCGTGAGGGCAAAAtacc                                                                                                                                                                   >  1:882197/1‑235 (MQ=39)

                                                                             ttttGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGGGGAACGTCAAACCACTTTCTCCTTTACAAAGCCTTGAGTAAAAACGCGCCCAAGGTTTCTGTACC‑AATCCCTCCCAACGTAAAAACcccccccgggaagacaaaaaaacac                                                                                                                                                                 >  2:902118/1‑218 (MQ=37)

                                                                             tttggtggccgggATAGGGTGTTT‑CCGCCGCCT‑CGGGAATTAACAATCAGTTCGTTCT‑CAACAATCTGCCCCGCCGGTACCCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCAccc                                                                                                                                                                 <  1:497318/228‑1 (MQ=255)

                                                                                  tAGGCCGGATAAGGCTTTT‑AC‑CCGCCTCCCGGAATAAACAATGCGTCCGTTGC‑AAAAAATCGGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGcc                                                                                                                                                            <  1:626252/239‑1 (MQ=255)

                                                                                  caGGCC‑GATAAGGGGTTTCACGCGGCAC‑CCGGAATGAAAAATGGTTAGGTTCT‑CAACAATCTGTACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGcc                                                                                                                                                            <  1:432473/238‑1 (MQ=255)

                                                                                       cGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATTAACAACTTCTGGCGGTGGAACGTCAACCCACTTTCTGCTTTACAAAGCATTGAGTAAAATCGCGCCAAAGGTTGCTGGACC‑ATACCCTCCCCACCGAAAACCGCCCAGCGTGAGAGCAAAAACACCCGCGCCCA‑Gc                                                                                                                                                       >  1:877683/1‑238 (MQ=255)

                                                                                       cGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGGACGTCAACCAACTTTCTGCTTAACAACGCCTTGAGTAAAATCGCGCCAAAGGTTTCTGTACC‑AATCCCTCCCAACGTAAAACAGCCCACCGTGGGAGGAAAATAACCCGCGCCCA‑Gcc                                                                                                                                                      >  1:398028/1‑238 (MQ=255)

                                                                                         gATAAGGGTTTT‑CCGCCGCTT‑CGGGCATAAACAGTGCGTAGGTTGC‑CAACAACCTGCCCCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCAc                                                                                                                                                     <  1:607224/238‑1 (MQ=255)

                                                                                         gATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATGCCGCCAAAGGTTGCTTGACC‑AATCCCTCCCAACGTAAAACCGCCCAGCCTGGGAGCCAAATAACCCGCGCCCacccct                                                                                                                                                     >  1:195503/1‑234 (MQ=255)

                                                                                                                            aatgatgaaggtgcaaaaactCTGCACCGCAGGGAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGtt                                                                                                                    <  1:1066719/219‑1 (MQ=255)

                                                                                                                             tGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGCAAAATCGCGCCCAAGGGTGCTGTACC‑AATCCCCTCCAACGTAAAACGACCCCCCGGCAGAGCAAAAACACCCCCGCCCACCCACCAAGG‑TACTCCCACCACCAACAAATTACCAttgt                                                                                                                  >  2:565019/1‑237 (MQ=255)

                                                                                                                                   cGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACCTAAAAACGCCCAGCGTGAGAGCAAAATCACCCCCGCCCA‑ACACTAACGTTACTTCGACCATAATCAAATTACCGATCTTGCTTa                                                                                                            >  1:830660/1‑239 (MQ=255)

                                                                                                                                       ttccaaacaTCTCCCCCGCCGGTAACCCCGGGGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCGGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAAt                                                                                                         <  2:398028/231‑1 (MQ=255)

                                                                                                                                            acattcgcccccgccgGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGaca                                                                                                    <  1:1052171/229‑1 (MQ=255)

                                                                                                                                             aattctcacccgccgGTACCCCCGGCGGTTTTCTGTTTATGGCTCCTAATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAAGCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGacac                                                                                                   <  1:1027641/230‑1 (MQ=255)

                                                                                                                                                   gctccgctgGTAGCCCCTGGTGTTTTCGTTTTATGGCTCCTGACAAACAACTTCTGGCTGGTGCACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGtttt                                                                                             <  1:434562/236‑1 (MQ=255)

                                                                                                                                                       ccccgGTAAACCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTAc                                                                                         <  1:842600/237‑1 (MQ=255)

                                                                                                                                                          cgGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCa                                                                                      <  1:428944/239‑1 (MQ=255)

                                                                                                                                                           ggTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCat                                                                                     <  1:1004357/239‑1 (MQ=255)

                                                                                                                                                                                   aTGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGGGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAACCGACACGGTTTTTTACCCATACCCttccccccccaccgccaaac                                                             >  2:892141/1‑221 (MQ=255)

                                                                                                                                                                                       ctctgattgaAAACCTTCGGGCGGT‑GAAAGCAAACCAATTTTCTGCTTAACAACGGATTTAGTAAAATCCCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGccc                                                        <  2:845642/232‑1 (MQ=255)

                                                                                                                                                                                       cTCCTGATTAACAATTTCTGGGGGGGGAACCGTAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTTTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGc                                                          <  1:415757/238‑1 (MQ=255)

                                                                                                                                                                                          cTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCCTTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACCGCGACCATAATCAAATTACCGTTCTCGCTTAAATCGGCACGGTTTTGTACCCAAATCCTTGAGCACGCGACGGCaaaaaccccaa                                                      >  2:653201/1‑230 (MQ=255)

                                                                                                                                                                                           tgatgaACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGaa                                                     <  1:881638/239‑1 (MQ=255)

                                                                                                                                                                                             atgaACAACTTCTTGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAaca                                                   <  1:33437/239‑1 (MQ=255)

                                                                                                                                                                                                aacaacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGCCAATCAGCCCGcacacac                                                >  1:845642/1‑238 (MQ=255)

                                                                                                                                                                                                aacaacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGACCA‑GCACTAACGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGCCACGGTTTTGTACCCATATCATTCCGCCTGCGACCGCAATCACCCCGAAAacaa                                                >  2:348828/1‑239 (MQ=255)

                                                                                                                                                                                                 acaacTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTACGGTTACTGCGACCATACTCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCCTAGCCTTGCGCCTGCGACGGCAATCAGCCCGAACacaac                                               >  1:45761/1‑239 (MQ=255)

                                                                                                                                                                                                    acTTCTGGGGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGGTCC‑AATCCCTCCCAACGTAAAGCCGTCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAAt                                            <  1:565019/239‑1 (MQ=255)

                                                                                                                                                                                                     cTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCCCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGTCCCTGCGCCGGCAATCAACCCGAACACAACAAta                                           >  2:821896/1‑238 (MQ=255)

                                                                                                                                                                                                            gggTGGAAACTCAACCAACTTTCTGCTTAACAACCCATTGAGTAAAATCGCCCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGca                                     <  2:830660/237‑1 (MQ=255)

                                                                                                                                                                                                               tGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGATCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACAAAACAATTGATGcccac                                  >  1:628500/1‑234 (MQ=255)

                                                                                                                                                                                                                    cGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAAt                             <  1:902118/238‑1 (MQ=255)

                                                                                                                                                                                                                            aacttttctcTTAACAACCCATTGAGTAAAATCGCGCCAAAGGTTTCTGTACC‑AATCCCCCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAGCCGCGGc                    <  2:33509/230‑1 (MQ=255)

                                                                                                                                                                                                                              cTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACcgcg                      <  1:803705/235‑1 (MQ=255)

                                                                                                                                                                                                                                 tccgattAACAACCCATTGAGTAAAATCGCGCCAAAGTTTGCTGGACC‑AATCCATCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGaa               <  1:326420/234‑1 (MQ=255)

                                                                                                                                                                                                                                 tCTGTTTAACACCGCATTGAATAAAACCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGGGAGAACAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTTCCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGaa               <  2:195503/239‑1 (MQ=255)

                                                                                                                                                                                                                                 tCTGCTTAACAACGCATGGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGaa               <  1:205226/239‑1 (MQ=255)

                                                                                                                                                                                                                                        aacaacGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCCCCACAAACCCCGCCCGAAATGGtcgg        >  2:803705/1‑236 (MQ=255)

                                                                                                                                                                                                                                        aacaacGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAACTGATGCACCACCAACAACCGCGGCCGGAAAGGAATg        >  2:669268/1‑239 (MQ=255)

                                                                                                                                                                                                                                              gCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAg  <  1:188141/239‑1 (MQ=255)

TGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATA‑TATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTT‑ACGCCGCAT‑CCGGCATGAACAATGCGTACGTTGT‑CAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATCGCGCCAAAGGTTGCTGTACC‑AATCCCTCCCAACGTAAAACCGCCCAGCGTGAGAGCAAAATCACCCGCGCCCA‑GCACTAAGGTTACTGCGACCATAATCAAATTACCGTTCTGGCTTAAATCGACACGGTTTTGTACCCATATCCTTGCGCCTGCGACGGCAATCAGCCCGAACACAACAATTGATGCACCACCAATAACCGCGGCCGGAATGGTATGAATCAG  >  NC_000913/1068338‑1068810

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