Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	742,106	0	A	C	43.8%	102.2 / 19.9	73	S151A (TCT→GCT)	dtpD	dipeptide and tripeptide permease D
Reads supporting (aligned to +/- strand):  ref base A (31/10);  new base C (0/32);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 3.90e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTAACGCCGCGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAG  >  NC_000913/741958‑742250                                                                                                                                                       |                                                                                                                                                 ttGTTAACGCCGCGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGAcc                                                                                                                                                >  2:1106381/1‑151 (MQ=255)    ttAACGCCGCGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGAcccc                                                                                                                                              >  2:1097434/1‑150 (MQ=255)    ttAACGCCGCGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGAcccc                                                                                                                                              >  2:132734/1‑150 (MQ=255)     tAACGCCGCGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCAc                                                                                                                                            >  1:192055/1‑151 (MQ=255)      aaCGCCGCGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACg                                                                                                                                           >  2:4750/1‑151 (MQ=255)         gccgcGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGtt                                                                                                                                         >  1:639379/1‑150 (MQ=255)         gccgcGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTa                                                                                                                                        <  2:253866/151‑1 (MQ=255)              gggTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACccgc                                                                                                                                   >  1:557568/1‑151 (MQ=255)                        gAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATc                                                                                                                         >  2:660998/1‑151 (MQ=255)                           aTGACGATTGCCACAAAAGAAAATGCCCAGCCCCGCGATCATCCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTACCC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGc                                                                                                                      <  1:132731/151‑1 (MQ=255)                              aCGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCg                                                                                                                     >  1:364239/1‑149 (MQ=255)                                 aTTGCCACATAAGAAAATGCCCAGCCCCGCGATCATGCCACCCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTACCC‑ACAGGCGATAGGTGCGATAATAGCCCCCACGTTACCCGCCGCATACATCAGCGAGaa                                                                                                                 <  2:582261/150‑1 (MQ=255)                                   tGCCACAAAAGAAAAGGACCAGCCCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGaaa                                                                                                                <  2:174875/149‑1 (MQ=255)                                      cacaTAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGcc                                                                                                           >  1:42885/1‑151 (MQ=255)                                         aTAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCAt                                                                                                         >  1:1145805/1‑150 (MQ=255)                                           aaGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCac                                                                                                      >  2:233320/1‑151 (MQ=255)                                            aGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCacg                                                                                                     >  2:815843/1‑151 (MQ=255)                                              aaaaTGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCacgac                                                                                                   >  1:200757/1‑151 (MQ=255)                                                  tGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCgg                                                                                               >  2:446565/1‑151 (MQ=255)                                                   gACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCgg                                                                                               >  2:1051156/1‑150 (MQ=255)                                                    aCCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGAt                                                                                             >  1:119224/1‑151 (MQ=255)                                                     ccggccccGCAACCGTGCCACCCGCCGCCGGCCAAAACCCC‑AGCGCCCGTCTGTCCTCTCCCTGGGCGAAACC‑ACAGGTTAAAGGTGCAAAAATAGCCCCCCCGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATc                                                                                            <  1:908061/146‑1 (MQ=255)                                                     ccAGCCCCGCGATCATCCCACCCGCCGCCAGCCAAAAGCCC‑ATCGCCCAGCTGTACTCTTCCGGGGCGTACCC‑ACAGGCGATAGGTGCGATAATAGCCCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGa                                                                                              <  2:531044/149‑1 (MQ=255)                                                         aCCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGtt                                                                                        >  1:117829/1‑151 (MQ=255)                                                               gATCGTGCCACCCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGCCCCCCCGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCa                                                                                  <  2:182337/151‑1 (MQ=255)                                                                   aTGCCACCCGCCGCCAGGCCAAACCCC‑ATCGCCCAGCTGTACTCTTCCGGGGCGTACCC‑ACAGGCGATGGGGGCGATAATAGCCCCCCCGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTc                                                                                   <  1:244956/146‑1 (MQ=255)                                                                      ccACCCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCGGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGCCCCCACTTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACa                                                                              <  2:915504/148‑1 (MQ=255)                                                                        acccgccgCCAGCCCAACCCCC‑ACCCCCCGGCTGTCCTCCCCGTGGGCGTCCCC‑ACAGGCGATAGGTGCAAAAATAGCCCCCCCGGAACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACa                                                                              <  1:1036885/144‑1 (MQ=255)                                                                        aaCCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCg                                                                         >  2:624751/1‑151 (MQ=255)                                                                             ccgccAGGCAAAGGCCC‑ATCGCCAAGCTGTACTTTCCGGGGGCGAAACC‑ACAGGCATTAGTTGCAATAATAGACCCCCCGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGa                                                                     <  1:634777/150‑1 (MQ=255)                                                                             ccgccAGCCCAAAGCCC‑ATCGCCCGGCTGTCCTCTCCCTGGGCGAACCC‑ACGGGCGATGGGTGCGATAATAGCCCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAg                                                                    <  1:233316/151‑1 (MQ=255)                                                                             cccccgggCCAGAGCCT‑CCCGCCCACCTTTACTCTACGTGGGCGGACCC‑ACAGGCGAAAGGTGCGATAAGAGACCCCGCTTTACCCGCCGCATGCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGa                                                                     <  2:133960/144‑1 (MQ=255)                                                                              cgcccagCCAAATCCCAATCGCCCAGCTCCCCT‑TTCCCGGGCCTACCCAACGGGCGATGGGTGAGATAATAGCCCCCATCTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTc                                                                          <  1:1063136/139‑1 (MQ=255)                                                                              ccccAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCGGGGCGAAACC‑ACGGGCGATAGGTGCGATAATAGCCCCCACGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGc                                                                   <  2:458487/149‑1 (MQ=255)                                                                                ccAGGCAAAGGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACGGGCGATGGGGGCAAAAATAGACCCCACGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAg                                                                 <  2:873253/151‑1 (MQ=255)                                                                                 cgggCCAAGGCCC‑GCCGCCAGCCTGTCCTCTTCCTGGGCGTAACC‑CCTGGCGAGGGGTGCAATAACGCCCCCCACGTGACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAg                                                                 <  1:800241/148‑1 (MQ=255)                                                                                   ggCCAAAGCCC‑CGCCCCCAGCGGTCCTCCACTTGGGCGTAACC‑ACGGGCGATCGCTGGAACAAAAGCCCCCGCTTGGCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCgcgc                                                                   <  1:4750/146‑1 (MQ=255)                                                                                   ggCAAAACCCC‑ATCGCCCAGCTGTCCTCTTCCGGGGCGTACCC‑ACAGGCGATAGGGGCGATAATAGCCCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACa                                                              <  2:931702/151‑1 (MQ=255)                                                                                      cAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                                           >  1:878456/1‑151 (MQ=255)                                                                                       aagcccc‑ACCGCCCACCTGTCCTCTCCCTGGGCGTACCC‑CCAGGGAATGGGTGCGATAAGAGCCCCCCCGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCtt                                                          <  1:958951/147‑1 (MQ=255)                                                                                       aaaGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCtt                                                          >  2:546637/1‑151 (MQ=255)                                                                                       aaaCCCC‑ACCGCCCAGCTTTCCTCTCCCTGGGCGTACCC‑AAGGGCGATAGGGGGAAAAATAGCCCCCCCGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCt                                                           <  1:747413/150‑1 (MQ=255)                                                                                        aaGCCC‑ATCGCCCAGCTGTACTCTCCCTGGGCGTAACC‑ACAGGCGATAGGGGCGAAAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTa                                                         <  2:165645/151‑1 (MQ=255)                                                                                            cc‑ATCGCCCAGCTTTCCTCTTCCTGGGCATACCC‑ACGGACGATGGGTGCAAAAATGGCCCCCACTTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt                                                     <  1:463797/151‑1 (MQ=255)                                                                                            cc‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCACCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt                                                     <  2:79409/151‑1 (MQ=255)                                                                                               accGCCCAGCTGTCCTCTTCCGGGGCGTACCC‑CCGGGCAATAGGTGCAATAATAGCCCCCCCGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGtt                                                     <  2:216926/147‑1 (MQ=255)                                                                                                   cccAGTTGTACTCTTCCTGGGCGTACCC‑ACAGGCAATAGGGGCGATAATAGCCCCCCCGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGAttt                                               <  2:746193/151‑1 (MQ=255)                                                                                                      aGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTaa                                             >  2:1018410/1‑150 (MQ=255)                                                                                                           tACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGc                                       >  2:171285/1‑151 (MQ=255)                                                                                                               ctcccTGGGCGTACCC‑ACAGGCGATGGGTGCAATAATAGCCCCCCCTTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGc                                       <  2:18725/144‑1 (MQ=255)                                                                                                               ctTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCAt                                    >  1:25930/1‑150 (MQ=255)                                                                                                                ttCCTGGGCGAAACC‑CAAGGCGATAGGTGCGAAAAAAGCCCCCCCTTTACCCGCCGCATCCATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGc                                       <  2:150697/146‑1 (MQ=255)                                                                                                                tcccggggCGTAACC‑ACAGGCGATAGGTGCGATAATAGCCCCCCCTTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACACGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACa                                          <  1:467844/138‑1 (MQ=255)                                                                                                                 tCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAg                                  >  1:650361/1‑150 (MQ=255)                                                                                                                 tCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGc                                 >  1:641145/1‑151 (MQ=255)                                                                                                                  ccgggggttAACC‑ACGGGGGATAGGGGCAAAAATGGCCCCCCCGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGcc                                      <  2:558901/137‑1 (MQ=255)                                                                                                                  ccgggggGTACCC‑ACGGGCGATGGGTGCGATAATAGCCCCCGCGTTCCCCGCCCCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCa                                     <  1:771770/143‑1 (MQ=255)                                                                                                                  ccTGGGCGAACCC‑AAAGGCGATGGGTGCGAAAAAAGCCCCCCCTTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAg                                  <  2:13150/149‑1 (MQ=255)                                                                                                                   cTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCg                               >  1:267027/1‑151 (MQ=255)                                                                                                                    tGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGc                              >  1:598342/1‑151 (MQ=255)                                                                                                                           aaCC‑ACAGGGGATAGGTGCGATAATAGCCCCCAGGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGAt                       <  1:357152/151‑1 (MQ=255)                                                                                                                            acccccgggCGATGGGTGCAAAAATAGCCCCCACGTCCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACg                         <  1:719856/142‑1 (MQ=255)                                                                                                                            acccaatgggaagaGGTACAATAATAGCCCCCCCGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGa                        <  1:1013361/137‑1 (MQ=255)                                                                                                                             cc‑ACAGGCGATAGGTGCGATAAAAAACCCCCCGTACCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACg                         <  2:6769/147‑1 (MQ=255)                                                                                                                                      gATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAgg             >  2:338132/1‑151 (MQ=255)                                                                                                                                       aTAGGGGCGATAATAGCCCCCCCGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAgg             <  2:641166/150‑1 (MQ=255)                                                                                                                                         gggggcaataataGCCCCCCCTTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACCCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGAc          <  1:605863/144‑1 (MQ=255)                                                                                                                                         aGGTGCGATAATAGACCCCACGTTACCCGCTGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGAc          >  1:589669/1‑151 (MQ=255)                                                                                                                                          ggggcaaaaataGCCCCCCCTTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGa           <  1:1051094/141‑1 (MQ=255)                                                                                                                                          ggggCGAAAATAGACCCCACGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCa               <  1:329955/142‑1 (MQ=255)                                                                                                                                             gCGATAATAGCCCCCACGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGAt      <  2:1155250/151‑1 (MQ=255)                                                                                                                                                 taataGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATag    >  2:677566/1‑149 (MQ=255)                                                                                                                                                  aataGCCCCCACGTTCCCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATagag  <  1:239344/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 TTGTTAACGCCGCGGGTATGAGTGAAATGACGATTGCCACATAAGAAAATGACCAGACCCGCGATCATGCCAACCGCCGCCAGGCCAAAGCCC‑ATCGCCCAGCTGTACTCTTCCTGGGCGTAACC‑ACAGGCGATAGGTGCGATAATAGACCCCACGTTACCCGCCGCATACATCAGCGAGAAACCGCCATCACGACGCGGATCGGTTGGCTCATACAGCTCGCCGAGCAGACAGCTTACGTTAGATTTAAACAGGCCATAGCCGCAGACGATAATCGCCAGGGACAGATAGAG  >  NC_000913/741958‑742250

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