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

TACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGG  >  NC_000913/1620970‑1621263                                                                                                                                                     |                                                                                                                                                  tACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcc                                                                                                                                                 >  1:692658/1‑151 (MQ=255)       cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATc                                                                                                                                            <  2:278164/150‑1 (MQ=255)       cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCt                                                                                                                                           <  2:1152521/151‑1 (MQ=255)       cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCACCATACGACCACTACCAGTAGCGCCAACACCCCATc                                                                                                                                            >  2:498947/1‑150 (MQ=255)        ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTGTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtt                                                                                                                                          <  2:531494/151‑1 (MQ=255)         gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                         <  2:947894/151‑1 (MQ=255)         gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                         >  1:689395/1‑151 (MQ=255)           cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtttt                                                                                                                                        <  1:866707/150‑1 (MQ=255)           cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                       <  2:44435/151‑1 (MQ=255)            gTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                       >  2:30605/1‑150 (MQ=255)              gCGACAAAAAAGATAGCCGGAAAAGCGTCCAGCGTAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTACGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCTGTAGCGCCAACTCCCCATCTTTTcgc                                                                                                                                     <  2:339982/150‑1 (MQ=255)                    aaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTc                                                                                                                              >  2:286483/1‑151 (MQ=255)                         gATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCg                                                                                                                         <  2:26685/151‑1 (MQ=255)                           tAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGcc                                                                                                                       >  1:283047/1‑151 (MQ=255)                            aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCct                                                                                                                      >  1:661799/1‑151 (MQ=255)                                ggAAAGGCGACCAGCATAAAGCGAAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGTCCAAAAAATTTACCCCCCATACGACCACTACCAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGa                                                                                                                  <  2:1183570/151‑1 (MQ=255)                                ggAAAAGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGa                                                                                                                  <  1:460783/151‑1 (MQ=255)                                   aaaGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAAtt                                                                                                               <  1:85123/151‑1 (MQ=255)                                        gACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATc                                                                                                          >  1:862776/1‑151 (MQ=255)                                               aaaaaGCGAAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACAACAAAATTTACCCCCCATACGCCCCCTACAAGTAGCGCAACCCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                   <  2:716523/149‑1 (MQ=255)                                                  aaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                >  2:1086022/1‑151 (MQ=255)                                                    gcgcAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                               <  1:692051/150‑1 (MQ=255)                                                        aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                           >  2:1011335/1‑150 (MQ=255)                                                        aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                           >  1:1372588/1‑150 (MQ=255)                                                          aCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACACAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                        >  2:463064/1‑151 (MQ=255)                                                           ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                        <  1:1021698/150‑1 (MQ=255)                                                             ggcacgCATCAGTGTGGTCATGTTATGAAGCCCCACTTTGATCACCACAAAATTTACCCCCTACCCACCCTCTCCAATTACCGCTAACCCCCCTTTTTTTCACGCAATCTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                                     <  2:1428366/146‑1 (MQ=255)                                                             ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCATTTTGATAACCACAAAATTTAGCCCCCATCCGCCCACTACCAGTAGCGCAACCCCCCCTTCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                                     <  1:812360/151‑1 (MQ=255)                                                             ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                       >  2:269783/1‑149 (MQ=255)                                                               ccAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                   <  1:15553/151‑1 (MQ=255)                                                               ccAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                    <  1:1158783/150‑1 (MQ=255)                                                               ccAGAACCATCGGTGCCTTGTTATACAGCCCCTTTTTATTGCCCACAAAATTTACCCCCCATACGCCCACCCCCAGTAGCTCAACCCCCCCATCTTTCCTCTACCTTTTGCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                   <  2:156458/151‑1 (MQ=255)                                                                 agcacaaGCGGTGGCAGGTTAGGAACCCCCACTTTGAGGACCACAAAATTTAGCCCCCATCCGACCCCTCCCAGGAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                    <  2:637780/142‑1 (MQ=255)                                                                 aGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTGGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGa                                                                                 <  1:92688/151‑1 (MQ=255)                                                                  gCAACAGGGGTGCCATGTTATGAACCCCCATTTTGATACCCACAAAATTTAGCCCCCATACGACCATTCCAATTAGCGCAAACCCCCCTTCTTTTCGAGATATTCGCCCGCCTCTGAATTTCATCTTTGGTAAGCAATCAACTTAGCTg                                                                                  <  2:1138872/149‑1 (MQ=255)                                                                             ggCATGTTATGACCCCCCATTTTGATGACCACAAAATTTACCCCCCATACGACCACTCCCAGTGGCCCCAACCCCCCTTTTTTGCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                       <  2:27954/149‑1 (MQ=255)                                                                              gCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCtt                                                                    <  1:101832/151‑1 (MQ=255)                                                                                aTGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                   >  2:1233679/1‑150 (MQ=255)                                                                                    tATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                >  2:245202/1‑149 (MQ=255)                                                                                    tATGAAGCCCCACTTTGATGACAACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTa                                                                  >  2:1273401/1‑147 (MQ=255)                                                                                        aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                            <  1:455701/149‑1 (MQ=255)                                                                                        aaGCCCCACTTTGATGACCACAAAATTTACCCCCCATACCCCCCCCCCCATTGGCGCCAACCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                           <  2:243599/150‑1 (MQ=255)                                                                                          gCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                        >  2:31157/1‑151 (MQ=255)                                                                                          cccccACTTTGATGACCACAAAATTTAGCCCCCATACGCCCCCTCCCAGTAGCGCCACCCCCCCATCTTTCCACGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                        <  2:369487/150‑1 (MQ=255)                                                                                           ccccccTTTGATGACCACAAAATTTAGCCCCCATACGCCCCCTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAt                                                         <  1:1119174/149‑1 (MQ=255)                                                                                           ccccACTTTAATGACCACAAAATTTACCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                        <  1:1023117/150‑1 (MQ=255)                                                                                             ccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTGTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGt                                                     <  2:1106737/151‑1 (MQ=255)                                                                                              cACTTTGATAACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCCCCACCCCCCCTGTTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCttt                                                    <  2:488954/151‑1 (MQ=255)                                                                                                  ttGATGACCAAAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGt                                                 <  1:1159364/150‑1 (MQ=255)                                                                                                   tgatgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCg                                               >  1:434648/1‑151 (MQ=255)                                                                                                     atgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGcc                                             >  1:208877/1‑151 (MQ=255)                                                                                                         ccACAAAATTTACCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCTTCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGt                                          <  2:1023183/150‑1 (MQ=255)                                                                                                          cacaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGACTTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                         >  2:89202/1‑150 (MQ=255)                                                                                                           acaAAATTTACCCCCCAAACGACCACTCCCATTAGCGCCAACCCCCCATTTTTCCGCAACATGCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAGCTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                         <  2:467066/149‑1 (MQ=255)                                                                                                              aaaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                    >  2:595175/1‑151 (MQ=255)                                                                                                              aaaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                    >  1:254792/1‑151 (MQ=255)                                                                                                               aaTTTAGCCCCCACACGCCCCCTCCCAGTAGCGCTAACCCCCCATCTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                   <  2:286037/151‑1 (MQ=255)                                                                                                                aTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                  >  2:200249/1‑151 (MQ=255)                                                                                                                  ttAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                 >  1:513506/1‑150 (MQ=255)                                                                                                                    aGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                              >  2:255545/1‑151 (MQ=255)                                                                                                                      cccccATACGCCCACTTCCAGTAGCGCCAACACCCCATCTTTACGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTaa                            <  1:1176505/151‑1 (MQ=255)                                                                                                                       ccccATACGCCCTCGTCCAGTGGCGCAAACACCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTGCGACGGCATTAAt                           <  2:850465/151‑1 (MQ=255)                                                                                                                       ccccATACGACCCCTACCAGTGGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                           <  1:269847/151‑1 (MQ=255)                                                                                                                       ccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                           >  2:136685/1‑151 (MQ=255)                                                                                                                        cccATCCGACCACTCCCAGTAGCCCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATg                          <  2:152829/151‑1 (MQ=255)                                                                                                                           aTACGCCCACTCCCAGTAGCGCCAACCCCCCATTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc                       <  2:693696/151‑1 (MQ=255)                                                                                                                              cgcccccacccAGTAGCGCCACCCCCCAATTTTTTCTAAACTTTCTCCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                     <  2:1144599/141‑1 (MQ=255)                                                                                                                                 ccACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                 <  1:1361141/151‑1 (MQ=255)                                                                                                                                  ccctcccAGTATCTCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCACCTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                 <  2:338272/145‑1 (MQ=255)                                                                                                                                    cTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc              >  2:594377/1‑151 (MQ=255)                                                                                                                                       ccAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc              <  1:1115185/148‑1 (MQ=255)                                                                                                                                       ccAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGt               <  1:245259/147‑1 (MQ=255)                                                                                                                                              cgcCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtat         <  1:381268/146‑1 (MQ=255)                                                                                                                                               gcCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc    <  1:495775/150‑1 (MQ=255)                                                                                                                                                ccAACACCCCACCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCg   <  1:1399531/150‑1 (MQ=255)                                                                                                                                                 cAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg  <  2:640108/150‑1 (MQ=255)                                                                                                                                                 cAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg  >  1:781040/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGG  >  NC_000913/1620970‑1621263

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