Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	46.6%	90.0 / 14.9	60	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (27/4);  new base C (0/27);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.23e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993769‑994066                                                                                                                                                       |                                                                                                                                                      tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                     >  2:645656/1‑151 (MQ=255) tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                     >  1:521752/1‑151 (MQ=255)  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacca                                                                                                                                                    <  1:190989/151‑1 (MQ=255)  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                     >  1:721488/1‑150 (MQ=255)                 tGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                     >  2:373483/1‑151 (MQ=255)                       atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                               >  1:164485/1‑151 (MQ=255)                           ttCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                           >  2:253749/1‑151 (MQ=255)                            tCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                           >  1:284102/1‑150 (MQ=255)                             cgcgTGTCTTCCTAAATTCCAGCCGGGGGTGTGGTGCCCGCTAACACACCGCCGGCGTTTGGCTTTCCCAGCCCTG‑CCAGCGGGCGCAGCAGGTTATTTTTC‑CCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                         <  1:331404/151‑1 (MQ=255)                                       ccgaaaTTCCACCCGGGGGTGTGGTCCCCGCAAACACACCGCCGGGGGTGGGCTTTCCCAGCCCTG‑CCAGCAGGCGCAGCGGGGTATTTTTC‑CCACCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                               <  1:583199/147‑1 (MQ=255)                                       ccTGATTTTCAGCCGGGGGTGTGGCCCCCGCAAACCCTCCCCCTGGGGTTGGCTTTCCCAGCCCGC‑CCGGCGGCCGCAGCGGGTTCTTTTTC‑CCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                                  <  2:277611/148‑1 (MQ=255)                                       ccTGATTTCCAGCCAGGGGTGGGGTGCCCGCTAACACACCGCCTGGGGTTGGCGTTTCCAGCCCGG‑CCAGCAGGCGCAGCAGGGTTCTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                               <  1:688071/151‑1 (MQ=255)                                       ccTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                                 >  1:890961/1‑149 (MQ=255)                                       ccTAAATTTCACCCGGGGGGGGGGTCCCCTCTACCACCCCGCTGGGGGTGGGTTTTCCCAGCCCCC‑CCAGCAGCCGCAGCGGGTTCCTTCTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                                 <  2:909437/149‑1 (MQ=255)                                           atttccagccaggGGTGGGGGCCCCGCAAACACACCGCCTGGGTTTGGCTTTCCCAGCCCTG‑CCAGCGGGCGCAGCAGGTTCCTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                            <  1:939509/145‑1 (MQ=255)                                                 agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                     >  2:493265/1‑151 (MQ=255)                                                    cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                  >  2:208990/1‑151 (MQ=255)                                                     agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                  >  2:535478/1‑150 (MQ=255)                                                        ggtgggggccccGCAAACCCCCCGCCGGCGGTGGGCTTTCCCAGCCCGC‑CCGGCGGGCGCAGCGGGTTATTTTTC‑CCCCCACCGCTGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCa                                                                                                       <  2:718895/133‑1 (MQ=255)                                                             gggccccGCAAACCCGTCGCCTGGGGTTGGCTTTCCCAGCCCTG‑CCAGCGGGCGCAGC‑GGGTTTTTTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                         <  1:650050/147‑1 (MQ=255)                                                                  ccGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                     >  2:155847/1‑150 (MQ=255)                                                                       accccgccGCCTGGGGTGGGCTTTCCCAGCCCGC‑CCAGCGGGCGCAGCAGGTTTTTTTGC‑CCACCACGGCTGGGCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                               <  2:181108/144‑1 (MQ=255)                                                                       aaCACGTCGCCTGGGGTTGGCTTTTCCAGCCCTG‑CCAGCAGGCGCAGCGGGGTCTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                                <  1:101457/150‑1 (MQ=255)                                                                       aaCACACCGCCTGCGGTTGGC‑TTTCCCGCCCCGCCCAGCAGGCGCAGCAGGTTCTTTTTC‑CCCCCACCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                   <  2:811858/147‑1 (MQ=255)                                                                        acacATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                              >  1:701729/1‑151 (MQ=255)                                                                         cacaTCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                             >  1:69985/1‑151 (MQ=255)                                                                          acaTCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                            >  1:600680/1‑151 (MQ=255)                                                                           caccgccggggtttGGCGTTCCCGACCCGG‑CCAGCAGGCGCAGCAGGTTCCTTTTC‑CCCCCACCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                            <  2:382576/138‑1 (MQ=255)                                                                                ccTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                      >  1:658905/1‑151 (MQ=255)                                                                                  tGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                     >  2:628651/1‑150 (MQ=255)                                                                                   gCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                   >  2:227464/1‑151 (MQ=255)                                                                                   gCGGTTGGCGTTCCCAGCCCTG‑CCGGCAGGCGCGGCGGGGTCTTTTTG‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                                    <  1:123026/150‑1 (MQ=255)                                                                                     gtttGGCTTTCCCAGCCCTG‑CCAGCGGGCGCAGCAGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                      <  1:898050/144‑1 (MQ=255)                                                                                     ggTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                                 >  1:564406/1‑151 (MQ=255)                                                                                        tGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                              >  1:693805/1‑151 (MQ=255)                                                                                              ttcccagcccgcccAGCGGGCGCAGAAGGGTCTTTTTC‑CCCCCCCGGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTAAGc                                                              <  2:113208/134‑1 (MQ=255)                                                                                              tccaagcccTG‑CCGGCAGGCGCAGCGGGTTTTTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                         <  1:681120/143‑1 (MQ=255)                                                                                                ccAGCCCTC‑CCAGCGGCCGCAGCAGGGTGTTTTCG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                       <  1:836839/150‑1 (MQ=255)                                                                                                ccAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                      >  2:118099/1‑151 (MQ=255)                                                                                                ccAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                      >  1:141655/1‑151 (MQ=255)                                                                                                 ccaccccgccgagcagGCGCAGGAGGTTTCTTTCG‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                          <  2:729449/136‑1 (MQ=255)                                                                                                  aGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                      >  2:635870/1‑149 (MQ=255)                                                                                                      ccgccaagaagGCCCAGCGGGTTATTTTGC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                     <  1:263883/138‑1 (MQ=255)                                                                                                          ccagagggggggggggggTCCTTTTT‑CAACCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                             <  1:340084/136‑1 (MQ=255)                                                                                                              cagGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                         >  2:179719/1‑151 (MQ=255)                                                                                                               agGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCTGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                         >  1:805971/1‑150 (MQ=255)                                                                                                               agGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                         >  1:327333/1‑150 (MQ=255)                                                                                                                ggcgacgcgggttctttttG‑CCCCCCCGGTTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                                       <  1:41826/136‑1 (MQ=255)                                                                                                                ggcccaacgggctACTTTCC‑CCACCCCCGCGGCGCCCCTCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                        <  1:247560/138‑1 (MQ=255)                                                                                                                ggCGCAGCGGGTTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt                                          <  2:190993/148‑1 (MQ=255)                                                                                                                ggCGCAGCGGGTTACTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                        <  1:682247/150‑1 (MQ=255)                                                                                                                    cagcggGTTATTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                      <  1:128135/148‑1 (MQ=255)                                                                                                                     agcagGGTTATTTTC‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                    <  2:858121/149‑1 (MQ=255)                                                                                                                        ggggtttttttG‑CCCCCACCCTTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                                <  1:253766/143‑1 (MQ=255)                                                                                                                        agGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                                 >  1:28311/1‑149 (MQ=255)                                                                                                                          gttACTTTTG‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                                   <  2:863671/143‑1 (MQ=255)                                                                                                                                     ccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                   <  1:645728/145‑1 (MQ=255)                                                                                                                                      cccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg                  <  1:1369/146‑1 (MQ=255)                                                                                                                                      cccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg                  <  2:860181/146‑1 (MQ=255)                                                                                                                                        cccccGCTGCGGCCCGCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGATCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg                  <  1:942863/146‑1 (MQ=255)                                                                                                                                         caccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc               >  2:315023/1‑151 (MQ=255)                                                                                                                                          accGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct              >  1:663407/1‑151 (MQ=255)                                                                                                                                           ccGGTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc             <  1:262954/151‑1 (MQ=255)                                                                                                                                             gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct           >  1:35035/1‑151 (MQ=255)                                                                                                                                                    ccccccGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt          <  2:32771/145‑1 (MQ=255)                                                                                                                                                       accaccGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAt  >  2:893098/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      TAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993769‑994066

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