Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	42.2%	102.1 / 28.5	71	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (38/3);  new base C (0/30);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 5.69e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993764‑994063                                                                                                                                                         |                                                                                                                                                       cACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGc                                                                                                                                                       >  1:339017/1‑151 (MQ=255) cACGGTAGCAGACGCGCATCCAGACACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGAGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGc                                                                                                                                                       >  1:83029/1‑151 (MQ=255)      tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCccacc                                                                                                                                                  >  2:865056/1‑151 (MQ=255)          aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCg                                                                                                                                              >  2:314542/1‑151 (MQ=255)            aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGcc                                                                                                                                            >  2:424780/1‑151 (MQ=255)             cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCca                                                                                                                                           <  2:630291/151‑1 (MQ=255)                 cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACaa                                                                                                                                        >  2:563232/1‑150 (MQ=255)                  aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCCGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCcaca                                                                                                                                         >  1:699461/1‑148 (MQ=255)                       gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                                 >  2:312562/1‑151 (MQ=255)                              cacgcgcgTCTTTCCCTGTATTCCCGCCGGGGGGGGGGTGCCCGCAAACACCCCCCGTGCGGTGGGCTTTCCCCCCCCTCCCAGCGGGCGCAGGGGGTTCGTTTCG‑CCCCCC‑CG‑GCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGa                                                                                                                          <  1:925750/147‑1 (MQ=255)                                        tCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                                >  2:972549/1‑151 (MQ=255)                                          ttCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACACACTGACCTGCCGGAATATGTAAATcc                                                                                                              >  1:168552/1‑151 (MQ=255)                                          ttCCTAATTTTCAGCCTGGGGGGTGGCCCCCGCAACCACACCGCGGGGGGTTGCCGTTCCCAGCCCCGCCAGCGGCCGCAGCAGGGTCTTTTCC‑CCCCCA‑CG‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                     <  1:385131/144‑1 (MQ=255)                                            cctaattttCAGCCGGGGGTGTGGTCCCCGCTAACACCCCCCCTGGGGTTGGCTTTCCCAGCCCTCCCAGCGGGCGCAGCGGGGTACTTTTC‑CCCCCC‑CC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                                <  2:822153/141‑1 (MQ=255)                                            ccTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              >  1:494766/1‑149 (MQ=255)                                              tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                          >  2:422951/1‑151 (MQ=255)                                                attttCAGCCAGGGTGGTGGCCCCCGCTAACCCACCGCCGGGGTTGGGCTTTCCCAGCCCTGCGCGCGGGCGCAGC‑GGGTTTTTTTGCCCCCCA‑CC‑GCTGGGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                            <  1:987829/145‑1 (MQ=255)                                                aaTTTCATCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTACATCCAGTTgg                                                                                                        >  1:355991/1‑151 (MQ=255)                                                 aTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       >  1:607775/1‑151 (MQ=255)                                                  tttCAGCCGGGGGGGTGGTCCCCCCTACCACCCCGCCGGGGTTTGGCTTTCCCAGCCCTGCCCGCAGGCGCAGCAGGTTACTTTGG‑CCCCCC‑CC‑GCTGCGCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              <  1:838113/143‑1 (MQ=255)                                                  tttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                      >  1:54935/1‑151 (MQ=255)                                                   ttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  2:784301/1‑151 (MQ=255)                                                   ttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  1:434267/1‑151 (MQ=255)                                                   ttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  2:3279/1‑151 (MQ=255)                                                     cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  1:132534/1‑149 (MQ=255)                                                     cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  2:1073774/1‑149 (MQ=255)                                                        ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                                >  2:978266/1‑151 (MQ=255)                                                            cGGTGTGGTGCCCGATAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                             >  1:462992/1‑150 (MQ=255)                                                             ggtgggggccccGCAAACACATCCCCGGGGTTGGGCTTTTCCGGCCCGGCAAGGGGGCCCGGGAGGGTTCTTTAC‑CCCCCG‑CG‑GCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGata                                                                                            <  2:306874/141‑1 (MQ=255)                                                             ggTGGGGGGCCCCCAAACACCTCGCCGGGGTTTGGCTTTCCCAGCCCCGCCAGCAGGCGCAGCGGGGTACTTTGG‑CCCCCA‑CC‑GCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGata                                                                                            <  1:513344/150‑1 (MQ=255)                                                                gtggcccccGCAAACCCTCCGCCTGGGGTTGGCTTTCCCAGCCCTCCCAGCGGGCGCAGCGGGGTCTTTTTG‑CCCCCC‑CC‑GCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                          <  2:210933/143‑1 (MQ=255)                                                                  gggccccGCAAACACCCCGCCTGGGTTTGGCTTTCCCAGCCCTCCCGGCAGGCGCAGCAGGGTCTTTTCC‑CCCCCC‑CC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                      <  2:225254/147‑1 (MQ=255)                                                                  ggcccccGCAAACCCACCGCCGGGGGTGGGCTTTCCCAGCCCTGCCAGCGGGCGCAGAGGGTTTTTTTTG‑CCCCCC‑CC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         <  2:330869/144‑1 (MQ=255)                                                                            aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            >  2:841730/1‑151 (MQ=255)                                                                              cacaTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          >  2:58077/1‑151 (MQ=255)                                                                                         cGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATCTTTTCCGCGTAATGTTTGCTTACTGCAtt                                                               >  1:51223/1‑151 (MQ=255)                                                                                           gTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              >  2:180503/1‑150 (MQ=255)                                                                                            ttGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              >  1:86403/1‑149 (MQ=255)                                                                                              ggtttttCCAGCCCTCCCAGCAGCCCCAGCAGGTTACTTTTC‑CCCCCC‑CC‑CCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                            <  2:1008904/145‑1 (MQ=255)                                                                                                  ttcccAGCCCTGCCGGCAGGCGCAGCGGGGTTTTTTTC‑CCCCCC‑CC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                         <  2:1037896/145‑1 (MQ=255)                                                                                                  tccccggcccggCCAGCAGGCGCAGCGGGGTTTTTTTG‑CCACCC‑CC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                            <  1:487564/134‑1 (MQ=255)                                                                                                     ccAGCCCTCCCAGCAGGCGCAGCAGGGTTTTTTCC‑CCCCCC‑CC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  2:1092113/151‑1 (MQ=255)                                                                                                      cAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                  >  2:758168/1‑151 (MQ=255)                                                                                                       aGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                 >  2:392468/1‑151 (MQ=255)                                                                                                         cccgcccggcagGCCCAGCAGGGTCTTTTTG‑CCCCCC‑CC‑GCGCCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               <  2:879957/143‑1 (MQ=255)                                                                                                         cccTGCCAGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCC‑CC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               <  2:921542/150‑1 (MQ=255)                                                                                                          cccgCCAGCGGGCGCAGCGGGTTTTTTTCC‑CCCCCC‑CC‑GCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                       <  1:589532/139‑1 (MQ=255)                                                                                                             gCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                           >  2:1086550/1‑151 (MQ=255)                                                                                                               cggcagGCGCAGCAGGTTTTTTTTG‑CCCCCC‑CC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                           <  2:151471/147‑1 (MQ=255)                                                                                                               cagcagGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                         >  2:500083/1‑151 (MQ=255)                                                                                                                   agGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      >  1:741637/1‑150 (MQ=255)                                                                                                                    ggCGCAGGGGGGTTTTTTTC‑CCCCCC‑CC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                                    <  2:989527/151‑1 (MQ=255)                                                                                                                     gcgcAGAGGGGTTTTTTTC‑CCCCCC‑CC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  1:177321/135‑1 (MQ=255)                                                                                                                      cgcAGCGGGTTTTTTTTC‑CCCCCC‑CC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                     <  2:464237/148‑1 (MQ=255)                                                                                                                      ccccgcaggttttttttC‑CCCCCC‑CC‑GCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  <  1:127848/138‑1 (MQ=255)                                                                                                                       cgcgggagggtcTTTTG‑CCCCCC‑CCGGCGGCGCCCCCCCGCCGCCACAAACGGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                          <  1:644860/132‑1 (MQ=255)                                                                                                                           cagGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                              >  2:208367/1‑150 (MQ=255)                                                                                                                            gggttatttttC‑CCCCCC‑CG‑GCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg                            <  1:1029286/144‑1 (MQ=255)                                                                                                                             ggttttttttC‑CCCCCC‑CC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg                            <  2:772287/144‑1 (MQ=255)                                                                                                                             ggttctttttC‑CCCCCC‑CC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                               <  2:704019/141‑1 (MQ=255)                                                                                                                             ggttatttttG‑CCCCCC‑CC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  <  2:494846/138‑1 (MQ=255)                                                                                                                             gggtttctttccccaccc‑CC‑GCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                   <  1:319679/132‑1 (MQ=255)                                                                                                                             gggTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                           >  1:990919/1‑151 (MQ=255)                                                                                                                              gggttttttG‑CCCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                              <  2:60257/143‑1 (MQ=255)                                                                                                                              ggTTCTTTCC‑CCCCCA‑CC‑GCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                          <  2:932012/151‑1 (MQ=255)                                                                                                                               gtttctttgcccaccg‑CC‑GTTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                                <  2:48840/135‑1 (MQ=255)                                                                                                                                 ctttttG‑CCCCCC‑CG‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                               <  1:699604/141‑1 (MQ=255)                                                                                                                                 aCTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       >  1:1047537/1‑151 (MQ=255)                                                                                                                                          cacca‑cC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg               >  2:463385/1‑151 (MQ=255)                                                                                                                                                      gCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGGGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc      >  1:258299/1‑150 (MQ=255)                                                                                                                                                      gCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc      >  1:947580/1‑150 (MQ=255)                                                                                                                                                         gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  1:889214/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993764‑994063

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