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

AGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993773‑994066                                                                                                                                                   |                                                                                                                                                      aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCaccacc                                                                                                                                                  >  1:979072/1‑150 (MQ=255)   aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGcc                                                                                                                                               >  2:210443/1‑151 (MQ=255)         aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAAc                                                                                                                                         >  1:330348/1‑151 (MQ=255)           ccagAAACATCATCCGTGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACACCGCGTGCGGTTGGCTTTCCCAGCCCTGCCAGCAGCCGCGGCGGGGTGTTTTT‑GCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACt                                                                                                                                        <  1:944645/147‑1 (MQ=255)           ccTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                       >  1:805932/1‑151 (MQ=255)                   atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                               >  1:729106/1‑151 (MQ=255)                     catTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                              >  1:1046779/1‑150 (MQ=255)                     catTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                              >  2:609257/1‑150 (MQ=255)                     catTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                              >  2:367708/1‑150 (MQ=255)                         cgcgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                         >  1:79782/1‑151 (MQ=255)                                 ttCCTAATTTTCAGCCGGGGGTGTGGTGCCCGCTACCACACCCCCTGCGGTTGGCTTTCCCAGCCCTCCCAGCAGGCGCAGCGGGGTTCTTTT‑CCCC‑CCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAat                                                                                                                           <  2:181024/141‑1 (MQ=255)                                       attttCAGCCGGGGGTGTGGGCCCCGCTAACACACCGCCGGGGGTGGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGTTCTTTTT‑CCCC‑CCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                           <  2:308976/149‑1 (MQ=255)                                       atttccagccaggGTTGTGGCCCCCGCAAACCCACCGCCGGGGTTTGGTTTTCCCAGCCCTCCCAGCGGGCGCAGCGGGGTCTTTTT‑CCCC‑CCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                         <  2:635161/132‑1 (MQ=255)                                          ttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                        >  1:582016/1‑151 (MQ=255)                                                cagggggggggtgCCCGCAACCACAGCCCCTGCGTTTGGCTTTTCCAGACCTGCCGGCGGGCGCAGCAGGGTTCTTTT‑CCCC‑CCCCCGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                  <  2:549945/142‑1 (MQ=255)                                                 agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                  >  2:669405/1‑150 (MQ=255)                                                  ggggTGGGGTCCCCGCAACCACTCCGCTTGGGTTTGGCTTTCCCAGACCGGCCGGCAGGCGCAGAGGGTTTTTTTT‑GCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                      <  1:758101/143‑1 (MQ=255)                                                    ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                              >  2:838345/1‑151 (MQ=255)                                                    ggTGTGGTCCCGGCTACCACTCCGCTGGGGGTGGGCTTTCCAAGACCCCCGGGAAGGCGCAGCA‑GGTTCTTTT‑TCCACCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                              <  2:881263/151‑1 (MQ=255)                                                       gtgggccccGCAAACCCCCCGCCTGGGGTTGGCTTTCCCAGCCCTGCCGGCGGGCGCAG‑AGGGTTCTTTTCCCCC‑CCCCCGCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                            <  1:396563/144‑1 (MQ=255)                                                         ggtgCCCGCTAACACACCGCCTGCGTTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTATTTTT‑CCCC‑CCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                          <  2:47495/150‑1 (MQ=255)                                                           tgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGAACTGCCAGCAGGCGCAGCAGGGTAATTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                        >  1:938179/1‑150 (MQ=255)                                                             cccccaaaCCCACCGCCGGCGGTGGGCTTTTCCAGCCCCGCCAGCAGGCGCGGCGGGGTTCTTTG‑CCCA‑CCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                                       <  2:430860/143‑1 (MQ=255)                                                             cccGCTAACACACCGCCTGCGGTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGGTTTTTTT‑GCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                      <  2:408424/150‑1 (MQ=255)                                                                gCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                                  >  2:417424/1‑151 (MQ=255)                                                                   accacaCCGCCGGCGGTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTTCTTTT‑CCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                                 <  2:361723/147‑1 (MQ=255)                                                                    acacATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                              >  2:971513/1‑151 (MQ=255)                                                                       caTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt                                                                           >  2:404192/1‑151 (MQ=255)                                                                          cGCCTGCGGTTGGTTTTCCCAGCCCGGCCGGCAGGCGCGGCGGGTTTTTTTC‑CCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                          <  2:184257/149‑1 (MQ=255)                                                                             cTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                     >  2:138379/1‑151 (MQ=255)                                                                             cTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                     >  2:666839/1‑151 (MQ=255)                                                                               gCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                   >  1:376721/1‑151 (MQ=255)                                                                                cGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                                  >  2:525557/1‑151 (MQ=255)                                                                                 ggTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                                 >  2:658596/1‑151 (MQ=255)                                                                                   ttGGCGTTCCCGGCCCGCCCCGCGGGCGCAG‑AGGGTTCTTTTGGCCC‑CCCCCGTGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                         <  1:229057/141‑1 (MQ=255)                                                                                     gtttttcccaccccTCCCCGCGGGCGCAGCGGGTTGTTTTC‑CCCC‑CCCCCGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                              <  1:721459/138‑1 (MQ=255)                                                                                     ggCTTTCCCAGCCCTGCCAGCGGGCGCGGCGGGTTTCTTTT‑GCCA‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                               <  2:93991/149‑1 (MQ=255)                                                                                     gctttccccaccccTCCCAGCGGGCCCAGCGGGTTCCTTTC‑CCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                                 <  1:627830/135‑1 (MQ=255)                                                                                      gCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                             >  1:63396/1‑150 (MQ=255)                                                                                        gTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                           >  1:499710/1‑150 (MQ=255)                                                                                         ttcccagcccccccAGCGGGCGCAGCAGGTTTCTTTT‑CCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                                 <  1:747457/131‑1 (MQ=255)                                                                                         ttcccaccccggCCGGCGGGCGCGGCAGGGTTTTTTG‑CCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                               <  2:894490/134‑1 (MQ=255)                                                                                            ccAGCCCTGCCGGCAGGCGCAGCAGGAATTTTTT‑CCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                        <  2:345999/149‑1 (MQ=255)                                                                                                 ccTCCCAGCAGGCGCAGCGGGTTACTTTT‑CCCC‑CCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                                 <  2:1035392/151‑1 (MQ=255)                                                                                                        gcagGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt                                          >  2:537304/1‑151 (MQ=255)                                                                                                           ggCGCAGCGGGTTATTTTG‑CCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                        <  2:60862/150‑1 (MQ=255)                                                                                                              gcagcaGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                    >  2:524828/1‑151 (MQ=255)                                                                                                               cagcagGGTCTTTTT‑GCCC‑CCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                    <  2:80903/150‑1 (MQ=255)                                                                                                               cagcagGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                    >  1:153624/1‑150 (MQ=255)                                                                                                                 gcggggTATTTTC‑CCCC‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                                  <  2:219746/147‑1 (MQ=255)                                                                                                                    ggtttcTTTT‑CCCA‑CCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                                <  2:1042122/144‑1 (MQ=255)                                                                                                                     ggTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                              >  1:237387/1‑150 (MQ=255)                                                                                                                          tttt‑CCCC‑CCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                            <  1:158512/147‑1 (MQ=255)                                                                                                                           ttt‑CCCA‑CCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                         <  2:294602/149‑1 (MQ=255)                                                                                                                             tcgCCC‑CCACCGCGGCGTCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa                       <  2:806055/148‑1 (MQ=255)                                                                                                                                 ccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                   <  1:643541/145‑1 (MQ=255)                                                                                                                                   ccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg                 <  1:368180/147‑1 (MQ=255)                                                                                                                                              ggCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt          <  2:303805/147‑1 (MQ=255)                                                                                                                                               gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt     >  2:24001/1‑151 (MQ=255)                                                                                                                                                ccccccccCGCCACAAACGGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc         <  1:461522/146‑1 (MQ=255)                                                                                                                                                  caccacCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAt  <  2:828997/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      AGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATGAT  >  NC_000913/993773‑994066

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