Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	44.0%	110.1 / 19.6	77	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (26/16);  new base C (0/33);  total (26/51)
Fisher's exact test for biased strand distribution p-value = 1.98e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTT  >  NC_000913/1620970‑1621267                                                                                                                                                      |                                                                                                                                                      tACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACAcc                                                                                                                                                     >  1:1499524/1‑151 (MQ=255)  aCTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACAccc                                                                                                                                                    >  1:1048311/1‑151 (MQ=255)  aCTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACAcc                                                                                                                                                     <  1:91514/150‑1 (MQ=255)      tCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCAt                                                                                                                                                 <  2:519591/150‑1 (MQ=255)       cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATc                                                                                                                                                <  1:429435/150‑1 (MQ=255)        ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                              >  1:548736/1‑151 (MQ=255)           cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                           <  2:777883/151‑1 (MQ=255)             tGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                         >  2:445629/1‑151 (MQ=255)              gCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTcgcg                                                                                                                                        <  2:751431/151‑1 (MQ=255)              gCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                         <  2:1435613/150‑1 (MQ=255)                  cAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACa                                                                                                                                     <  2:1484780/150‑1 (MQ=255)                     aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                                 >  1:899315/1‑151 (MQ=255)                       aaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTAc                                                                                                                               >  1:1437455/1‑151 (MQ=255)                        aGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTAcc                                                                                                                              <  1:1293570/151‑1 (MQ=255)                            aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCct                                                                                                                          >  1:368072/1‑151 (MQ=255)                                  aaaaGCGCCCAGCATAAAGCGAAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGC‑CCACTCCCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAAt                                                                                                                    <  1:1411100/151‑1 (MQ=255)                                  aaaaGCGACCAACATAAAGCCCAACCCGGCCAGCATCAGCGGTGGCATTTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAAt                                                                                                                    <  1:473903/151‑1 (MQ=255)                                        gACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGACCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCAt                                                                                                               >  2:1414911/1‑150 (MQ=255)                                           cAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCttt                                                                                                           >  2:1434151/1‑151 (MQ=255)                                            aGCAAAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGCCCACAAAATTTAGCCCCCATCCGC‑CCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtttt                                                                                                          <  2:929088/151‑1 (MQ=255)                                             acATAAAGCGCAACCCGCCCAGCTGCCCGGGTGGCTTTTTCTTAACCCACACTTTGTCGCCCAAAAATTTTACCCCCCACCCTC‑CCCCTCCCAGTAGCCCCACCACCCCATCTTTGCGCGACATCCTACCGCCTCTGAATTTCATCTTTTg                                                                                                         <  2:1532710/150‑1 (MQ=255)                                               aTAAAGCGCAACCCGCCCAGCCTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTCCCAGTAGCGCCAACACCCCATCTTTCCGCGACTTTCTCCCGCCTCTGAATTTCATCTTTTGTa                                                                                                       <  2:1052299/151‑1 (MQ=255)                                                 aaaGCGAAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGC‑CCACTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                        <  2:533915/148‑1 (MQ=255)                                                 aaaGCGAAACCCGGCCAGCATCAGCGGGGGCATGTTATAAAGCCCCATTTTGATGACCACAAAATTTACCCCCCATACGA‑CCCCTCCCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                       <  1:1126083/149‑1 (MQ=255)                                                 aaaCCGCAAGCGGGCAGGCATCAGCGGGGGCATGTTAGGAGCCCCCCTTTTGATACCCAAAAAATTTACCCCCCATACGA‑CCGCCCCGTGGAGCCCCACCCCCCCATCTTTTCCGCACATTCTACCGCCTCTGAATTTCAGCTTTTg                                                                                                         <  2:1241610/147‑1 (MQ=255)                                                      gcAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                                >  1:672905/1‑151 (MQ=255)                                                        aacccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                               <  1:47850/147‑1 (MQ=255)                                                        aaaCCGGCCAGCATCAGCGGTGGCATTTTATGAAGCCCCATTTTGATGACCACAAAATTTACCCCCCATACGA‑CCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                               <  1:503310/150‑1 (MQ=255)                                                             ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCCCTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTGCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                                         <  2:413101/151‑1 (MQ=255)                                                             ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                           >  1:1065104/1‑149 (MQ=255)                                                               ccAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                       >  2:1516152/1‑151 (MQ=255)                                                                cAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      >  2:1468480/1‑151 (MQ=255)                                                                 accATCAGCGGTGGCTTGTTATGAAGCCCCTTTTTAAGGACCAAAAAATTTAGCCCCCATCCGA‑CCCCTCCCAGGAGCGCCAACCCCCCATCTTTCCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGa                                                                                     <  2:1079334/149‑1 (MQ=255)                                                                 aGCATCAGGGTTGGCATGTTATGAACCCCCATTTTGATGACCAAAAAATTTAGCCCCCATCCGA‑CCCCTCCCAGTAGCGCAACCCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      <  2:9837/150‑1 (MQ=255)                                                                 aGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      >  1:264911/1‑150 (MQ=255)                                                                  gCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      >  1:451109/1‑149 (MQ=255)                                                                   cttCAGCGGTGGATTTTGAAAAAGCCCCATTTTGGTGCCTACAAATTTGACCCCCCATGCGA‑CCCCTACGAGTCGCGCCAACTCCCCTTGTTGTTGCGAAATTCTGCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                            <  2:557620/140‑1 (MQ=255)                                                                    aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGa                                                                                     <  1:400835/148‑1 (MQ=255)                                                                      cAGCGGTGGCATGTTATGAAGCCCCCTTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                                   <  2:1097155/148‑1 (MQ=255)                                                                      cAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATGCGA‑CCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAttt                                                                                 <  1:74138/150‑1 (MQ=255)                                                                      cAGCGGGGGCATGTTAGGAGCCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGC‑CCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      <  1:1120286/145‑1 (MQ=255)                                                                         cGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                             >  1:245700/1‑151 (MQ=255)                                                                           gTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                           >  2:1519211/1‑151 (MQ=255)                                                                            tGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTc                                                                          >  1:833878/1‑151 (MQ=255)                                                                             ggCATTTTATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGC‑CCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                         <  2:869810/151‑1 (MQ=255)                                                                                     aggaggCCCCATTTGGTTGACCACAAAATATGACCCCCA‑ATGACCCACGCCGGGTGACGCAAAAACCCCACTTCTTCGCAAATTACTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                    <  2:1149515/143‑1 (MQ=255)                                                                                     aTGAAGCCCCTCTTTGATGCCCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCCCCAACCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                                 <  2:1198320/151‑1 (MQ=255)                                                                                       gAAGCCCCATTTTGATGACCAAAAAATTTACCCCCCATACGC‑CCCCACCCAGTAGCCCCATCCCCCCTTTTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACa                                                                   <  2:1013877/147‑1 (MQ=255)                                                                                       gAAGCCCCACTTTGATGCCCAAAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCCCCAACGCCCCATCTTTTCGCCACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAg                                                                  <  1:1526628/148‑1 (MQ=255)                                                                                       gAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGC‑CCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGCTATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                                <  2:353686/150‑1 (MQ=255)                                                                                        aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCCCTACCAGTAGCTCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGa                                                              <  2:1525989/151‑1 (MQ=255)                                                                                                 tttGATGACCACAAAATTTAGCCCCCATCCGC‑CCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTa                                                       <  1:956001/149‑1 (MQ=255)                                                                                                 tttGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGt                                                     >  2:298067/1‑151 (MQ=255)                                                                                                  ttGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTc                                                    >  2:378405/1‑151 (MQ=255)                                                                                                    gatgaCCACAAAATTTAGCCCCCATCCGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGc                                                  <  1:857849/151‑1 (MQ=255)                                                                                                         ccACAAAATTTAGCCCCCATACGA‑CCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCg                                                <  1:115579/148‑1 (MQ=255)                                                                                                         ccACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                             >  1:150175/1‑151 (MQ=255)                                                                                                         ccACAAAATTAAGCCCCAATCCGC‑CCACTCCCAGTAGCCCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCgg                                               <  2:164387/149‑1 (MQ=255)                                                                                                         ccAAAAAATTTACCCCCCATCCGC‑CCCCTACAAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                             <  2:1101312/151‑1 (MQ=255)                                                                                                           acaAAATTTAGCCCCCATACGC‑CCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTa                                            <  1:351417/150‑1 (MQ=255)                                                                                                             aaaaTTTACCCCCCATACGC‑CCACTACCAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacg                                          <  2:68404/150‑1 (MQ=255)                                                                                                              aaaTTTAGCCCCCATCCGC‑CCACTACCATTAGCCCCAACCCCCCATCTTTTCGCGATTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                        <  2:18938/151‑1 (MQ=255)                                                                                                                 tttACCCCCCATACAC‑CCACTACCAGTAGCCCCAACCCCCCATCTTTTCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                      <  1:1498471/150‑1 (MQ=255)                                                                                                                     gCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                  >  1:1182139/1‑150 (MQ=255)                                                                                                                            atccgcccACTCCCAGGAGCGCCCACCCCCCATCTTTCCGCGACACTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                             <  2:565636/143‑1 (MQ=255)                                                                                                                             aCGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                          >  1:238325/1‑150 (MQ=255)                                                                                                                                  ccacacccAGTACCGCCAACCCCCCATCTTTTCGCAACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                       <  2:883729/143‑1 (MQ=255)                                                                                                                                  ccACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                     <  1:1502427/151‑1 (MQ=255)                                                                                                                                   cACTACCAGTAGCGCCACCACCCCTTCTTTTCGCGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAg                    <  2:377765/151‑1 (MQ=255)                                                                                                                                   cACTACCAGTAGCCCCAACCCCCCTTCTTTTCGCGACATTCTGCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                          <  2:938943/145‑1 (MQ=255)                                                                                                                                       aCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGc                >  2:189553/1‑151 (MQ=255)                                                                                                                                        ccAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc                  <  1:1493459/148‑1 (MQ=255)                                                                                                                                         cAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt               >  1:317530/1‑150 (MQ=255)                                                                                                                                               cgcCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACt          <  2:1524530/149‑1 (MQ=255)                                                                                                                                               cgcCAACCCCCCATCTTTCCCCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc        <  1:600041/151‑1 (MQ=255)                                                                                                                                                   aaCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGa     <  2:572670/150‑1 (MQ=255)                                                                                                                                                   aaCACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAt    >  2:145138/1‑151 (MQ=255)                                                                                                                                                    acacCCCATCTTTGCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATAc           <  2:548714/143‑1 (MQ=255)                                                                                                                                                     cacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAttt  >  2:208242/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGA‑CCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTT  >  NC_000913/1620970‑1621267

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: ‑