Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,378,802	0	A	C	43.4%	141.3 / 16.8	77	V48G (GTA→GGA)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (41/2);  new base C (0/33);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.66e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378653‑4378948                                                                                                                                                             |                                                                                                                                                    ccccTCGAGCAATAGCGTCGCCCCCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑a                                                                                                                                                     <  1:453381/150‑1 (MQ=255)  cccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑acc                                                                                                                                                   >  2:849305/1‑151 (MQ=255)  cccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑acc                                                                                                                                                   >  2:756444/1‑151 (MQ=255)     tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑accg‑cc                                                                                                                                               >  1:843295/1‑151 (MQ=255)         gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑accg‑ccaccg                                                                                                                                           >  1:939955/1‑151 (MQ=255)         gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑accg‑ccaccg                                                                                                                                           >  2:951285/1‑151 (MQ=255)         gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑accg‑ccaccg                                                                                                                                           <  2:834428/151‑1 (MQ=255)            aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑accg‑ccaccgcc                                                                                                                                         >  1:1052901/1‑150 (MQ=255)             tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCAt                                                                                                                                       >  2:1132728/1‑151 (MQ=255)               gCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATAc                                                                                                                                     <  2:265477/151‑1 (MQ=255)                cGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATAcc                                                                                                                                    >  2:1128551/1‑151 (MQ=255)                cGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATAcc                                                                                                                                    >  2:1190893/1‑151 (MQ=255)                 gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATAccc                                                                                                                                   >  2:304590/1‑151 (MQ=255)                  tCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCg                                                                                                                                  >  1:89434/1‑151 (MQ=255)                   cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGCACCTAA‑CTCAAACAACGGTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCg                                                                                                                                  >  1:875986/1‑150 (MQ=255)                   cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCgg                                                                                                                                 >  2:1158185/1‑151 (MQ=255)                      caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGAt                                                                                                                              >  2:382757/1‑151 (MQ=255)                         caAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCtt                                                                                                                           >  2:193329/1‑151 (MQ=255)                             cccgcaataaatttttttGCGCACCAACCCTAA‑CAAAAACAACTTTTGTGGGGGGGGGGGTTTTTTTTGGGCGATGCCCGCATCCCCCAAGAAAAAG‑TAA‑AAGG‑TTTTCCC‑TT‑AAAAATTT‑CCCC‑CCACCGCCATACCCGGGATc                                                                                                                             <  1:475713/131‑1 (MQ=255)                                 ccAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGct                                                                                                                   >  1:1034766/1‑151 (MQ=255)                                   aGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 >  1:499645/1‑151 (MQ=255)                                    gTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 >  2:431153/1‑150 (MQ=255)                                    gTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 >  2:809154/1‑150 (MQ=255)                                      aatttttttGCTGACCG‑ACCTAACCTAAACCAACTTTTGTTGTGTGGCGGGTTGTTTTTTGGCGATGCCCGCATACCCCAAGGTAAAG‑AAATAAGT‑TTCCCCC‑TA‑AAAAT‑TC‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 <  2:528717/143‑1 (MQ=255)                                      aaaTTTTTTGC‑GGCCGAACCCAACCTAAAAAAAGTTTTGTTGTGGGAGGGGCGGTTTTTGGGCGAGGCCCGCAT‑GGCCCAGGAAAAA‑TAAAAAGGTTTTCCCC‑TA‑AAAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                              <  2:823915/151‑1 (MQ=255)                                           ttttGCTGACCGACCCTAT‑CTAAAACAACGTTTTCTGTGGGAGGGGTTTTTTTTTGGCGATGCCCGAATACCCCCTGGTAAAG‑AAATAAGT‑TTTCCCC‑TG‑AAAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                          <  2:215653/150‑1 (MQ=255)                                            tttGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                        >  1:1008696/1‑151 (MQ=255)                                                     cGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                               >  1:959031/1‑151 (MQ=255)                                                      gAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                              >  2:581945/1‑151 (MQ=255)                                                       acccaaa‑CACAAACAACGTTTGGTGTGTGGCGGGCTGTTTTTGGGCGATGTCCGCATACCCCAAGGAAAAG‑TAATAAGT‑TTTCCCC‑TA‑ATAAA‑TC‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                     <  2:398542/138‑1 (MQ=255)                                                            aa‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                           >  1:1019534/1‑148 (MQ=255)                                                             a‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                       >  2:1015572/1‑151 (MQ=255)                                                                  aaaCATCTTTTGCTTTGTGAGGGGCTTTTTTTTGGCGGTGCCCCCAGGCCCCCAGAAAAAG‑AAAAAAGT‑TTTCCCC‑TA‑AAAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                       <  1:416854/147‑1 (MQ=255)                                                                  aaaCACCTTTTGCTGTGTGAGGGGCTGTTTTTTGGCGATGCCCCCATACCCCCAGAAAAAG‑AAATAGGT‑TTTCCCC‑TA‑AAAAA‑TC‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                          <  1:849291/144‑1 (MQ=255)                                                                  aaaCAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    >  1:1162052/1‑150 (MQ=255)                                                                  aaaCAACGTTTGCGTTGTGCCGGGTTTTTTTTTGGCGATTCCCCCATACCCCCAGGTAAAG‑TAATAAGT‑TTTCCCT‑TG‑AAAAA‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    <  2:311149/150‑1 (MQ=255)                                                                   aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                  >  1:977305/1‑151 (MQ=255)                                                                      aacTTTTGCTGTGTGAGGGGCTGTTTTTTGGCGATGCCCGCATACCCCAAGGAAAAT‑TAAAAAGT‑TTTCCCT‑TA‑ATAAA‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                               <  2:260274/151‑1 (MQ=255)                                                                       acGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                              >  1:1170001/1‑151 (MQ=255)                                                                          tttGCTGTGTGAGGGGCTTTTTTTGGGCGATTTCCGCATACCCCCAGGAAAAA‑TAATAGGT‑TTTCCCC‑TG‑AAAAT‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                             <  2:804310/149‑1 (MQ=255)                                                                          tttGCTGTGGGACGGGCTGTTTTTTGGCGATGCCCGCATACCCCCGGGAAAAT‑AAATAGGG‑TTTCCCC‑TA‑AAAAT‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                             <  2:606128/149‑1 (MQ=255)                                                                            tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCACAGCCCCAGGTAAAG‑GAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGAAATTGTGCGATGCACAATATCGTTg                                                                             >  1:681813/1‑147 (MQ=255)                                                                              cTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑AAATAAGG‑TTTACCC‑TG‑ATACA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                       >  2:673427/1‑151 (MQ=255)                                                                                  gtgACGGGTTTTTTTTGGGCGATGTCCGCATCCCCCCAGAAAAAG‑TAATAATT‑TTTCCCC‑TA‑ATAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                   <  1:595744/151‑1 (MQ=255)                                                                                  gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGACATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                     >  1:510720/1‑149 (MQ=255)                                                                                       gggTTGTTTTTTGGCGATGCCCGCAAACCCCCGGGAAAAT‑TAAGAAGT‑TTTCCCG‑TA‑AAAAT‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                  <  2:1034807/147‑1 (MQ=255)                                                                                       gggCTGTTTTTTGGCGATCCCCGCAGACCCCGAGAAAAAG‑AAAT‑AGG‑TTTTCCC‑TT‑AAAAA‑TTCCCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                 <  1:214627/148‑1 (MQ=255)                                                                                        gtttgtttttgGGCGATGCCCGCAGAGCCCAAGGTAAAG‑TAAAAGGT‑TTTCCCC‑TA‑AAAAT‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                             <  1:1147170/145‑1 (MQ=255)                                                                                         gctgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                            >  2:842772/1‑151 (MQ=255)                                                                                           tcttttttGGCATTGCCCACATACCCCCAAAAAAAT‑AAAGAATT‑TTTACCC‑TA‑ATAAT‑TC‑GCCG‑CCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                             <  1:839982/145‑1 (MQ=255)                                                                                             ttttttGGCGATCCCCGAATACCCCGAGAAAAAG‑AAAAAAGT‑TTCCCTC‑TA‑ATAAT‑TC‑CCCG‑CCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaaa                                                         <  1:524116/149‑1 (MQ=255)                                                                                              tttttGGCGCTGCCCGCAGTCCCCCAGCAAAAG‑TAATAGGT‑TTTCCCC‑TG‑ATAAA‑TC‑CCCG‑CCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                        <  2:105963/150‑1 (MQ=255)                                                                                               tttgGGCGATCCCCGCATACCCCCAGGAAAAG‑TAATAGGT‑TTTCCCC‑TA‑AAAAT‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                       <  2:1057376/150‑1 (MQ=255)                                                                                                  tGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                    >  2:397235/1‑150 (MQ=255)                                                                                                   gtcGATGCCCACATCTCCCAAGGTAAAA‑TAA‑AAGG‑GTTCCCC‑TGCAAAAA‑TT‑CCCG‑CCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                    <  1:49675/147‑1 (MQ=255)                                                                                                   ggCGTTGTCCGCATCCCCCCAGGAAAAG‑TAATAAGT‑TTTCCCC‑TA‑ATAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                    <  2:586927/149‑1 (MQ=255)                                                                                                      gATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  2:73289/1‑150 (MQ=255)                                                                                                       aTGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  2:553788/1‑149 (MQ=255)                                                                                                       aTGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  1:679262/1‑149 (MQ=255)                                                                                                         tccccgaacaccccAAGGAAAGATAATAAGT‑TTTCCCC‑TA‑ATAAC‑TC‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                               <  1:405668/139‑1 (MQ=255)                                                                                                          cccgaagaccccaAGGAAAAG‑TAATAAGT‑TTTCCCC‑TA‑AAAAT‑TC‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                           <  1:291683/142‑1 (MQ=255)                                                                                                           ccGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                          >  2:144497/1‑151 (MQ=255)                                                                                                           ccGAATACCCCCAGTTAAAA‑TAATAAGT‑TTTACCT‑TT‑ATAAA‑TA‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                             <  1:335713/148‑1 (MQ=255)                                                                                                             gCAAAGCCCCAGATAAAG‑AAATAGGG‑TTTCCCC‑AG‑ATAAT‑TA‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                             <  1:376568/146‑1 (MQ=255)                                                                                                                tAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                      >  1:191426/1‑150 (MQ=255)                                                                                                                atgCCCCAAGAAAAA‑TTAAAAGG‑TTTTCCC‑TT‑AAAAATTT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                      <  2:360870/149‑1 (MQ=255)                                                                                                                   cccaAGGAAAAG‑TAATAGGG‑TTTACCC‑TA‑AAAAT‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                    <  2:468840/149‑1 (MQ=255)                                                                                                                    cccAGTAAAAG‑AAAAAAGT‑TTTACCT‑TA‑ATAAA‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                 <  2:22524/151‑1 (MQ=255)                                                                                                                      cAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                >  1:504763/1‑150 (MQ=255)                                                                                                                      aggatAAAG‑AAATATGT‑TTTCCCT‑TA‑ATAAA‑TT‑CCCG‑CAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                <  1:431150/146‑1 (MQ=255)                                                                                                                       aGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               >  2:139444/1‑150 (MQ=255)                                                                                                                       aGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               >  1:16260/1‑150 (MQ=255)                                                                                                                       aGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               >  2:396022/1‑150 (MQ=255)                                                                                                                               agaaaaaGG‑TTTTCCCTTA‑AAAAA‑TC‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                        <  2:433671/146‑1 (MQ=255)                                                                                                                                taataaGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                         >  1:1038974/1‑148 (MQ=255)                                                                                                                                   taaGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcg                   >  1:504963/1‑151 (MQ=255)                                                                                                                                    aggg‑TTTCCCC‑TA‑AAAAT‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcg                   <  1:539733/148‑1 (MQ=255)                                                                                                                                       g‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcgt               >  1:570372/1‑151 (MQ=255)                                                                                                                                         tttccgc‑TG‑AAAAA‑TC‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                      <  2:1143141/137‑1 (MQ=255)                                                                                                                                         tttCCCC‑TA‑AAAAA‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                 <  2:135877/148‑1 (MQ=255)                                                                                                                                             tccccctaaaaa‑TT‑ACCGCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt              <  1:64665/141‑1 (MQ=255)                                                                                                                                             ccc‑TA‑AAAAA‑TC‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGa          <  1:196161/151‑1 (MQ=255)                                                                                                                                                    aTAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt      >  2:660675/1‑150 (MQ=255)                                                                                                                                                     tAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGCAATTAATAAGGCGCAGAGCGTCGTTTTGAACa       >  2:311549/1‑148 (MQ=255)                                                                                                                                                        a‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg  >  1:891636/1‑150 (MQ=255)                                                                                                                                                             |                                                                                                                                                    CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAA‑CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAG‑TAATAAGG‑TTTACCC‑TG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378653‑4378948

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