Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	385,797	0	T	G	41.2%	122.2 / 16.5	69	V189G (GTG→GGG)	tauA	taurine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (9/31);  new base G (28/0);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 2.39e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTA  >  NC_000913/385647‑385943                                                                                                                                                       |                                                                                                                                                   cggcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgt                                                                                                                                                    <  1:858780/151‑1 (MQ=255)   gcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtgt                                                                                                                                                  <  1:412429/151‑1 (MQ=255)     aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCCAGGGGtt                                                                                                                                                 >  2:470805/1‑150 (MQ=255)        gAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAcc                                                                                                                                             <  1:243949/151‑1 (MQ=255)         aaaCACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGGCAGGGGGTGAcc                                                                                                                                             >  1:476058/1‑150 (MQ=255)           acacTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGAt                                                                                                                                          <  1:359834/151‑1 (MQ=255)               tGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAAAAGGCAAGGGGTTGAGCGGTTCt                                                                                                                                       >  2:752134/1‑150 (MQ=255)                ggggCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCCAGGTGTTGACCGATTCt                                                                                                                                       >  2:125869/1‑149 (MQ=255)                        aaaCCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGGGTTGACCGATTCTGAACAgg                                                                                                                                >  1:790371/1‑148 (MQ=255)                        aaaCCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCTGGAAAAAGACGGCAAGGGGGTGACCGATTCTGAACAGGt                                                                                                                               >  2:294835/1‑149 (MQ=255)                            ccGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCa                                                                                                                         <  1:657314/151‑1 (MQ=255)                               ggCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTg                                                                                                                      <  1:702878/151‑1 (MQ=255)                                gCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGAAAAAAAAGGGAAGGGGTTTACCGATTCTGAACAAGGCGGGGCGTgg                                                                                                                     >  2:399501/1‑151 (MQ=255)                                 cAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCTGGAAAAAAAAGGGCAGGGGTTGACCGATTCTGAACAGGGCGGGCAGgggg                                                                                                                    >  2:666047/1‑151 (MQ=255)                                  aaGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGAAGGCAAGGGGTTGACCCGTTCTGAACAGGGCGGGCAGTgggg                                                                                                                   >  2:183433/1‑151 (MQ=255)                                   aGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCGTGGGAAACACGGGAAGGGGTTGAACGAGTTTTAACAAGGCggggcgggggg                                                                                                                   >  1:105611/1‑142 (MQ=255)                                      ggAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGGCAGGGGTTGGCCGATTCTGAACAGGTCggggcgggaggcg                                                                                                                 >  1:1047838/1‑139 (MQ=255)                                      ggAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAAACGGCAAGGGGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcg                                                                                                                 >  2:664798/1‑149 (MQ=255)                                          aTTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCaa                                                                                                           <  2:903402/151‑1 (MQ=255)                                            tgtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAAACGGGAAGGTGGTGACCGATTTTGAACAAGTCGGGCCGGGGGGCGCGCCa                                                                                                            >  2:810738/1‑148 (MQ=255)                                             gtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGGAAAAGAAGGCAAGGGGTTGACCGATTCTGAACAGGGCGGGCCGGGGGGGGCGccca                                                                                                           >  1:536431/1‑146 (MQ=255)                                               gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGCAAGGGGGTGACCGATTCTGAACAGGGCGGGCGGGGGGgggc                                                                                                                >  1:49391/1‑141 (MQ=255)                                               gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCt                                                                                                       <  1:925171/150‑1 (MQ=255)                                                aaCCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTgg                                                                                                     <  2:629897/151‑1 (MQ=255)                                                    tGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCTGGGAAAAAACGGGAAGGGGGTGGCCGATTTTGGAAAGGGCGGGGCGGGGGGCGCGGCAACGCTgg                                                                                                     >  1:1016535/1‑147 (MQ=255)                                                      cAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCt                                                                                               <  1:342057/151‑1 (MQ=255)                                                      cAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGGGGTGACCGATTCTGAACAGGTCGGGCAGGGGGGCGCGCCAACGCTGGACg                                                                                                  >  2:571917/1‑148 (MQ=255)                                                       aGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTTACGCCCTGGGAAAAGACGGGAAGGGGTTTTCCGGTTCTGAAAAGGGCGGGGAGTGGGGGGCGCCAACGCTGGGCGTCTg                                                                                              >  1:213871/1‑151 (MQ=255)                                                       aGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTTACGCCCTGGGAAAAGACGGCAAGGGGGTGACCGGTTCTGAACAGGGCCGGCAGGGGGGGGCGCCCACGCTGGACGTCt                                                                                               >  1:255860/1‑150 (MQ=255)                                                        gccgccCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTgg                                                                                             <  2:623854/151‑1 (MQ=255)                                                         ccgcccgcGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGGGTTGACCGATTCTGAACAGGGCGGGCAGGGGGGCGCGCCAACGCTGGGCGGCTggg                                                                                            >  1:999605/1‑151 (MQ=255)                                                              cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGGAAGGGGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGGCAAAGCTGGACGTCTGGgtggtg                                                                                       >  1:737766/1‑151 (MQ=255)                                                              cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGACCAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGgtggtg                                                                                       <  1:485227/151‑1 (MQ=255)                                                                  aTTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCa                                                                                   <  1:60636/151‑1 (MQ=255)                                                                         cTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTTACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGGCGGGCAGTGGGGCGCGCCAACGCTGGACGGCTGGGGGGTGGGCAAAAAttt                                                                            >  1:177203/1‑151 (MQ=255)                                                                                gCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGGGGTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGGGGTGCGCAAAGATTTTGCCg                                                                       >  2:855046/1‑149 (MQ=255)                                                                                 cAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGGTAACGGCCTGGAAAAAGACGGCAAGGAGTTGAACGATTCTTGACAAGGCGGGGCGGGGGGCGCGGCAACGCTGGACGGCTGGGGGGGGCGCAAAGATTTTGCCggg                                                                     >  2:180635/1‑148 (MQ=255)                                                                                   gCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGaaa                                                                  <  2:363356/151‑1 (MQ=255)                                                                                    cGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAAc                                                                 <  1:226469/151‑1 (MQ=255)                                                                                     ggggAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACa                                                                <  1:548577/151‑1 (MQ=255)                                                                                     ggggAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCTGGGAAAAGACGGCAAGGTGTTGACCGGTTCTGAAAAGGTCGGGCCGGGGGGGGCGCCAACGCTGGGCGTCTGGGTGGGGCGCAAAGATTTTGCCGAGaaa                                                                  >  1:887280/1‑149 (MQ=255)                                                                                      gggAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACAt                                                               <  1:780995/151‑1 (MQ=255)                                                                                       ggAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTTAAGCCCTGGGAAAAGACGGCAAGGGGTTGACCGGTTCTGAACAGGGCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACAt                                                               >  1:236354/1‑150 (MQ=255)                                                                                       ggAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATc                                                              <  2:137121/151‑1 (MQ=255)                                                                                        gagaTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGAAAAAGACGGCAAGGGGTTGACCGATTCTGAACAGGGCGGGCAGTGGGGCGCGCCAACGCCGGACGTCTGGGTGGGGCGCAAAGATTTTTCCGAGAAAAATcc                                                             >  2:145010/1‑151 (MQ=255)                                                                                          gaTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCt                                                            <  2:598555/150‑1 (MQ=255)                                                                                               tGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGAAAAAAACGGCAAGGGGTTGACCGATTCTGAACAGGGCGGGCAGGGGGGCGCGCCCACGCCGGACGTCTGGGGGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAgg                                                        >  2:1023388/1‑149 (MQ=255)                                                                                                  tGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAAACGGCAAGGGGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGGGGTGCGCCAAGATTTTTCCGAGAAACATCCTGAGGTCGTg                                                   >  2:99075/1‑151 (MQ=255)                                                                                                   ggTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTg                                                   <  1:706913/150‑1 (MQ=255)                                                                                                     tGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAAACGGCAAGGGGGTGACCGATTCTGAACAGGGCGGGCAGGGGGGCGCGCCAACGCTGGACGTTTGGGGGGGGCGCAAAGATTTTTCCGAGAAACATCCTGGGGTCGTGaa                                                 >  2:38451/1‑150 (MQ=255)                                                                                                      gCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAg                                               <  2:702767/151‑1 (MQ=255)                                                                                                      gCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAg                                               <  2:832405/151‑1 (MQ=255)                                                                                                      gCTTATGTCTGGGCACCGGCGGCTAACGCCCCGGAAAAAAACGGCAAGGGGTTGACCGATTCTGAACAGGTCGGGGCGTGGGGCGCGCCAACGCTGGGCGTCTGGGGGGGGCGCAAAGATTTTTCCGCGAAACATCCTGAGGGCGTgga                                                 >  1:491311/1‑147 (MQ=255)                                                                                                        ttATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGCAAGGGGTTGACCGGTTCTGAACAGGGCGGGCAGTGGGGCGCGCCAACGCTGGGCGTCTGGGGGGGGGGCAAAGATTTTTCCGGGAAACAACCTGGGGTCGTGGAAg                                               >  1:215750/1‑149 (MQ=255)                                                                                                         tATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGc                                              <  2:979392/149‑1 (MQ=255)                                                                                                          aTGTCTGGGCCCCGGCGGTTAACTCCCTGTAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCAAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGt                                            <  1:1044850/150‑1 (MQ=255)                                                                                                               tGGGCACCGGCGGTTAACGCCCCGGAAAAAGACGGCCAGGGGTTGGCCGAGTCTGAACAGGGCGGGCAGTGGGGCGCGCCCAAGCGGGACGTCTGGGGGGTGCGGAAAGATTTTGCCGAGCAACATCCCGAGGTCGTGGAAGCGGTCGc                                        >  1:463846/1‑149 (MQ=255)                                                                                                                          ggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATc                           <  2:641084/151‑1 (MQ=255)                                                                                                                          ggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATc                           <  1:874159/151‑1 (MQ=255)                                                                                                                           gTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGAGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCg                          >  2:1022636/1‑151 (MQ=255)                                                                                                                                 gCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTc                    >  2:310335/1‑151 (MQ=255)                                                                                                                                  cccTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTc                    >  1:803529/1‑150 (MQ=255)                                                                                                                                  cccTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCa                   <  1:551942/151‑1 (MQ=255)                                                                                                                                      ggAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCaa               <  2:733677/151‑1 (MQ=255)                                                                                                                                      ggAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCa                <  2:653695/150‑1 (MQ=255)                                                                                                                                      ggAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCa                <  2:286907/150‑1 (MQ=255)                                                                                                                                         aaaaGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCg            >  1:779695/1‑151 (MQ=255)                                                                                                                                               cGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACAtt      <  2:920984/151‑1 (MQ=255)                                                                                                                                               cGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACAtt      >  1:28472/1‑151 (MQ=255)                                                                                                                                                ggCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTg     <  1:1048609/151‑1 (MQ=255)                                                                                                                                                   aaGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTa  <  1:358118/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTA  >  NC_000913/385647‑385943

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