Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	827,507	0	T	G	45.1%	116.1 / 11.9	71	D492A (GAC→GCC)	ybhF	putative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (9/30);  new base G (32/0);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 2.49e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGC  >  NC_000913/827357‑827654                                                                                                                                                       |                                                                                                                                                    gTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGt                                                                                                                                                     <  2:578455/151‑1 (MQ=255)     cGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCaa                                                                                                                                                  >  1:351295/1‑150 (MQ=255)      gCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCaa                                                                                                                                                  >  1:752498/1‑149 (MQ=255)         ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCg                                                                                                                                             <  1:408886/151‑1 (MQ=255)           cGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCaaaggccg                                                                                                                                            >  2:1216550/1‑144 (MQ=255)               tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGt                                                                                                                                        <  1:925735/150‑1 (MQ=255)               tAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCggga                                                                                                                                          >  2:526457/1‑146 (MQ=255)                aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTc                                                                                                                                       <  1:41311/150‑1 (MQ=255)                     cccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcg                                                                                                                                 <  1:500729/151‑1 (MQ=255)                       cgcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtc                                                                                                                               <  2:211853/151‑1 (MQ=255)                        gcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGTGAACGCCGGAAGTCGGTtcgtcg                                                                                                                              >  1:264977/1‑151 (MQ=255)                        gcgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCCGAAGCCgggttgt                                                                                                                                >  2:729774/1‑144 (MQ=255)                          ggTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGTTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCgag                                                                                                                            <  2:577430/151‑1 (MQ=255)                              cacCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAGTCGGTTCGTCGCGaaa                                                                                                                         >  1:1154133/1‑150 (MQ=255)                               acCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACa                                                                                                                       <  2:530743/151‑1 (MQ=255)                                 cAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGa                                                                                                                     <  1:122256/151‑1 (MQ=255)                                        aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTTGTCGAGAAAAAGAATGTcc                                                                                                               >  1:296177/1‑150 (MQ=255)                                        aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGCGCGGGGGGGAACGCCGGAAAGCGGGTACTCGGGAAACAGAAAGTcc                                                                                                               >  2:666419/1‑150 (MQ=255)                                          gCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCgg                                                                                                             <  2:569327/150‑1 (MQ=255)                                           cGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt                                                                                                           <  2:541086/151‑1 (MQ=255)                                           cGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCCACGCCGGAAGACGGTTTGTCGAGAAACAGAATGTACgg                                                                                                             >  2:1168915/1‑149 (MQ=255)                                            ggTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTc                                                                                                          <  1:239779/151‑1 (MQ=255)                                            ggTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAATGGCCGGtt                                                                                                           >  2:401252/1‑150 (MQ=255)                                             gTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCCGGAGTCGGGTCGTCGAGAAACAGAATGTCCgggt                                                                                                           >  2:879839/1‑147 (MQ=255)                                               cGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAACAGAATGTCCCGTTTAt                                                                                                        >  1:1173539/1‑150 (MQ=255)                                                 cAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCa                                                                                                     <  1:1246261/151‑1 (MQ=255)                                                       tCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGc                                                                                                <  2:1118320/150‑1 (MQ=255)                                                        ccGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGa                                                                                              <  2:457192/151‑1 (MQ=255)                                                         cGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGTCCACGGCGGCAGATGGTCTCTCCAGCAAGATAATGGTCGGTTCCTGCAtgttcg                                                                                               >  2:952668/1‑144 (MQ=255)                                                              cATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGGCGGGTCGTCGGGAAAAAGAATGTCCGGGTCATGCATCAGCGAAAAgg                                                                                         >  1:363878/1‑150 (MQ=255)                                                              cATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGAGGGGGGGGTAACCCCCGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCCGCGa                                                                                              >  2:1015015/1‑145 (MQ=255)                                                                tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTGCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAg                                                                                          <  2:396169/147‑1 (MQ=255)                                                                  cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagc                                                                                    >  2:1041973/1‑151 (MQ=255)                                                                     aaaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCTGCGAACAGgccagcggc                                                                                 >  2:1094102/1‑149 (MQ=255)                                                                        gtgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGAAGGCGGGTCGTCCAGAAAAAGAAAGTTCCGGTCATGCATCAACGGACAGgcccgcgcc                                                                                 >  1:1097856/1‑148 (MQ=255)                                                                         tgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGCATTAGCGAACAGGCCAGCGCCAGcc                                                                             >  2:840166/1‑151 (MQ=255)                                                                                 gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGttt                                                                      <  2:514587/150‑1 (MQ=255)                                                                                           gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATg                                                           <  1:460024/151‑1 (MQ=255)                                                                                              aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCa                                                        <  1:879751/151‑1 (MQ=255)                                                                                              aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCa                                                        <  2:486674/151‑1 (MQ=255)                                                                                                gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCa                                                        <  2:94718/149‑1 (MQ=255)                                                                                                 ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTTGTCGAGAAACAGAATGTTCGGGTTATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAACTAATGGCAGtt                                                     >  2:110685/1‑151 (MQ=255)                                                                                                  cTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCCGCGAACAGGGCCGCGCCAGCCGCTGTTTTAAAACTAAAGGCAggt                                                     >  1:201575/1‑148 (MQ=255)                                                                                                   tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCCGTTTTAAAACTAAAGGCAGtt                                                     >  1:478900/1‑149 (MQ=255)                                                                                                    ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTc                                                    <  1:113652/149‑1 (MQ=255)                                                                                                     tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCAt                                                  <  2:988868/150‑1 (MQ=255)                                                                                                      ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAg                                                       >  2:1177077/1‑144 (MQ=255)                                                                                                        ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCGGTTTTAAAACTaaatggcgtttattgg                                               >  1:973977/1‑136 (MQ=255)                                                                                                        ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGCAGACGGTTTGTTCAGAAACAGAAAGTCCGGGTCATGCCTCAACGAACAGGCCCGCGGCAGCCGCTGGTTAAAACCTaaaaggagattttcgg                                               >  1:730205/1‑136 (MQ=255)                                                                                                        ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGACGGTTCGTCGAGAAACAGAAAGTCCGGGTCATTCATCAGCGAACAAGGCCGCGCCAGCCGCTGTTTAAAACCTAATGGCAGGTCAt                                                  >  1:1099269/1‑147 (MQ=255)                                                                                                         tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGTGGGTTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGGTTTAAAACTAATTGCAGTTCATctggg                                             >  2:124098/1‑147 (MQ=255)                                                                                                          aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGt                                              <  2:1089107/149‑1 (MQ=255)                                                                                                          aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCCACGCCGGAAGGCGGGTCGTCGAGAAACAGAAAGTCCGGTTCATGCATCAGCGGAAAAGGCAGCGCCCGCCGCTGTTTTAAACCTAAAGGCAGGTCAt                                                  >  1:1185530/1‑145 (MQ=255)                                                                                                           ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAAAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCCGTTTAAAAACTAATGGCAGTTCCTCGGTgg                                            >  1:984900/1‑150 (MQ=255)                                                                                                             aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAAAAGGCCAGCGCCCGCCGCCGGTTTAAACCTAATGGGAGGTCATCGGGGGCg                                          >  1:478555/1‑150 (MQ=255)                                                                                                               gCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTgg                                       <  1:221862/151‑1 (MQ=255)                                                                                                                cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCCGAAGTCGGGTTGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAACAGGGCAGCGCCAGCCGCCGTTTTAAAACTAATGGCAGGTCATCGGGGGCGTggg                                      >  2:684410/1‑151 (MQ=255)                                                                                                                 tGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGTGGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCCTCAGCGAACAGGGCCGCGCCAGCCGCTGTTTAAAAACTAAAGGCAGTTTATCGGGGGCGTggg                                      >  2:604725/1‑150 (MQ=255)                                                                                                                   ttGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCCGTTTAAAAACTAAAGGCAGTTCATCGGGGGCGTgggggg                                   >  2:690500/1‑151 (MQ=255)                                                                                                                    tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGGTCGTCGAGAAAAAAAAGGTCCGGGGCATGCATCAGCGAAAAGGCCAGCGCCCGCCCCTGTTTTAAAACTAAAGGCAGTTCATCgggggggggggggg                                   >  1:46808/1‑150 (MQ=255)                                                                                                                         tgCAGCCAAAATTCACGGCGGGGGAGGGGGTGAACGCCGGAAGGCGGGTCGTCGAGAAACAGAATGTCCGGGTTCTGCCTCAGCGAACAGGCCCGCGCCAGCCGCTGTTTTAAACCTAATGGGAGTTGTTCCGGGGGGGGGGGGGCGATAc                             >  1:1190491/1‑151 (MQ=255)                                                                                                                                 aaaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAgg                     <  2:371274/151‑1 (MQ=255)                                                                                                                                     tCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAATAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCg                  <  1:526409/150‑1 (MQ=255)                                                                                                                                         ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGc             >  1:1202683/1‑151 (MQ=255)                                                                                                                                         ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGc             <  2:351330/151‑1 (MQ=255)                                                                                                                                           cggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCt           <  2:1173666/151‑1 (MQ=255)                                                                                                                                            gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTc          >  1:925156/1‑151 (MQ=255)                                                                                                                                               tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCt       >  2:270424/1‑151 (MQ=255)                                                                                                                                                 aGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTc      >  1:911944/1‑150 (MQ=255)                                                                                                                                                   ggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg   <  2:752563/151‑1 (MQ=255)                                                                                                                                                    gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGc  >  2:60616/1‑151 (MQ=255)                                                                                                                                                    gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGc  <  2:177236/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    GTGCCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATGC  >  NC_000913/827357‑827654

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