Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	47.5%	95.7 / 22.3	60	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (29/2);  new base C (0/28);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 8.39e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993766‑994063                                                                                                                                                      |                                                                                                                                                       cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGccc                                                                                                                                                     <  1:216836/151‑1 (MQ=255)  ggTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGccc                                                                                                                                                     >  1:195989/1‑150 (MQ=255)  ggTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGccc                                                                                                                                                     >  1:605672/1‑150 (MQ=255)        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)                                                                                                                                                      |                                                                                                                                                       CGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCA‑CCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993766‑994063

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