Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,499,609	A→C	22.9%	A54A (GCA→GCC)	yhfL →	small lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,499,609	0	A	C	22.9%	107.0 / 5.0	61	A54A (GCA→GCC)	yhfL	small lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (10/36);  new base C (11/4);  total (21/40)
Fisher's exact test for biased strand distribution p-value = 5.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAACATCATATGTTGTTATATTCATCATGCATTGTCA  >  NC_000913/3499393‑3499845                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         ttATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAAt                                                                                                                                                                                                                                  <  1:142789/238‑1 (MQ=255)                            attaatATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAAt                                                                                                                                                                                                                                  <  1:621338/211‑1 (MQ=255)                                      aataCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTTCGGGCCCACTTCAC‑AGCAAGTCTCTCAC‑AG‑AAAACTGGT‑AAACC‑CCGGTTTGTTTGCCCA‑CCC‑CCCC‑CT‑ggt                                                                                                                                                                                      >  1:748088/1‑238 (MQ=255)                                        taCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑GTGCGGTCCCACTGCCC‑ACGAAATCTCTGAT‑AG‑AAAAC‑GCCATAACA‑CGGCTTATATTCCCCA‑CCG‑cccctggcatt                                                                                                                                                                                      >  1:276888/1‑231 (MQ=255)                                             aGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTTCGGTCCCTCTCCCC‑ACACAAGCTCTTAC‑AG‑AAAAC‑GGC‑TAAAA‑ACCGGTTAATTTCCCC‑ACC‑CCCC‑CT‑TGG‑TTTTctcc                                                                                                                                                                             >  1:263292/1‑236 (MQ=255)                                             aGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGa                                                                                                                                                                             <  1:598280/239‑1 (MQ=255)                                             aGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGAGGTCCTACTGCAC‑AGTAAGTCTCTCTC‑AG‑AAAAC‑GGC‑TCAAA‑CCGCGTTAATTTCCCC‑CCC‑CCCC‑CT‑GGG‑TTTTctca                                                                                                                                                                             >  1:292916/1‑236 (MQ=255)                                              gggTATTTAAATAAAAAAATTTAATAAAATTGCAATAGTAGGTCCAGTAATGGCTTATTTAACTCGACTGAATGGTCCTTGTGTAATCCGGTATCAGAGCTGGACTTACGACTGCCTGCTTCTCCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCAT‑CT‑GGA‑TTTTCCGAg                                                                                                                                                                            <  2:1129087/239‑1 (MQ=255)                                               ggAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTTGGGTCCTACTGCCCAAGAATATTTCTCAC‑AA‑AAATA‑GG‑‑TGAAA‑ACGGGTTAATTTCCCG‑GCG‑CCCC‑CTGGGA‑TTTTCCggg                                                                                                                                                                            >  2:620975/1‑237 (MQ=255)                                               ggAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑GTGCGGTCCCACTGACC‑AACAAGTTCCTGAC‑AG‑AAAAC‑CGG‑TAACA‑CCCGGTTCTTTGCCCC‑ACG‑CCCC‑CT‑GGACTTTTCCggg                                                                                                                                                                            >  1:790905/1‑237 (MQ=255)                                                   tGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTTCGGTGCCACTTCCC‑CGCAAGTCTCTATC‑CA‑CAATC‑ACC‑TCAAA‑CACGTTTAATTGTCCGCAGC‑CCCA‑CT‑GGA‑TTTTCCCAggcgg                                                                                                                                                                        >  1:661939/1‑236 (MQ=255)                                                   tGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTCCTGCCC‑AGTAAGGCTCTGTC‑AG‑AAAAC‑GGC‑TAACACCCGGTTTATTTTCCC‑‑ACGCCCCC‑CT‑GGT‑TTTTcccgggggg                                                                                                                                                                        >  1:128606/1‑232 (MQ=255)                                                   tGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGGCCTACTACAC‑AGCAAGTGTCTCAC‑CA‑CAAAC‑GGC‑TAAAA‑CCCCGGTAATTTTCCC‑ACG‑CCCT‑CT‑GGG‑tttttccaaggggg                                                                                                                                                                       >  1:181996/1‑229 (MQ=255)                                                             cAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTTGCGG‑CCCACTTACC‑AGTAAGTCTCTTAC‑AGAAAAAC‑GGC‑AACAA‑CCGGTTTTATT‑CCCC‑CCG‑CCCC‑CTGGGT‑TTTTCCGAGGGGGTATttttttt                                                                                                                                                              >  1:1264361/1‑238 (MQ=255)                                                                aTTTTTTTAATTTTCACTGGTAGGTGCAGTACTGGCTTCGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAT‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATg                                                                                                                                                          <  2:530601/239‑1 (MQ=255)                                                                     ttaAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑ATGCGGTCCTACTCCCC‑AGTAAATCTTTGAC‑CG‑GAAAA‑GGC‑TAACA‑CCCGTTTAATTTCCCC‑ACC‑GCCC‑CT‑GGG‑TTTTTCGCGGGCGGAGTTTTTTCTATTTTctt                                                                                                                                                     >  1:322958/1‑239 (MQ=255)                                                                                             tACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCCC‑CGATAATTTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGGTTAATTTCCCG‑CGG‑CCCC‑CT‑GGG‑TTTTTCCCGGGGGGGTTTTTTTTTTTGTCCTTC‑TTTCCCCCCCCCCCAAAAAAAtt                                                                                                                            >  1:106786/1‑239 (MQ=255)                                                                                                  gCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTTCGGTCCCACCTCAC‑ACTAAATCTCTTAC‑CA‑GAATC‑AGC‑CAACA‑CCCGCTTAATTTTCCG‑GCC‑CCCC‑CC‑TGA‑TTTTTCCACGGCGTCATTTTTTTTTTTTTCTTC‑TTTCCCCCCCCCCTcaaacaaatttt                                                                                                                         >  2:144019/1‑229 (MQ=255)                                                                                                  gCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTTCGGGCCTACTGCAC‑CGTAAGTCTCTGTC‑AG‑AAAAC‑GGC‑CTACA‑CCCGTTTAATTGCCCC‑GCG‑CCCC‑CT‑GGC‑TTTTTCCAGGGCGGATTTGTTTTTTTGTCATTC‑TTTTCCCCCCCCCTCAAAAAAttttt                                                                                                                         >  1:414282/1‑236 (MQ=255)                                                                                                       gTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTTGTGTCCTTCTCCCC‑AATAAGTCTCTGCC‑CG‑AAAAT‑GGG‑TAACA‑CCGGGGTAATTGCCCGCACG‑‑CCT‑CT‑GGAGTTTCCCGAGGGGGCTTTTGTTTTTTTG‑CCTTCTTTCCCCCCCCCCCGAAAAATTTTTTCCTTTTa                                                                                                                   >  1:117935/1‑239 (MQ=255)                                                                                                       gTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAATCTCTTACAAA‑AAATG‑GCC‑AAACC‑CCGTTTTATTTCCCCG‑CGC‑CCCC‑CT‑GGT‑TTTTCC‑AGGGGGGATTTTTTTTTATTTCCTTC‑TTTGCCCCCCCCCTCAAAAAATTTTTCCTTTaa                                                                                                                  >  1:982356/1‑239 (MQ=255)                                                                                                        ttAACGTCATGTTCTGGTCATATTGAAAACCGGTATAAGAACTGCTCTTACGGCTTCCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTaaa                                                                                                                 <  2:790905/239‑1 (MQ=255)                                                                                                         tAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTaaa                                                                                                                 <  1:923716/238‑1 (MQ=255)                                                                                                                 tgcaCTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTCCCGCAC‑AGCAAATCTTTTCC‑AG‑AAAAC‑AGC‑TTACA‑CCCACTTGATTGCCCG‑CCC‑CCCT‑CT‑CGG‑TTTTCCGAGGGGGCGGTTTTTTTTTTTTTCTTCTTTCCCCCCCCCCGCAAAAAAATTTTCTTTTTAATTTAAAttt                                                                                                         >  2:266657/4‑239 (MQ=255)                                                                                                                  tttcTGGTCATATTGAAAACCGTTATAAGAACTTCTCTTACGGCTTCCTGCTTCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGa                                                                                                       <  2:414282/236‑1 (MQ=255)                                                                                                                   tACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGa                                                                                                                                                                             <  1:213764/169‑1 (MQ=255)                                                                                                                   tACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGAt                                                                                                      <  1:266657/239‑1 (MQ=255)                                                                                                                    aCTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGAt                                                                                                      <  1:54940/238‑1 (MQ=255)                                                                                                                          cATTTTGAAAACCGTGATAAGAACTGCTCTTACGGCTACCTTCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTAc                                                                                                <  2:263292/238‑1 (MQ=255)                                                                                                                           atatTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGTGGTCCCACTGCCC‑CTGAAGTTACTGAC‑AG‑AAAAC‑GGC‑TCACA‑CCCGTTTTATTTCTCT‑CCG‑CCCTCCT‑GGT‑TTTTCCAAGGG‑GGATTTGTTTCTTTGAGCTTT‑TTTGCCCCGCCCGCGCAAAAATATTTACTTTTAATTCAATTTTTTAACTACa                                                                                               >  1:251274/1‑238 (MQ=255)                                                                                                                             atTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGGCCTAATGCAC‑AGTAAGTCTCTGAC‑AG‑AAATA‑GGC‑TAACA‑CCCGTTTGATTGCCCG‑ACG‑CCCT‑CT‑GGG‑TTTTCCGATGGGGGATTTGTTTTTCTGTCCTTC‑TTGTCCCCGCCCGTACAAATATTTTTACTTTAAATTTCATTTGTTAACAACATc                                                                                             >  1:1138901/1‑238 (MQ=255)                                                                                                                                   aaCCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGc                                                                                      <  1:1225070/239‑1 (MQ=255)                                                                                                                                     ccGTTATAAGAACTGCTCTTACGACTACCTGCTGCCCCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCca                                                                                    <  1:144019/239‑1 (MQ=255)                                                                                                                                          aTAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑Ac                                                                                                                                                                                                   <  1:313549/128‑1 (MQ=255)                                                                                                                                          aTAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGAt                                                                                                      <  1:1051739/216‑1 (MQ=255)                                                                                                                                          aTAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTAc                                                                                                <  1:244579/222‑1 (MQ=255)                                                                                                                                          aTAAGAACTGCGCTTACGACTTCCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCca                                                                                    <  1:879183/234‑1 (MQ=255)                                                                                                                                                aCTGCTCCTGCTACTACCTGCTGCACCCGGCAAGTTCTATTTCTAAAATCATTGGCGG‑TTGCCGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAACC‑GGC‑TAACT‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATaa                                                                                                    <  1:99542/212‑1 (MQ=255)                                                                                                                                                     tctTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCcccgc                                                                                                                                         <  1:328563/170‑1 (MQ=255)                                                                                                                                                               gccTGCCGCCACCGGCAATTTCTATGTCAAAAATCATTGGCGG‑TTGCGGGTCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTtgt                                                          <  1:675723/238‑1 (MQ=255)                                                                                                                                                                      gcACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGTCACAGCGCGATCCTGTGCGAAATTTTGTGATCTtc                                                   <  1:2026/239‑1 (MQ=255)                                                                                                                                                                       cACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTTCGGGCCTACTGCAC‑AGTAAGTCTCTTAC‑AG‑AAAAC‑GGC‑TAACA‑CCCGGTTTATTTCCCC‑ACC‑CCCC‑CT‑GGA‑GTTTTCCAGGGGGGATTTGTTTTTTTTTCCCTC‑TTTTCCCCCCCCCCAAAAAAAATTTTTCTTTTAATTTAAATTTATAACTACATTCATTTTCCCCACCGCCATCCCTTTCGGAAATTTGTGATCTtc                                                   >  1:264811/1‑238 (MQ=255)                                                                                                                                                                        ccccGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCCGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTtcct                                                 <  1:620975/238‑1 (MQ=255)                                                                                                                                                                        aCCCGGCAATTTCTATTTCCAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTtcct                                                 <  1:790250/239‑1 (MQ=255)                                                                                                                                                                           cGGCAATTTCTATTTTTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCca                                              <  1:573017/239‑1 (MQ=255)                                                                                                                                                                           cGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCca                                              <  1:234096/239‑1 (MQ=255)                                                                                                                                                                             gctatttttttttctAAAATCATTGGCGG‑TTTCGGTCCTACTGCAC‑AGTAAGTCTTTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCcaca                                            <  2:982356/236‑1 (MQ=255)                                                                                                                                                                             gCAATTTCTTTTTCTAAAATGATTGGCGG‑TTGCGGGCCTACTGCAC‑ATTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCTTATTTGTTTATATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCcaca                                            <  2:276888/239‑1 (MQ=255)                                                                                                                                                                             gCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑Ac                                                                                                                                                                                                   <  1:309406/93‑1 (MQ=255)                                                                                                                                                                             gCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCcaca                                            <  1:753481/239‑1 (MQ=255)                                                                                                                                                                             gCAATTGCTAGTTCTAAAATCATGGGCGG‑TTGCGGGCCTACCGCAC‑AGTAAGTCTCGGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAAGTTGATAACTACATCCATTAGCCACATCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCcaca                                            <  1:584920/239‑1 (MQ=255)                                                                                                                                                                              cAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGGCTCTGAC‑AG‑AAATC‑GGC‑TAACT‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCcaca                                            <  1:170024/238‑1 (MQ=255)                                                                                                                                                                                  ttctatttctaAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACa                                       <  1:638748/239‑1 (MQ=255)                                                                                                                                                                                   tttttttctaAAATCTTTGGCGT‑TTGCGGTCCTACTGCCC‑AGTAAGTTTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACAt                                      <  2:322958/235‑1 (MQ=255)                                                                                                                                                                                        ttcttAAATTCTTGGCGG‑TTTGCGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGTTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAAcat                                 <  2:661939/239‑1 (MQ=255)                                                                                                                                                                                         tctaAAATCCTGTGCGG‑TTGCGGGCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAAcatc                                <  2:292916/239‑1 (MQ=255)                                                                                                                                                                                            aaaaTCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACC‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAACATCata                             <  1:593032/239‑1 (MQ=255)                                                                                                                                                                                              aaTCATTGGCGG‑TTTGGGTCCTTCTGCCC‑AGTAAGTCTTTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAACATCATAtg                           <  2:117935/239‑1 (MQ=255)                                                                                                                                                                                                     ggCGG‑TTGCGGTCCTACTCCCC‑AGTAAGTCTCGGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAACATCATATGTTGTtat                    <  2:1323562/239‑1 (MQ=255)                                                                                                                                                                                                                 tCCTACTGGCC‑AGTAAGTCCTTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAACATCATATGTTGTTATATTCATCATg          <  2:106786/238‑1 (MQ=255)                                                                                                                                                                                                                        gCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAACATCATATGTTGTTTTATTTATCCTGCAtttgca  >  1:659157/1‑235 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         TTATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGG‑TTGCGGTCCTACTGCAC‑AGTAAGTCTCTGAC‑AG‑AAATC‑GGC‑TAACA‑CCGGTTTAATTGCCCG‑ACG‑CCCT‑CT‑GGA‑TTTTCCGAGGGCGTATTTGTTTCTATGTCCTTC‑TTGCCCCCGCCCGTAAAAATATTTTCACTTTAAATTCAATTTGATAACTACATCCATTAGCCACAGCGCGATCCTGTGCGAAATTTTGTGATCTTCCTCCACATTACATAACATCATATGTTGTTATATTCATCATGCATTGTCA  >  NC_000913/3499393‑3499845

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