Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	49.3%	82.5 / 31.5	76	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (29/9);  new base C (0/37);  total (29/47)
Fisher's exact test for biased strand distribution p-value = 3.95e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGA  >  NC_000913/1620968‑1621264                                                                                                                                                        |                                                                                                                                                   ggTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAaca                                                                                                                                                    <  2:1151844/151‑1 (MQ=255)  gTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAacac                                                                                                                                                   >  2:286226/1‑151 (MQ=255)    aCTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACAccc                                                                                                                                                 >  2:262752/1‑151 (MQ=255)         cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCt                                                                                                                                            <  1:1067797/151‑1 (MQ=255)          ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                           >  2:525143/1‑151 (MQ=255)              gTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcg                                                                                                                                       <  1:184814/151‑1 (MQ=255)                gCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                      <  2:694702/150‑1 (MQ=255)                    cAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAt                                                                                                                                 >  1:815479/1‑151 (MQ=255)                        aaaTATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                              >  2:1145211/1‑150 (MQ=255)                        aaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTa                                                                                                                             <  1:409930/151‑1 (MQ=255)                            aTAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGc                                                                                                                         >  1:1051628/1‑151 (MQ=255)                            aTAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGc                                                                                                                         <  2:55911/151‑1 (MQ=255)                              aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCct                                                                                                                       >  1:100606/1‑151 (MQ=255)                                 cGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTg                                                                                                                    >  2:1104672/1‑151 (MQ=255)                                       aGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTc                                                                                                              >  2:244652/1‑151 (MQ=255)                                              aGCAAAAAGCGCAAACCGCCCAGCATCAGCGGTGGCATGTTATGAAGCCCCATTTTGATGACCACAAAATTTAGCCCCCAC‑CCGACCACTCCCAGTAGCGCCACCCCCCCTTCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtt                                                                                                         <  1:38222/149‑1 (MQ=255)                                                 aTAAAGCGCAACCCGGCCAGCATCAGCGGTGGCTTTTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                     <  2:241422/150‑1 (MQ=255)                                                 aTAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGCCCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                    <  2:1065793/151‑1 (MQ=255)                                                 aTAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                     >  2:1256001/1‑150 (MQ=255)                                                  tAAAGCGAAACCCGGCCAGCATCAGCGGGGGCATTTTATGAAGCCCCCCTTTGATGCCCCCAAAATTTAGCCCCCAT‑CCGACCCCTACCAGTAGCCCCAACCCCCCATCTTTGCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTg                                                                                                      <  2:23778/148‑1 (MQ=255)                                                   aaaGCCCAACCCGGCCAGCACCAGCGGTGGCTTGTTATGAACCCCCATTTTGATGCCAAAAAATTTTACCCCCAAT‑CCGCCCCCTCCAAGTGGCGCCAACCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                     <  1:851297/148‑1 (MQ=255)                                                   aaaCCGCAACCCGGCCCGCATCAGGGTTGGTTTGTTAAGAAGCCCCACTTTGAGGACAAAAAAATTTAGCCCCCAA‑ACGCCCCCACCCAGTGGCGCCAACCCCCTATCTTTGCGCGACATTCATCCGCCTCTGAATTTCATCTTTTGTAAg                                                                                                  <  2:988615/151‑1 (MQ=255)                                                    aaGCGCAACCCGGCCAGCACCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCAT‑CCGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                 <  2:343569/151‑1 (MQ=255)                                                        gcAACCCGGCCAGCCTCACCGGGCGAATGTTTTGACCCCCCTTTTTGTTGCCCACAAATTTAAGCCCCCAG‑CCGACCGCTCCCAGTGGCCCCAATCCCCCTTCCTTACGCCCCATACTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                                <  2:1296233/148‑1 (MQ=255)                                                          aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                            >  1:152021/1‑150 (MQ=255)                                                           aaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAAc                                                                                          >  1:1210513/1‑151 (MQ=255)                                                             ccGGCCAGCATCAGGGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                        <  2:935669/151‑1 (MQ=255)                                                             ccGGCCAGCATCAGCGGTGCCATTTTATGAACCCCCTCTTTGATGACCACAAAATTTACCCCCCAT‑ACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                        <  1:571896/151‑1 (MQ=255)                                                             ccGGCCAGCATCAGCGGGGGCATGTTATGAAGCCCCTTTTTGTTGCCCACAAAATTTACCCCCCAT‑CCACCCCCTCCCAGTAGCCCCAACCCCCCACCTTTTCGCAACATTCTGCCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                             <  1:1081347/146‑1 (MQ=255)                                                              cGGCCAGCACCAGCGGTGGCATGTTATAAACCCCCATTTTAATGCCCAAAAAATTAACCCCCAAT‑ACGCCCCCTCCCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTa                                                                                       <  1:425205/151‑1 (MQ=255)                                                                 ccAGCATCAGCGGTGGCTTGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                     <  1:17100/150‑1 (MQ=255)                                                                   aGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                   >  1:1123776/1‑150 (MQ=255)                                                                      aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAtt                                                                               >  2:309672/1‑151 (MQ=255)                                                                          ggggTGGAATTTTATGAAGCCCCTCTTTGATGCCCAAAAAATTTAGCCCCCAT‑ACGCCCCCTCCCAGTAGCGCCACCCCCCCATCTTTCCGCGATTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTa                                                                             <  1:527601/147‑1 (MQ=255)                                                                            gggggCATTTTATGAACCCCCCCTTTGTTACCCACAAATTTTAGCCCCCAT‑CCGCCCACTCCCAGTAGCGCCACCTCCCCTTTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTa                                                                             <  1:548352/144‑1 (MQ=255)                                                                             gTGGCATTTTATGAACCCCCACTTTGATGACCAAAAAATTTACCCCCCAT‑ACGACCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                        <  2:1047545/151‑1 (MQ=255)                                                                                   ttttttaatgCCCCTCTTTATTGCCCACAAATTTTAGCCCCCAT‑ACGCCCACTCCCGGTAGCCCTAACCCCCCATTTTTTCGATTCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTc                                                                       <  1:1182818/137‑1 (MQ=255)                                                                                    gTTATCAACCCCCCCTTAACTGACCACAAAA‑TTATCCCCCGTGACACCCCATCCCAGTAGCCCCTACACCCCACCTTTTCGCTAAATTGTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATGTACGTTTCTTTaa                                                                  <  1:190537/150‑1 (MQ=255)                                                                                     ttaggaccccccACTTTGATGCCCAAAAATTTTACCCCCCAT‑CCCCCCACTCCCATTGGCGCCAACCCCCCATCTTTGCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTa                                                                   <  2:122154/140‑1 (MQ=255)                                                                                          aaGCCCCCTTTTGAAGACCAAAAAATTTAGCCCCCAT‑ACGCCCACTACAAGTAGCCCCAACCCCCCATCTTTCCGCGATATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                              <  2:820358/148‑1 (MQ=255)                                                                                          aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                            >  1:353021/1‑150 (MQ=255)                                                                                            gCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                         >  1:1093070/1‑151 (MQ=255)                                                                                            cccccATTTTGATGACCACAAAATTTAGCCCCAAT‑ACGCCCACTCCCAGTAGCCACAACCCCCCATTTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                             <  1:1211795/146‑1 (MQ=255)                                                                                             ccccgcttgaaggaCCACAAAATTTACCCCCCAT‑CCACCCCCAACAAGTACCGCCACCACCCCATCTTTCCGAGCTTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                            <  1:43023/135‑1 (MQ=255)                                                                                             ccccACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                         <  2:763552/150‑1 (MQ=255)                                                                                              cccACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                       >  2:376630/1‑151 (MQ=255)                                                                                              cccACTTTGAGGACCAAAAATTTTAGCCCACAT‑CCGAAAAAAACAGGAGGCGCTAACCCCCCTTTTTTACGCAAAATTCTTCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                         <  1:882862/149‑1 (MQ=255)                                                                                               ccACTTTAAGGACCAAAAAATTTACCCCCCAT‑ACGACCACTCCAAGTAGCGCAAACCCCCCATCTTTTCCCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGt                                                      <  2:613271/151‑1 (MQ=255)                                                                                                   tttGATGCCCACAAAATTAAGCCCCAGT‑CCGACCACTACAAGTCGCGCCCGCACCCAATTTTTTCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTa                                                    <  2:275107/149‑1 (MQ=255)                                                                                                   tttGATGACCACAAAATTTACCCCCCAT‑ACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTa                                                    <  2:1283569/149‑1 (MQ=255)                                                                                                     tgatgaCCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCg                                                >  1:1203823/1‑151 (MQ=255)                                                                                                          aCCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTGGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGt                                           >  1:1254881/1‑151 (MQ=255)                                                                                                           ccACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                          >  2:292178/1‑151 (MQ=255)                                                                                                           ccACAAAATTTACCCCCCAT‑ACGCCCACTCCCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                          <  1:1247578/151‑1 (MQ=255)                                                                                                           ccAAAAATTTTAGCCCCCAT‑CCCCCCACTCCCATTAGCCCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCgg                                            <  1:1279189/149‑1 (MQ=255)                                                                                                            caaaaaaTTTACCCCCCAT‑ACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTa                                         <  2:1226003/148‑1 (MQ=255)                                                                                                               aaaaTTTGCCCCCCAT‑CCGACCACTCCCATTACCCCCAACCCCCCTTCTTTCCGCTACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTa                                         <  1:1179199/148‑1 (MQ=255)                                                                                                               aaaaTTTAGCCCCCAT‑ACGCCCACTCCCAGTAGCCCCAACCCCCCATCTTTCCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacga                                      <  1:244650/151‑1 (MQ=255)                                                                                                               aaaaTTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                        >  1:369490/1‑149 (MQ=255)                                                                                                                 aaTTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                    >  1:1161856/1‑151 (MQ=255)                                                                                                                 aaTTCAGCCCCCAC‑ACCCCCACTACCAGAGGCCCAACCCCCCCATCTTTTCGAAACATTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                    <  1:1001965/151‑1 (MQ=255)                                                                                                                    ttAGCCCCCAT‑CCGCCCACTACCAGTGGCGCAAACCCCCCATCTTTTCCCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                  <  1:885502/150‑1 (MQ=255)                                                                                                                     tAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAt                                >  2:534691/1‑151 (MQ=255)                                                                                                                      aGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                               >  2:677916/1‑151 (MQ=255)                                                                                                                        cccccAT‑CCGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTCCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                               <  1:1145139/149‑1 (MQ=255)                                                                                                                          cccAT‑ACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATg                           <  2:1016857/151‑1 (MQ=255)                                                                                                                           ccAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                          >  2:978458/1‑151 (MQ=255)                                                                                                                                        tACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg              >  2:463298/1‑151 (MQ=255)                                                                                                                                         cccAGTAGCCCCACCCCCCCTTTTTTTCGCGACTTTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGc             <  1:525128/150‑1 (MQ=255)                                                                                                                                          ccAGTCGCCCCAACCCCCCATCTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg              <  1:618055/149‑1 (MQ=255)                                                                                                                                          ccAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt            <  2:836779/151‑1 (MQ=255)                                                                                                                                           cAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt            >  2:102771/1‑150 (MQ=255)                                                                                                                                            aGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtat          >  2:198961/1‑151 (MQ=255)                                                                                                                                                 cgcCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATAc        <  2:448636/148‑1 (MQ=255)                                                                                                                                                  gcCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCg    >  2:1211171/1‑151 (MQ=255)                                                                                                                                                  gcCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCg    >  2:126193/1‑151 (MQ=255)                                                                                                                                                   ccAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtata         <  1:1236296/145‑1 (MQ=255)                                                                                                                                                     aaCCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTGACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGa  <  1:483555/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   GGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAT‑ACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGA  >  NC_000913/1620968‑1621264

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