Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,717,014	0	T	G	41.6%	151.8 / 12.8	77	G156G (GGT→GGG)	yiaD	multicopy suppressor of bamB; outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (10/35);  new base G (32/0);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 3.80e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGT  >  NC_000913/3716865‑3717161                                                                                                                                                       |                                                                                                                                                    aTCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGggtg                                                                                                                                                    <  2:209582/151‑1 (MQ=255)  tCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTACGTGATTTGTTATAACGACAACACgggggg                                                                                                                                                   >  2:554392/1‑151 (MQ=255)     tCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAACGTGGTTGGGTTTACCCACAGCACggggggg                                                                                                                                                  >  1:438125/1‑148 (MQ=255)          atGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGAcc                                                                                                                                           <  2:1183559/151‑1 (MQ=255)          atGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGAcc                                                                                                                                           <  2:318609/151‑1 (MQ=255)          atGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGGGATTTGGTTTACCCGCAACACGGGGGGGCACGAcc                                                                                                                                           >  1:814501/1‑151 (MQ=255)            gCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTg                                                                                                                                         <  2:646476/151‑1 (MQ=255)             ccGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAAACGGATTTGGTATACCCCCACCACCGGGGGGGACGAccc                                                                                                                                          >  2:436977/1‑148 (MQ=255)               gAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGGGATTGGTTATACCGCCAGCCCGGGGGGGCACCACCCGaa                                                                                                                                       >  2:669933/1‑150 (MQ=255)                 aCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACAt                                                                                                                                    <  2:696796/151‑1 (MQ=255)                 aCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGGGGTCACGACCCGAACa                                                                                                                                     >  2:965044/1‑150 (MQ=255)                 aCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAACGGGGTTGGGTTTAACCGCAGCACCGGGGGGCCCGAcccgga                                                                                                                                       >  1:654105/1‑144 (MQ=255)                    aTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTTAAGGGGTTTGGTTTAACCG‑CAAGACGGGGGGGGCCGACCTGGACAAGc                                                                                                                                  >  1:389647/1‑150 (MQ=255)                    aTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTTACGAGGTTGGTTATAACGACAGCACGGGGGGGTACGACCTGAACacggg                                                                                                                                 >  1:438822/1‑147 (MQ=255)                      gtgACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGtc                                                                                                                               <  2:814507/151‑1 (MQ=255)                      gtgACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGATTGGTTATACCCACAGCACGGGGGGGCACGACCTGAACATg                                                                                                                                   >  2:536533/1‑147 (MQ=255)                        gACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGtct                                                                                                                              <  2:1104097/150‑1 (MQ=255)                        gACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGtct                                                                                                                              <  1:116049/150‑1 (MQ=255)                          ccTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGGTTGGTTATACCGACAGCACGGGGGGTCACGACCCGAACATGCGtctctc                                                                                                                           >  2:581182/1‑151 (MQ=255)                          ccTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGCTAAAGTGATTTGTTATACCCACAGCACCGGGGGTCACGACCTCGACATGGGtc                                                                                                                               >  2:554611/1‑147 (MQ=255)                              cGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAACGTGATTTGTTTTAACGACAGCACCGGGGGTCACGACCCGAACATGCGTCTCTccc                                                                                                                         >  2:874812/1‑149 (MQ=255)                                 cagcagcagCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCaa                                                                                                                    <  1:362819/151‑1 (MQ=255)                                 cagcagcagCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCaa                                                                                                                    <  2:1134590/151‑1 (MQ=255)                                 cagcagcagCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGGGATTGGGTTTACCGACAGCACGGGGGGGCACGACCCGAACATGCGTCTCTccc                                                                                                                         >  2:30940/1‑146 (MQ=255)                                 cagcagcagCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGGTAAAGGGATTGGGTTTATCCACAGGACGGGGGGGGACGACCTGAACAAGCGGCCCTccc                                                                                                                         >  1:1245629/1‑146 (MQ=255)                                  agcagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTTGTTATACCGCCAGCACCGGGGGGCACGACCTGAACATGCGTCTCTCCCCGCaa                                                                                                                    >  2:1015185/1‑150 (MQ=255)                                    cagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGt                                                                                                                 <  1:776409/151‑1 (MQ=255)                                       cagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTTTAACCACAGCACGGGGGGGCACGACCCGAACATGCGTCTCTCCCCGccaactgcg                                                                                                              >  1:273584/1‑143 (MQ=255)                                            cgACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTc                                                                                                         <  1:835854/151‑1 (MQ=255)                                              aCCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCg                                                                                                       <  2:894556/151‑1 (MQ=255)                                              aCCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGGGGTTTGGTTTAACGACAGCACGGGGGGGCACGACCTGAAAATGCGTCTCTCCCAGCAACGCGCggggt                                                                                                          >  1:808995/1‑145 (MQ=255)                                                   gAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGcc                                                                                                  <  2:452823/151‑1 (MQ=255)                                                    aaaCCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCa                                                                                                 <  1:683243/151‑1 (MQ=255)                                                      aCCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAgc                                                                                               <  2:812530/151‑1 (MQ=255)                                                       ccGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAgcg                                                                                              <  1:641644/151‑1 (MQ=255)                                                        cggcggGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAgcgc                                                                                             <  1:519916/151‑1 (MQ=255)                                                            gggcgCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGGTATACCGACAGCACGGGGGGGCACCACCTGAAAAAGCGTCTCTCCCCGCAACGGGCGGATTCCGTTGCCCGCGCGt                                                                                           >  2:602904/1‑149 (MQ=255)                                                            gggcgCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGGTTTGTTATACCGACAACACCGGGGGGCACGACCTGAACATGCGTCTCTCCCAGCAACCTGCGGGTTTCGTTGCCCGCGCGtt                                                                                          >  1:1179589/1‑150 (MQ=255)                                                             ggcgCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGGTATACCGACAACACGGGGGGTCACCACCCGAACATTCGTCTCTCCCCGCAACGGGCGGGTTTCGTTGCCAgcg                                                                                              >  1:865694/1‑145 (MQ=255)                                                              gcgcTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGGTTGGTTATAACCACAGCACGGGGGGTCACGACCTTAACATGCGTCTCTCCCCGCAACGCGCGGATTTCGTTGCCcgcgcggt                                                                                          >  1:1097298/1‑146 (MQ=255)                                                              gcgcTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGGGATTGGTTATAACGACAGCACGGGGGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCCGCGCGTTggt                                                                                       >  2:1014243/1‑149 (MQ=255)                                                               cgcTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATc                                                                                      <  2:937823/151‑1 (MQ=255)                                                                   aaCACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCAccc                                                                                  <  1:184402/151‑1 (MQ=255)                                                                   aaCACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGATTGGGTATACCGACAGCACGGGGGGGCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCCGCGCGTTGATCAccc                                                                                  >  2:279710/1‑151 (MQ=255)                                                                   aaCACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGGGATTGGGTATACCGACAACACGGGGGGGCACGACCTGAACAAGCGGCTTTCCCAGCAAAGTGCGGGTTCCGTTGCCCGCGCGTTGATCAccc                                                                                  >  2:977849/1‑151 (MQ=255)                                                                       cccTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCTCCCAgg                                                                               <  2:872263/150‑1 (MQ=255)                                                                        ccTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGc                                                                             <  1:749596/151‑1 (MQ=255)                                                                              cGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGGTATAACGACAACACGGGGGGGCACGAACCGAACAAGCGTCTCTCCCACCAAAGCGCGGGTTCCGTTGCCAGCGCGTTGATCACCCCGGGCggg                                                                          >  2:603820/1‑146 (MQ=255)                                                                               ggCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTAAGGGATTGGGTTTAACCACAGCAAGGGGGGGCACCGCCTGAACATGCGTCTCTCCCCACAACGTGCGGGATCCCTTGGCAGCGGGGTGGTAACCCAAGGCgggg                                                                         >  2:525256/1‑148 (MQ=255)                                                                               ggCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACg                                                                      <  1:581485/151‑1 (MQ=255)                                                                                     gCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTTAAGTGGTTGGTTTTACCGACAGCACGGGGGGTCACCACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCCGCCCGTTGATCACCCAAGGCGCGGAagccagcc                                                                >  2:377934/1‑151 (MQ=255)                                                                                     gCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGATTTGGTATACCGACAGCACGGGGGGTCACGACCTGAACATGCGTCTCTCCCCGCAACGTGCGGATTTCGTTGCCCGCGCGTTGATCACCCCGGGCGGGGACGCCAGcc                                                                >  1:816103/1‑151 (MQ=255)                                                                                       aaTGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGGTATACCGACAGCACGGGGGGTCACGACCTGAAAATGCGTCTCTCCCCGCAACGGGCGGGTTCCGTTGCCCGCGCGTTGATCACCCCGGGCGTGGGCGCCAGCCGc                                                              >  1:961568/1‑151 (MQ=255)                                                                                        aTGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGATTGGTTATACCCACAGCACGGGGGGGCACGCCCCGAACATGCGTCCCTCCCAGCAAAGGGCGGATTCCCTTTCCAGCCCGTTGATAACCCCAGGCGCGGGCGCCAGCCGc                                                              >  2:3326/1‑150 (MQ=255)                                                                                           gTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAACACGGGGGGGCACGACCTGAACATGCGTCTCTCCCCGCAACGTGCGGATTCCCTTGCCAGCGCGGTGATCACCCACGGCGGGGACGCCAGCCGCATcc                                                          >  1:513400/1‑151 (MQ=255)                                                                                           gTACTGAAAGAGTATCCGAAAACGGC‑GGTTAAAGTGGTTGGGTTTACCCACAGCACGGGGGGGCACGACCTGAACATGCGTCTCTCCCCGCAACGTGCGGATTCCGTTGCCCGCGCGTTGATCACCCCGGGGGGGGACGCCAGCCGccacc                                                          >  1:1243472/1‑147 (MQ=255)                                                                                                gAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACt                                                     <  2:1109598/151‑1 (MQ=255)                                                                                                 aaaGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCCAAAGCACGGGGGGTCACGACCTGAACATGCGTCTCTCCCCGCAAAGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGGGGACGCCAGCCGCCTCCGGACTc                                                    >  2:1261666/1‑151 (MQ=255)                                                                                                  aaGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCa                                                   <  2:232500/151‑1 (MQ=255)                                                                                                   agagTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAg                                                  <  1:206510/151‑1 (MQ=255)                                                                                                        aTCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGcc                                              <  2:1267955/150‑1 (MQ=255)                                                                                                                    c‑ggTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTa                                  >  1:478575/1‑150 (MQ=255)                                                                                                                          aaCGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCaa                            <  2:1039069/151‑1 (MQ=255)                                                                                                                            cGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATc                          >  2:740241/1‑151 (MQ=255)                                                                                                                             gTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCg                         >  1:721709/1‑151 (MQ=255)                                                                                                                             gTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCg                         <  2:199129/151‑1 (MQ=255)                                                                                                                               gATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGcc                       >  1:1151039/1‑151 (MQ=255)                                                                                                                                aTTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCa                      <  1:1091801/151‑1 (MQ=255)                                                                                                                                    gTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCaa                  <  2:924560/151‑1 (MQ=255)                                                                                                                                      tataCCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAAc                 <  1:539998/150‑1 (MQ=255)                                                                                                                                        taCCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGc              >  1:1091468/1‑151 (MQ=255)                                                                                                                                         aCCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGc              >  2:229357/1‑150 (MQ=255)                                                                                                                                          ccGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCa             <  1:353993/150‑1 (MQ=255)                                                                                                                                          ccGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAc            <  2:958703/151‑1 (MQ=255)                                                                                                                                            gACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCg          <  2:454387/151‑1 (MQ=255)                                                                                                                                             aCAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGc         <  2:143102/151‑1 (MQ=255)                                                                                                                                              cAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCGGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGc         >  1:634075/1‑150 (MQ=255)                                                                                                                                               aGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGc         >  1:370838/1‑149 (MQ=255)                                                                                                                                                 cACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAg       >  1:568229/1‑149 (MQ=255)                                                                                                                                                  aCGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGa      >  1:550972/1‑149 (MQ=255)                                                                                                                                                     ggtggtCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGt  >  2:1202607/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                    ATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGC‑GGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGT  >  NC_000913/3716865‑3717161

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