Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	58.3%	51.1 / 76.1	60	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (19/6);  new base C (0/35);  total (19/41)
Fisher's exact test for biased strand distribution p-value = 8.66e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-09

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993768‑994064                                                                                                                                                           |                                                                                                                                                    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:93614/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                 <  1:1373058/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                 <  2:1066793/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                 <  2:1142422/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                 <  1:324505/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)                                                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                                                                                 >  1:201966/1‑151 (MQ=255)                                                                     cGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                 >  2:1033036/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)                                                                                                                                                       cccccccccCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt   <  2:347958/150‑1 (MQ=255)                                                                                                                                                       cccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt   <  1:997922/150‑1 (MQ=255)                                                                                                                                                       ccccACAACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt   <  1:1295755/150‑1 (MQ=255)                                                                                                                                                        ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg  <  2:1105035/151‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                    GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993768‑994064

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