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

ACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTT  >  NC_000913/1620971‑1621267                                                                                                                                                     |                                                                                                                                                      aCTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAccc                                                                                                                                                    <  2:139805/151‑1 (MQ=255)    ttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCAt                                                                                                                                                 >  2:1279222/1‑151 (MQ=255)     tCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATc                                                                                                                                                >  1:1264710/1‑151 (MQ=255)       ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACGCCCCATCt                                                                                                                                               <  1:292448/150‑1 (MQ=255)            tGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                         >  1:71676/1‑151 (MQ=255)               gACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAc                                                                                                                                      >  2:313940/1‑151 (MQ=255)                aCAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAc                                                                                                                                      >  2:649869/1‑150 (MQ=255)                  aaaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAtt                                                                                                                                   >  2:136202/1‑151 (MQ=255)                   aaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTc                                                                                                                                  >  1:745874/1‑151 (MQ=255)                    aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                                 <  2:469197/151‑1 (MQ=255)                         aTAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGc                                                                                                                            >  2:265370/1‑151 (MQ=255)                          tAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGcc                                                                                                                           >  2:642225/1‑151 (MQ=255)                           aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGcc                                                                                                                           >  1:213092/1‑150 (MQ=255)                           aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCct                                                                                                                          <  1:705103/151‑1 (MQ=255)                           aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCct                                                                                                                          <  2:78532/151‑1 (MQ=255)                             ccGAAAAAGCGACAAACAAAAAGCGC‑AACCGGCCCGCCATCAGCGGTGGCATGTTATGAACCCCCATTTTGATGACCAAAAATTTTACCCCCCATCCGCCCCCCCCCAGTAGCCCAAACCCCCCATCTTTCCGCGACATTCTACCGCctct                                                                                                                        <  1:1247276/151‑1 (MQ=255)                             cagaaaaaGCGCCCAGCAAAAAGCGCAACCCGCCCAG‑CACCAGCGGTGGCATTTTATGAACCCCCCTTTTAATGACCAAAAATTTTACCCCCCATCCGCCCACTCCCAGTGGCGCAAACCCCCCATCTTTTCGCGACATTCTACCGCctct                                                                                                                        <  1:870118/147‑1 (MQ=255)                              cGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTg                                                                                                                       >  1:107536/1‑151 (MQ=255)                               ggAAAAGCGACCAGCATAAAGCGAAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAACCCCCCCTTTGATGACCAAAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGa                                                                                                                      <  2:309397/151‑1 (MQ=255)                                  aaggCGACCAGCATAAACCGAACACCGGCCAG‑CATCACGGGTGGCTTGTTAGAAACCCCCATTTTGAGCACCAAAAAATTTCCCCCCCATACCCCCCCCCCAGGAAGCGCCACCCCCCCACCTTTTCGCGACATTCTACCGCCTCTGAAtt                                                                                                                   <  1:951642/148‑1 (MQ=255)                                   aaGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAAttt                                                                                                                  >  1:815381/1‑151 (MQ=255)                                     gCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGCCCACAAAATTTAGCCCCCATACGCCCCCTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCa                                                                                                                <  2:1000943/151‑1 (MQ=255)                                         ccAGCATAAGGCGCAAACCGGCCAG‑CACCAGGGGTGGCATGTTTTGAAGCCCCCCTTTGATGACCACAAAATTTAGCCCCCATACGCCCCCTCCCAGTACCGCCACCACCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCt                                                                                                             <  2:1100734/150‑1 (MQ=255)                                         ccAGAATAAAGCGAAAACCGGCCAG‑CATCAGGGGTGGCATGTTATGAAGCCCCATTTTGATGACCACAAAATTTAGCCCCCATCCGACCCCTACCAGTAGCGCCAACACCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCt                                                                                                             <  2:421348/150‑1 (MQ=255)                                              aaaaaGCGCAAACCGGCCAG‑CATCAGGGGTGGCATGTTATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                        <  2:485417/148‑1 (MQ=255)                                               tAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                       >  1:304351/1‑150 (MQ=255)                                               tAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGCTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                       >  1:1039448/1‑150 (MQ=255)                                                aaaGCCAAACCCGCCCAG‑CATCAGCGGTGGCTTGTTAGGAACCCCCTCTTTGATTCCCACAAAATTTACCCCCCATACGACCCCTCCCATTAGCGCAAACCCCCCATCTTTTCCCGACTTTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                        <  2:48832/148‑1 (MQ=255)                                                    cgcAACCCGGCCAG‑CACCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCaa                                                                                                  <  2:411574/150‑1 (MQ=255)                                                     gcAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAt                                                                                                 >  1:671842/1‑150 (MQ=255)                                                           cGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                           >  1:1294324/1‑150 (MQ=255)                                                            ggCCAG‑CATCAGCGGCGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                           >  1:979512/1‑149 (MQ=255)                                                                   cATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                                   >  2:97942/1‑151 (MQ=255)                                                                    accAGCGGTGGCATGTTATGAACCCCCATTTTGATGCCCACAAAATTTACCCCCCATACGCCCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                        <  1:864182/143‑1 (MQ=255)                                                                          gggggctttttATGAACCCCCACTTTGATGCCCAAAAAATTTACCCCCCATACGCCCACTCCCAGTAGCCCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                            <  2:613136/142‑1 (MQ=255)                                                                          ggTGGCATGTTATGAACCCCCCCTTTGATACCCACAAAATTTACCCCCCATACGCCCCCTCCCAGTAGCCCAACCCCCCCACCTTTCCGCAACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                             <  1:613596/150‑1 (MQ=255)                                                                           ggggCATTTTATGAACCCCCCCTTTGATGACCACAAAATTAACCCCCAAAACAACCACTACCAGTAGCCCAACCCCCCCATCTTTGCGCGACTTTTTCCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                           <  2:124010/149‑1 (MQ=255)                                                                           gTGGCTTTTTAGGAAACCCCACTTTGATGCCCAAAAAATTTACCCCCCATCCTCCCCCCACCAGTAGCGCCACCCCCCATTCTTTCCCCGACATTTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                           <  1:977938/151‑1 (MQ=255)                                                                             ggCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                         <  2:924079/151‑1 (MQ=255)                                                                                 tGTTATGAACCCCCACTTTGATGACCACAAAATTTACCCCCCATACGCCCACTCCCAGTAGCGCCACCCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                         <  1:13535/147‑1 (MQ=255)                                                                                   ttATGAAGCCCCACTTTGATGCCCACAAAATTTACCCCCCATACGCCCCCTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTaa                                                                     <  1:19740/149‑1 (MQ=255)                                                                                   ttATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                    <  2:1041884/150‑1 (MQ=255)                                                                                   ttATAAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                    <  1:526359/150‑1 (MQ=255)                                                                                    tATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACa                                                                   <  2:152343/150‑1 (MQ=255)                                                                                        acgCCCCATTTTGAAGCCCCAAAAGTTTACCCCCCATAGGACCCCTTCCAGAGGCGCCAACTCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGa                                                              <  2:815413/149‑1 (MQ=255)                                                                                        aaccccccttttGATGCCCACAAAATTTACCCCCCATACGCCCCCTCCCATTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                                <  2:619668/140‑1 (MQ=255)                                                                                        aaGCCCCATTTTGAGAACCACAAAATTAAGCCCCCAACCGCCCACTCCCAGTAGCGCCACCCCCCCATCTTTTCGCAACATCCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTaa                                                                     <  1:532206/144‑1 (MQ=255)                                                                                         aGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAt                                                             >  1:619485/1‑151 (MQ=255)                                                                                            cccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                          >  2:967210/1‑151 (MQ=255)                                                                                              cACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGtt                                                        >  1:1120727/1‑151 (MQ=255)                                                                                                 tttGATGACCAAAAAATTTAGCCCCCATACGCCCACTCCCAGTGGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTa                                                       <  2:961550/149‑1 (MQ=255)                                                                                                  ttGATGACCACAAAATTTACCCCCAATCCGCCCACTACCAGTAGCGCCAACACCCCATCTTTCCCCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGt                                                     <  1:517697/150‑1 (MQ=255)                                                                                                    gatgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGc                                                  >  1:1244532/1‑151 (MQ=255)                                                                                                         ccACAAATTTTAGCCCCCATACGACCACACCCAGTGGCGCAACCACCCTATCTTTCCGCGACATTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                             <  2:182191/151‑1 (MQ=255)                                                                                                         ccACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                             >  1:562209/1‑151 (MQ=255)                                                                                                           acaAAATTTACCCCCCATACGCCCACTCCCAGTAGCCCAAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCg                                                <  1:1279152/146‑1 (MQ=255)                                                                                                           aaaaatTTTACCCCCCATACGACCACTCCCAGTCGCGCCAACCCCCCATCTTTCCGCAACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGt                                              <  1:1064725/146‑1 (MQ=255)                                                                                                             aaatTTTAGCCCCAATACGCCCACTCCAAAAAGCGCCACCACCCCATATTTTCGCGATATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCAGGTTac                                           <  1:230625/149‑1 (MQ=255)                                                                                                              aaaTTTACCCCCCATACCCCCACTACCAGCAGCCCCAACCCCCCTTCTTTGCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                        <  1:124021/151‑1 (MQ=255)                                                                                                               aaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                       >  2:124412/1‑151 (MQ=255)                                                                                                               aaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                       >  2:505/1‑151 (MQ=255)                                                                                                                 tttAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACTGTTGATTCGTTAGTCGCCGGTTACGACGGc                                     >  2:724914/1‑151 (MQ=255)                                                                                                                        cccATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATg                              <  2:1175480/151‑1 (MQ=255)                                                                                                                         ccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                             >  1:1253064/1‑151 (MQ=255)                                                                                                                         ccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                             >  2:1150726/1‑151 (MQ=255)                                                                                                                              cGACCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaaa                        <  2:1208685/151‑1 (MQ=255)                                                                                                                                 ccACTACCGGTAGCGACACCCCCCCATCTTCTCGATCCTTTATACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                       <  1:121504/149‑1 (MQ=255)                                                                                                                                     tACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg                 >  2:998983/1‑151 (MQ=255)                                                                                                                                       ccAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt               >  2:1081980/1‑151 (MQ=255)                                                                                                                                        cAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCta              <  2:642748/151‑1 (MQ=255)                                                                                                                                         aGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtat             >  1:135719/1‑151 (MQ=255)                                                                                                                                            aGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACt          >  1:206087/1‑151 (MQ=255)                                                                                                                                              cgcCACCACCCCTTCTTTCCGCAACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc        <  1:791820/151‑1 (MQ=255)                                                                                                                                                ccAACCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCg       <  2:235612/150‑1 (MQ=255)                                                                                                                                                ccAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg      >  2:206697/1‑151 (MQ=255)                                                                                                                                                 cAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg      >  2:1061981/1‑150 (MQ=255)                                                                                                                                                   acacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtt   >  1:1295021/1‑151 (MQ=255)                                                                                                                                                   acacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtt   >  2:1234205/1‑151 (MQ=255)                                                                                                                                                    cacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAttt  >  1:242759/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAG‑CATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTT  >  NC_000913/1620971‑1621267

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