Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,665,702	0	C	A	48.6%	33.2 / 53.9	70	A37S (GCG→TCG)	gadX	acid resistance regulon transcriptional activator; autoactivator
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base A (0/34);  total (17/53)
Fisher's exact test for biased strand distribution p-value = 1.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAA  >  NC_000913/3665554‑3665850                                                                                                                                                        |                                                                                                                                                     tgatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCaa                                                                                                                                                    <  1:69741/151‑1 (MQ=255)  gatgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCaaa                                                                                                                                                   >  2:780421/1‑151 (MQ=255)    tgaACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATTTACTCGAATTTGG‑CTTGCATCCGCaaaaa                                                                                                                                                 >  2:639086/1‑151 (MQ=255)      aaCGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAAcc                                                                                                                                               <  2:145328/151‑1 (MQ=255)      aaCGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAAcc                                                                                                                                               <  1:269658/151‑1 (MQ=255)        cGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAg                                                                                                                                             >  1:366717/1‑151 (MQ=255)            tAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTc                                                                                                                                          <  2:650802/150‑1 (MQ=255)                 aagaaACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGc                                                                                                                                     <  2:729426/150‑1 (MQ=255)                       cTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTaaa                                                                                                                              >  1:23660/1‑151 (MQ=255)                         tttttATGTGCATGAAGATTTACTGCCTCCTCCTTGAGCATGGGGAGAAATAATAAAAGCGTGTCCCTT‑GAC‑ACAAAAAAAAAATTTCCAACCCACTTATCTACTCAATTTTGG‑CTTCCCTCCGAAAAAACCAGGTCCCCGCCATTAAAAt                                                                                                                            <  1:318512/151‑1 (MQ=255)                         tttttATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAAAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATAAAAAATTTTCAACACATTTATCTACTCGAATTTGG‑CTTCCATCCGAAAAAACCAGGTCACCGCCATTAAAAt                                                                                                                            <  2:175503/151‑1 (MQ=255)                         tttttATGGGCAAGAAGATTTAAGCCCTCCTCCGTGAGCATGGGGAAAAATAATAAAAGCGT‑TCCCTTAAAC‑ACAAATAAAAAATTTTCAACACACTTATCTCCTCAAATTTGT‑CTTGCACCCGCAAAACCCAGGTCACCGCCATTAAAAt                                                                                                                            <  1:782297/151‑1 (MQ=255)                                   aaagAAGATTAAATCCCTCCTCCTTGAGCATCGGGAAAAAAAATAAAAGCGTGTCCCTT‑AAC‑ACAAATAAAAAATTTCAAACACACTTATCTACTAAAATTTGG‑TTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                                                                                                  <  1:89195/148‑1 (MQ=255)                                       aaaaTTTAATCCCCCCCCCTTGAGCATCGGGAAAAAAAATGAAAGCGTGTCCCTT‑AAC‑CCGAATAAAAAATTTTCAACCCTTTTATCTCCTCGATTTTGG‑CTTGCACCCCCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                                                                                               <  2:177781/147‑1 (MQ=255)                                         gttttAAGGCCTCCTCCTTGAGCAAGGGGAAAAATAATGAAAGCGTTTCCTTT‑GAC‑ACGAAAAAAAAATTTGCAACACCCTTATCTGCGAAATTTGGG‑CTTCCAGCCGAAAAAACCGGGTCGCCGCCATTAAAATAGCGATATTCACCAt                                                                                                             <  2:169225/148‑1 (MQ=255)                                         gATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                                                                            >  2:109671/1‑151 (MQ=255)                                         gATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                                                                            >  1:76387/1‑151 (MQ=255)                                         gATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCAt                                                                                                             >  1:169191/1‑150 (MQ=255)                                          ttttAACCCCTCCTCCTTGAACATCGGGAAAAATAATGAAAGCGTCCCCCTT‑GAC‑ACAAAAAAAAATTTTACACCCCATTTATCTACTCAATTTTGG‑TTTGCACCCACAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCAt                                                                                                             <  1:485236/148‑1 (MQ=255)                                             tAATGCCTCCTCCTTGAGCATCGGAAAAAATAATGAAGGCGTTCCCCTT‑GCC‑ACGAAAAAAAATTTTTCAACACACTTATCTACAAGATTTTGG‑CTTCCACCCAAAAAACCCAGTTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                        <  1:22708/151‑1 (MQ=255)                                             tAATCCCTCCTCCTTGAGCATGGGGAAAAAAAATAAAGGCGTTTCCTTT‑GAC‑AAGAATAAAAAATTTTAAACACACTTTTCTACACAATTTGGG‑TTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                        <  1:707521/151‑1 (MQ=255)                                             tAATCCCCCCTCCTGGAGCATGGGGAAAAAAAAGAAAAGCGTGTCCTTT‑GAC‑ACAAATAAAAAATTTTCAACACCTTTATCTACTCAATTTTGG‑TTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                        <  2:482762/151‑1 (MQ=255)                                             tAATCCCCCCCCCTTGAGCATCGGGAAAAAAAATGAAAGCGTGCCCCTT‑GAC‑ACAAATAAAAAATTTTCAACACCCTTATCTACCCAAATTTGG‑CTTCCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACAATTAAcc                                                                                                        <  2:46835/151‑1 (MQ=255)                                                       tcctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                                                                              >  2:553593/1‑151 (MQ=255)                                                          ttGAGCATGGGGAAAAATAATGAAAGCGTGTCCCTT‑AAC‑ACAAATAAAAAATTTCAAACACCCTTATCTACCCGAATTTGG‑CTGGCACCCGAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatat                                                                                           <  2:555758/151‑1 (MQ=255)                                                            gAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatatat                                                                                         >  1:142774/1‑151 (MQ=255)                                                               cATCGGGAAAAAAAAAAAAGGGGTGTCCTTG‑GAC‑ACAAAAAAAAATTTTTCAACACTCTTATCTCCCCGAATTTGG‑CTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                        <  2:596701/149‑1 (MQ=255)                                                                 tCGGGAAAAATAATGAAGGCGTGTCCCTT‑AAC‑ACAAATACAAAATTTTCAACACATTTATCTCCTCGAATTTGG‑TTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATg                                                                                    <  2:646002/151‑1 (MQ=255)                                                                  cggaaaaaaTAATGAAGGCGTTTCCCTT‑GAC‑ACAAATAAAAAATTTTCAACACATTTATCTCCTCGAATTTGG‑CTTGCACCCGAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                                   <  2:50406/145‑1 (MQ=255)                                                                  cGGGAAAAATAATGAAAGCGTGTCCCTT‑GAC‑ACAAATAAAAAATTTCCAACACACTTATCTACAAGAATTTGG‑CTTCAATCCAAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                                   <  2:236816/151‑1 (MQ=255)                                                                   gggAAAAAAAAGAAACGCGTGTCCTTG‑GAC‑ACAAATAAAAAATTTTCAACACACTTATCTCCTCAAATTTGG‑CTCCCAACCGAAAAAACCGGGTAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                                   <  1:543089/150‑1 (MQ=255)                                                                   gagaaaaaTAATGAAGGCGTG‑CCCTT‑GACAACGAAAAAAAATTTTCCAACCCCCTTATCTCCACAATTTGGG‑TTTCCACCCACAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                                   <  1:264268/145‑1 (MQ=255)                                                                       gaAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTg                                                                               >  2:211543/1‑150 (MQ=255)                                                                             atGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCaa                                                                                                                                                    >  1:356392/1‑75 (MQ=255)                                                                             atGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCaa                                                                                                                                                    <  2:356575/75‑1 (MQ=255)                                                                                     ttcccccTT‑GAC‑CAAAAAAAAAAATTTACACCCCCTTTATCTACT‑AAATTTGGTCTTGCACCCAAAAAACCGGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                                <  1:19074/147‑1 (MQ=255)                                                                                      ttcccttt‑GAC‑AAAAAAAAAAAATTTCCACCCCCCTTATCTACTGAAATTTGG‑CTTGCACCCACAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                                <  2:109483/144‑1 (MQ=255)                                                                                      tgtCCTTT‑GAC‑ACAAATAAAAAATTTTAAACCCTCTTATCTCCTCAATTTTGG‑CTACCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                                <  2:212584/150‑1 (MQ=255)                                                                                      tgtCCTTG‑GAC‑ACGAATAAAAAATTTTCAACACATTTATCTACTAGATTTTGC‑CTTGCATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                                <  1:739741/150‑1 (MQ=255)                                                                                      tgtCCCTA‑GAC‑AAAAAAAAAAAATTTCCAACACACTTATCTCCTCGATTTTGC‑CTTGCACCCAAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                                <  1:395702/150‑1 (MQ=255)                                                                                      tgcccttt‑GAC‑AAAAATAAAAAATTTAAAACACACTTATCTCCTCAATTTTGG‑CTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTACGTCTTGCATACGCAATTAGAc                                                                <  2:633060/144‑1 (MQ=255)                                                                                         cccTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAAtt                                                            >  2:702432/1‑151 (MQ=255)                                                                                           ccgcgAC‑ACGAATAAAAAATTTCCATCAACTTTTTCAACAGTATTTTGG‑TTTCCACCCAAAAAACCGAGTTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTc                                                           <  1:652724/147‑1 (MQ=255)                                                                                            ctggAC‑ACAAAAAAAAATTTTTCAACCCTCTTATCTACTCAATTTGGG‑TTTGCACCCGAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTcc                                                          <  2:65513/148‑1 (MQ=255)                                                                                                 c‑acaAATAAAAATTTTCCACCCCACTTATCTACTCGAATTTGG‑CTTGCACCCGAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                                     <  2:678339/151‑1 (MQ=255)                                                                                                 c‑acGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                                     >  2:713836/1‑151 (MQ=255)                                                                                                  cccaaaTAAAAAATTTCAAACCCACTTATCTACCCGAATTTGG‑CTGGCATCCACAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                                     <  2:262768/147‑1 (MQ=255)                                                                                                  cacaaaaaaaaaaTTTAAAACACACTTATCTACTCGAATTGGG‑CTTGCACCCGAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                                     <  2:508304/142‑1 (MQ=255)                                                                                                    aaaaaaaaaaaTTTACACCACCTTAATCTATTGAATTTTGT‑CTTCCACCCCAAAAACCCAGTCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                                   <  1:503615/144‑1 (MQ=255)                                                                                                      aaTACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTg                                                 >  1:745639/1‑151 (MQ=255)                                                                                                        aaaaaaaTTTACAACACATTTATCTACTAAAATTTGG‑CTTCCACCCGCAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGAt                                               <  1:140238/148‑1 (MQ=255)                                                                                                                  tCAACACACTTATCTACTCGATTTGGG‑CTGGCATCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTg                                     <  2:556354/151‑1 (MQ=255)                                                                                                                    aaCACATTTTTCTACTCGATTTTGG‑TTTCCACCCGCAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGAc                                   <  2:463756/151‑1 (MQ=255)                                                                                                                       ccctttACCTACTCAAATTTGG‑TTTCCATCAAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTa                                <  1:684116/147‑1 (MQ=255)                                                                                                                          cTTATCTACTCGAATTTGG‑CGTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGCAGTGATTGCATAGTTGAATTAata                             >  1:791393/1‑151 (MQ=255)                                                                                                                           ttttctcCTCAAATTTGG‑CTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAatat                            <  2:200822/148‑1 (MQ=255)                                                                                                                           ttATCTACTAAATTTTGG‑CTTGCACCCGAAAAACCCAGCTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAatat                            <  1:377053/151‑1 (MQ=255)                                                                                                                                tccTCGATTTTGC‑CTTGCACCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACAt                       <  1:761899/149‑1 (MQ=255)                                                                                                                                taccagAATTGGG‑CTTGAACCCACAAAACCGAGGCCGCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACAt                       <  2:628860/146‑1 (MQ=255)                                                                                                                                      aaTTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACat                 >  1:471701/1‑151 (MQ=255)                                                                                                                                      aaTTGGG‑TTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACat                 <  1:785814/151‑1 (MQ=255)                                                                                                                                        tttGG‑CTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACatat               <  2:773126/151‑1 (MQ=255)                                                                                                                                          ttgctttGCATCCACAAAAACAAGGCATCCGCCATTAAATTAGCGATATCCACCATTAACCATAGAGAGAAGATATTTATGTCTTTCATACGCAATTAGACAATTCCCAGGTAGAGATTGCATAGTTGACTTAATATTACATAAGCATAtt              <  1:512281/146‑1 (MQ=255)                                                                                                                                          tggcttaccacccGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATAtt              <  1:3324/140‑1 (MQ=255)                                                                                                                                            g‑CTTGCATCCGAAAAAACAAGGTCGCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTAc            <  1:772530/150‑1 (MQ=255)                                                                                                                                            g‑CTTGCATCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTAc            <  2:125229/150‑1 (MQ=255)                                                                                                                                            g‑CTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTAc            <  1:713505/150‑1 (MQ=255)                                                                                                                                              tttGCACCCGAAAAAACAAGGACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTg          <  1:402582/150‑1 (MQ=255)                                                                                                                                                     ccGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTttatta   >  1:673882/1‑151 (MQ=255)                                                                                                                                                      cacAAAAACAAGTTAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTaa  <  1:114756/149‑1 (MQ=255)                                                                                                                                                      cGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTaa  <  2:455293/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                     TGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTT‑GAC‑ACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGG‑CTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAA  >  NC_000913/3665554‑3665850

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