Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	42.5%	121.7 / 21.3	80	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (31/15);  new base C (0/34);  total (31/49)
Fisher's exact test for biased strand distribution p-value = 3.72e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGG  >  NC_000913/1620968‑1621263                                                                                                                                                       |                                                                                                                                                   ggTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAaca                                                                                                                                                    <  1:961065/151‑1 (MQ=255) ggTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAaca                                                                                                                                                    <  2:36629/151‑1 (MQ=255)    aCTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAccc                                                                                                                                                 >  2:702462/1‑151 (MQ=255)      tttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcccc                                                                                                                                                >  2:700221/1‑150 (MQ=255)       ttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcccc                                                                                                                                                <  1:880305/149‑1 (MQ=255)         cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATc                                                                                                                                             <  2:244542/150‑1 (MQ=255)          ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                           <  1:1184719/151‑1 (MQ=255)          ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                           >  2:671227/1‑151 (MQ=255)          ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                           >  2:32347/1‑151 (MQ=255)           gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                          >  1:1008030/1‑151 (MQ=255)             cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                        >  1:1111064/1‑151 (MQ=255)                  gACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAc                                                                                                                                   >  2:588109/1‑151 (MQ=255)                     aaaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑Att                                                                                                                               <  1:996589/151‑1 (MQ=255)                        aaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTa                                                                                                                            <  2:506915/151‑1 (MQ=255)                              aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGcc                                                                                                                       <  2:191562/150‑1 (MQ=255)                                ccGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCctct                                                                                                                    >  1:295284/1‑151 (MQ=255)                                  ggAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTGGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGa                                                                                                                  <  2:1072877/151‑1 (MQ=255)                                  ggAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTg                                                                                                                   >  2:221449/1‑150 (MQ=255)                                   gAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGa                                                                                                                  >  2:284162/1‑150 (MQ=255)                                   gAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGa                                                                                                                  >  1:28892/1‑150 (MQ=255)                                   gAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGa                                                                                                                  >  1:423297/1‑150 (MQ=255)                                    aaaaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAAt                                                                                                                >  2:1232212/1‑151 (MQ=255)                                      aaGCGCCCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAAttt                                                                                                              <  2:1242359/151‑1 (MQ=255)                                      aaGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAAttt                                                                                                              <  2:702410/151‑1 (MQ=255)                                           aCCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCt                                                                                                         >  1:1034131/1‑151 (MQ=255)                                             cAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCttt                                                                                                       >  2:597111/1‑151 (MQ=255)                                              aGCAAAAGGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCCCTTTGATACCCAAAAATTTTAGCCCCCATACGCCCGCTACCAGTAGCGCCAACCCCCCTTCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCtttt                                                                                                      <  2:767206/151‑1 (MQ=255)                                                 aTAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCATTTTGATGACCACAAAATTTAGCCCCCATACGACCCCTCCAAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTg                                                                                                     <  2:1051898/149‑1 (MQ=255)                                                   aaaGCCAAACCCGGCCAGCATCAGCGGTGGCATTTTATGAAGCCCCACTTTGATACCCACAAAATTTAGCCCCCATACGCCCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAg                                                                                                 <  1:231170/151‑1 (MQ=255)                                                    aaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                >  1:602172/1‑151 (MQ=255)                                                         cAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                            >  2:82044/1‑150 (MQ=255)                                                          aacccGGCCAGCATCAGGGGTGGCTTTTTATGAAGCCCCATTTTGATGACCACAAAATTTACCCCCCATACGCCCACCCCCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                               <  1:542422/143‑1 (MQ=255)                                                          aacccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCaa                                                                                          <  2:376476/148‑1 (MQ=255)                                                          aacccGGCAAGCATCAGCGGTGGCTTTTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAAACGCCCACTACCAGTAGCGCCACCCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                           <  1:94796/147‑1 (MQ=255)                                                          aaaCCGGCCAGCATCAGGGGTGGCATGTTATGAAGCCCCCCTTTGATGACCACAAAATTTAGCCCCCATACGCCCCCTACCAGTAGCGCCAACCCCCCATCTTTCCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                               <  1:1174282/146‑1 (MQ=255)                                                             ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCTCTTTGATGCCCACAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCCCCAACCCCCCATCTTTTCCCGAC‑TTCCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                        <  1:1239216/150‑1 (MQ=255)                                                             ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                        <  1:229962/150‑1 (MQ=255)                                                             ccGCCCAGCATCAGCGGTGGCATTTTATGAACCCCCATTTTGATGACCACAAAATTTACCCCCCATCCGCCCACTACCAGTAGCCCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAAc                                                                                         <  2:1150488/149‑1 (MQ=255)                                                                 ccAGCATCAGCGGTGGCATGTTATGAACCCCCACTTTGATAACCACAAAATTTAGCCCCCATCCGCCCCCTCCCAGAAGCGCCAACCCCCCATCTTTTCGCGAA‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                   <  2:646854/151‑1 (MQ=255)                                                                  cAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                  >  1:444674/1‑151 (MQ=255)                                                                          gCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACt                                                                          >  1:951674/1‑151 (MQ=255)                                                                            gggggCAGGTAATGAAGCCCCACTTTGATAACCACAAAATTTAGCCCCCATGGGACCCCTCCTGGTAGCGCCACCCCCCCATCTTTTCGCGAC‑ATTCTGCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                         <  1:516194/147‑1 (MQ=255)                                                                            ggTGGCATGTTATGACCCCCCACTTTGATGACCAAAAAATTTACCCCCCATACGCCCCCTACCAGTACCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                        <  1:853456/151‑1 (MQ=255)                                                                                 cATGTTATGAACCCCCACTTTGATGCCCACAAAATTTAGCCCCCATCCGACCACTCCCAGTAGCCCCAACCCCCCATCTTTTCGCGAC‑TTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCtt                                                                    <  1:470178/150‑1 (MQ=255)                                                                                  aTGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGCCCACTACCAGTAGCCCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                   <  2:1004362/150‑1 (MQ=255)                                                                                   tGTTATGAAGCCCCACTTTGATGACCAAAAAATTTACCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTa                                                                  <  2:549652/150‑1 (MQ=255)                                                                                    gTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                >  1:710963/1‑151 (MQ=255)                                                                                     ttATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                >  1:243774/1‑150 (MQ=255)                                                                                     ttATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                >  1:464065/1‑150 (MQ=255)                                                                                     ttATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                >  1:554998/1‑150 (MQ=255)                                                                                         gAAGCCCCCTTTTGATGACCAAAAAATTTAGCCCCCATACGACCCCTACCAGTAGCCCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                           <  2:128615/151‑1 (MQ=255)                                                                                          aacccccccTTTGATGACCAAAAAATTTACCCCCAATCCGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                             <  1:852025/145‑1 (MQ=255)                                                                                          aacccccccTTTAATGACCAAAAAATTTAGCCCCCATACGACCCCTCCCATTAGCCCCAACCCCCCATCTTTTCGCGAC‑TTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                            <  2:346178/146‑1 (MQ=255)                                                                                             ccccACTTTGATGCCCAAAAAATTTAGCCCCCATCCGCCCCCTACCAGTAGCGCCAACCCCCCACCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                        <  1:231501/150‑1 (MQ=255)                                                                                             ccccACTTTGATGACCAAAAAATTTACCCCCCATCCGCCCCCTCCCAGTAGCGCCACCCCCCCATCTTTCCGCGAC‑TTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                        <  1:26285/150‑1 (MQ=255)                                                                                              cccACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAttag                                                      <  2:1092814/151‑3 (MQ=255)                                                                                               ccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGt                                                     >  1:437754/1‑151 (MQ=255)                                                                                                       atgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTc                                                <  2:371544/148‑1 (MQ=255)                                                                                                       atgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGcc                                             <  2:1035362/151‑1 (MQ=255)                                                                                                         gaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCgg                                           >  1:220764/1‑151 (MQ=255)                                                                                                             acaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                       >  1:61523/1‑151 (MQ=255)                                                                                                             acaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                       >  1:1171906/1‑151 (MQ=255)                                                                                                                aaaTTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                    <  2:1124622/151‑1 (MQ=255)                                                                                                                 aaTTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                   <  2:1178722/151‑1 (MQ=255)                                                                                                                  aTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                  >  1:492153/1‑151 (MQ=255)                                                                                                                   tttAGCCCCCATACGCCCCCTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                  <  1:47321/150‑1 (MQ=255)                                                                                                                   tttAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                  <  2:410483/150‑1 (MQ=255)                                                                                                                    ttAGCCCCCATTCTCCCACGACCAGCACCCCCACCCCCCCATTTTATC‑CAACTATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                 <  1:381148/150‑1 (MQ=255)                                                                                                                        cccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTa                             >  2:403587/1‑150 (MQ=255)                                                                                                                        cccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTa                             <  2:475132/150‑1 (MQ=255)                                                                                                                        cccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTa                             <  2:1002503/150‑1 (MQ=255)                                                                                                                                   cccctcccAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                   <  1:567740/143‑1 (MQ=255)                                                                                                                                   ccACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                 <  1:906321/151‑1 (MQ=255)                                                                                                                                   ccACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa                  >  2:483853/1‑150 (MQ=255)                                                                                                                                    cccTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAg                <  2:461319/149‑1 (MQ=255)                                                                                                                                       tcccATTAGCCCCAACCCCCCATCTTTTCGCGAA‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGt               <  1:939249/147‑1 (MQ=255)                                                                                                                                         ccAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑TTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGt               <  1:463499/147‑1 (MQ=255)                                                                                                                                         ccAGTAGCGCCAACCCCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc              <  2:654148/148‑1 (MQ=255)                                                                                                                                         ccAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGc            <  1:344810/150‑1 (MQ=255)                                                                                                                                                cgcCAACCCCCCATCTTTTCGCGAA‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCta          <  1:284161/145‑1 (MQ=255)                                                                                                                                                cgcCAACCCCCCATCTTTCCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc    <  2:1190928/151‑1 (MQ=255)                                                                                                                                                   cAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg  >  1:380069/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                   GGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAC‑ATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGG  >  NC_000913/1620968‑1621263

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