Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	52.2%	66.2 / 114.4	67	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (22/10);  new base C (0/35);  total (22/45)
Fisher's exact test for biased strand distribution p-value = 2.38e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-10

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993768‑994063                                                                                                                                                       |                                                                                                                                                       gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCcca                                                                                                                                                    >  1:826553/1‑150 (MQ=255) gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCcca                                                                                                                                                    >  1:698676/1‑150 (MQ=255) gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTT‑G‑CCACCA‑CCGCTGCGGCccac                                                                                                                                                   >  1:204577/1‑150 (MQ=255)  tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCccacc                                                                                                                                                  <  2:1277748/151‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                <  2:1142422/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                <  2:1066793/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                <  1:324505/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                <  2:93614/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                <  1:1373058/150‑1 (MQ=255)               tCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACaaa                                                                                                                                       >  1:957140/1‑149 (MQ=255)               tCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACaaa                                                                                                                                       >  2:855606/1‑149 (MQ=255)                ccTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACaaa                                                                                                                                       >  2:1047891/1‑148 (MQ=255)                   gAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                  >  2:174866/1‑150 (MQ=255)                   aaaccatcttCCGCGTTTTTTCCTAATTTT‑CGGCCGGGGGTGGGGTCCCCGTAAACACTCCGCCGGGGTTT‑GGCTTTCCCAGCCCTGCCAGCGGGCGCAGCAGGGTTTTTTT‑G‑CCACCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGAcc                                                                                                                                 <  2:785719/147‑1 (MQ=255)                                     cTTCCTAAATTT‑CAGCCGGGGGTGGGGCCCCCGTAAACAACCCGCCGGGGGTT‑GGCGTTCCCAGCCCCCCCAGGAGGCGAAGCGGGGTTTTTTT‑C‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                               <  2:1500447/151‑1 (MQ=255)                                       ttcttaaTTTCCAGCCGGGGGTGGGGTCCCCGCTAAC‑CCTCCCCCGGGGTTGGGCTTTCCCAGCCCCGCCAGCAGGCGAAGCGGGGTTTTTTT‑C‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              <  2:889309/146‑1 (MQ=255)                                       ttccggatatccagccgggGGGGGGGTGCCCGGTAACCAACCCCCTGGGGTG‑GCCTTTACCACACCGGCCAGCGGGCGCAGCAGGGTTCTTTA‑G‑CCCAAA‑CCGCTGCGGCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              <  1:48824/140‑1 (MQ=255)                                           gAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGccc                                                                                                                                                     <  1:92247/107‑1 (MQ=255)                                           gAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGccc                                                                                                                                                     >  2:92283/1‑107 (MQ=255)                                                tcccgccggggttgtgGGGCCCGCAACCCCACGCCCGGCGGTG‑GCCTTTCCCAGCCCCCCCAGCGGGCGCAGCGGGGTTTTTTG‑C‑CCACCC‑CCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     <  1:666456/139‑1 (MQ=255)                                                        gcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGGACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                             >  2:1012398/1‑151 (MQ=255)                                                                gtgCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                     >  2:569438/1‑151 (MQ=255)                                                                   cccGCAACCACACCCCCTGGGTTG‑GGCTTTCCCAGCCCGCCCGGCGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCC‑CCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                  <  2:124829/151‑1 (MQ=255)                                                                   cccGCAAACACACCGCCGGGGGTT‑GGCTTTCCCAGCCCGCCCAGCGGGCGCGGGGGGGTTTTTTT‑C‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                  <  2:699073/151‑1 (MQ=255)                                                                     cGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                >  2:1033036/1‑151 (MQ=255)                                                                     cGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                >  1:201966/1‑151 (MQ=255)                                                                     cGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTCGGTTCAGGACGATATTTTCCGCGTaa                                                                                >  1:1439885/1‑151 (MQ=255)                                                                       cTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                              >  2:1149307/1‑151 (MQ=255)                                                                           cacaTCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          >  2:659728/1‑151 (MQ=255)                                                                             caTCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGAGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          >  1:37995/1‑149 (MQ=255)                                                                                      cGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                               >  2:998411/1‑151 (MQ=255)                                                                                       cgttgggcttccccAGCCCCGCCAGCAGGCGCAGCAGGTTATTTTT‑C‑CCACCC‑CCGCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                <  1:1032528/138‑1 (MQ=255)                                                                                          tgggCTTTCCCAGCCCCCCCAGCAGGCGCAGCGGGTTTTTTTC‑C‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTCGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                            <  2:1440616/149‑1 (MQ=255)                                                                                                ttgcacgACCGGCCAGCGGGGGCAGGAGGGTTTTTTT‑A‑CCACCT‑CCTCGCCGCCACCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCATACTGCATTGAGcaacaa                                                      <  1:1127798/145‑1 (MQ=255)                                                                                                ttcccagcccccccAGCGGGCGC‑GGAGGGTTCTTTT‑GCCCCCCC‑CCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                      <  2:3624/139‑1 (MQ=255)                                                                                                ttcccAGCCCGGCCGGCAGGCGCAGCGGGTTTTTTTT‑C‑CCCCCC‑CCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                      <  2:1120423/148‑1 (MQ=255)                                                                                                 ttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                     >  2:1168755/1‑151 (MQ=255)                                                                                                 ttCAAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                     >  1:803529/1‑151 (MQ=255)                                                                                                   ccagcccccccAGCAGGCGCAGCGGGGTCTTTTT‑C‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  2:56704/142‑1 (MQ=255)                                                                                                   ccaccccTCCCAGCGGGCGCGGGAGGTTTTTTTT‑C‑CCCCCC‑CCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  2:46885/146‑1 (MQ=255)                                                                                                   ccAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCC‑CCGCTGCGGCCCACCACCGCCAAAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  1:1038182/151‑1 (MQ=255)                                                                                                     acccccgCCAGCGGGCGCGGCGGGTTTTTTTC‑C‑CCCCCC‑CCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                 <  2:708013/145‑1 (MQ=255)                                                                                                       ccccgCCGGCGGGCGCAG‑GGGGTTCTTTTCC‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                 <  2:929553/145‑1 (MQ=255)                                                                                                       cccccccagagggCGCAGCGGGGTTTTTTC‑C‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                 <  2:1053826/138‑1 (MQ=255)                                                                                                           gcgagcggGCCCAGCGGGGTCCTTTG‑C‑CCACCA‑CCGCGGCGCCCCACACCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                            <  1:877043/147‑1 (MQ=255)                                                                                                                cgggCGCAGCGGGTTTTTTTC‑C‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt                                       <  2:567129/148‑1 (MQ=255)                                                                                                                 cgggccaagcgggTTCTTTT‑C‑CCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt                                       <  2:1265551/140‑1 (MQ=255)                                                                                                                    cgcAGCAGGCGACTTTA‑C‑CACCCA‑CGGCTGCGCCCCACCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  <  1:658967/151‑1 (MQ=255)                                                                                                                    cccggggggtttttttt‑C‑CCCCCC‑CCGCCGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                   <  2:1233667/137‑1 (MQ=255)                                                                                                                     ccggcgggttACTTTC‑C‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 <  1:471432/142‑1 (MQ=255)                                                                                                                      cacgggggttttttt‑GCCCACCC‑CCGCGGCGCCCCCCCCCCGCCACAAACTGACCCGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 <  1:527497/140‑1 (MQ=255)                                                                                                                        gcagGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                               >  1:670958/1‑150 (MQ=255)                                                                                                                        gcagGGCACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCTTGCCCTGGTTCCGACGAg                               >  1:1123362/1‑150 (MQ=255)                                                                                                                               aCTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCT‑GTTCAGACGAGCAGTATTCa                      >  1:153814/1‑151 (MQ=255)                                                                                                                                  ttt‑C‑CCACCC‑CCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa                    <  2:581192/151‑1 (MQ=255)                                                                                                                                   tt‑G‑CCCCCC‑CCGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:204633/151‑1 (MQ=255)                                                                                                                                    tgg‑CCCCCA‑CCCCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGGCGAGCAGTATTCATAc                   <  1:340861/149‑1 (MQ=255)                                                                                                                                     ttcccccca‑ccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  1:569190/145‑1 (MQ=255)                                                                                                                                     tgcccaccc‑CCGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:193669/148‑1 (MQ=255)                                                                                                                                      ttcccccccccccggcgcCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:1367863/137‑1 (MQ=255)                                                                                                                                      ttcccccccccccTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:1492764/139‑1 (MQ=255)                                                                                                                                      tggccccccccccggcgcCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:841141/137‑1 (MQ=255)                                                                                                                                      tgcccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:1231257/142‑1 (MQ=255)                                                                                                                                      tccccccccccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  2:882240/142‑1 (MQ=255)                                                                                                                                      tccccccccccGCTCCGCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                   <  1:1173037/142‑1 (MQ=255)                                                                                                                                      ccccccca‑ccGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                  <  1:322791/146‑1 (MQ=255)                                                                                                                                       ccccccccccGCTGCCCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                  <  2:356001/143‑1 (MQ=255)                                                                                                                                       cccacca‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                 <  1:940890/150‑1 (MQ=255)                                                                                                                                        cccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                 <  1:236505/144‑1 (MQ=255)                                                                                                                                                       gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  2:987606/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993768‑994063

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