Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,543,440	T→C	7.3%	N30S (AAC→AGC)	cysP ←	thiosulfate/sulfate ABC transporter periplasmic binding protein CysP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,543,440	0	T	C	7.3%	198.1 / 5.4	68	N30S (AAC→AGC)	cysP	thiosulfate/sulfate ABC transporter periplasmic binding protein CysP
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (3/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAA  >  NC_000913/2543294‑2543588                                                                                                                                                   |                                                                                                                                                     ccTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTca                                                                                                                                                   <  1:694416/150‑1 (MQ=255) ccTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTca                                                                                                                                                   >  2:652624/1‑150 (MQ=255)    gTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTcagca                                                                                                                                                <  2:633449/150‑1 (MQ=255)    gTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTcagca                                                                                                                                                >  1:797878/1‑150 (MQ=255)      aaaaTCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGt                                                                                                                                              >  1:1795431/1‑150 (MQ=255)          tCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCg                                                                                                                                          <  1:1819025/150‑1 (MQ=255)            gCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGtt                                                                                                                                        >  1:1342239/1‑150 (MQ=255)            gCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGCTGAACAGTTcttat                                                                                                                                        <  1:498208/150‑5 (MQ=255)            gCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGCTGAACAGTTcttat                                                                                                                                        <  2:292232/150‑5 (MQ=255)               aGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGcc                                                                                                                                     <  2:1795605/150‑1 (MQ=255)               aGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGcc                                                                                                                                     <  1:1450325/150‑1 (MQ=255)                  gcCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGt                                                                                                                                  >  2:709806/1‑150 (MQ=255)                   ccTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTa                                                                                                                                 <  1:916225/150‑1 (MQ=255)                      gTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACAt                                                                                                                              >  1:1843739/1‑150 (MQ=255)                       ttttGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATg                                                                                                                             >  2:1167538/1‑150 (MQ=255)                                cccGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCcag                                                                                                                    <  2:797878/150‑1 (MQ=255)                                   ggCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCcagcag                                                                                                                 >  2:1376907/1‑150 (MQ=255)                                    gCAGGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCcagcagc                                                                                                                <  1:1576218/150‑1 (MQ=255)                                       tGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCaga                                                                                                             <  1:1441994/150‑1 (MQ=255)                                        gagaTTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCagag                                                                                                            <  1:1585680/150‑1 (MQ=255)                                           aTTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAg                                                                                                         <  2:1487587/150‑1 (MQ=255)                                           aTTGTTTAATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACCTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGCCCCGCCAGCAGCAGAGaac                                                                                                         >  2:596459/1‑149 (MQ=255)                                            ttGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGc                                                                                                        <  2:116677/150‑1 (MQ=255)                                              gTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGa                                                                                                      >  1:1845106/1‑150 (MQ=255)                                                tttATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGAcc                                                                                                    <  2:584057/150‑1 (MQ=255)                                                tttATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGAcc                                                                                                    <  1:13671/150‑1 (MQ=255)                                                   aTCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCCGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAgc                                                                                                 >  2:780662/1‑150 (MQ=255)                                                   aTCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTACGTTGCCTGTACATGGCCCCCCAGCAGCAGAGAAGCGACCAgc                                                                                                 >  2:63526/1‑150 (MQ=255)                                                     cGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAgcgc                                                                                               >  1:1360190/1‑150 (MQ=255)                                                     cGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAgcgc                                                                                               >  1:1988351/1‑150 (MQ=255)                                                            ttGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCgagtgag                                                                                        <  2:18270/150‑1 (MQ=255)                                                            ttGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCgagtgag                                                                                        <  1:6102/150‑1 (MQ=255)                                                                  ccgccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGCCCAGCGCGAGTGAGTTCttt                                                                                  >  1:138894/1‑150 (MQ=255)                                                                     ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGACCTGTTCAGCAGTTCCGTTGCCTGTACATGGCCGGCCAGCAGCAGAGAAGCGACCAGCGCGAGGGAGTTCTTTTTc                                                                               <  2:1021706/150‑1 (MQ=255)                                                                     ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                               <  2:1845106/150‑1 (MQ=255)                                                                     ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                               <  1:446182/150‑1 (MQ=255)                                                                     ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                               <  1:1376907/150‑1 (MQ=255)                                                                     ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                               >  1:1798715/1‑150 (MQ=255)                                                                     ccgTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTc                                                                               <  1:172426/150‑1 (MQ=255)                                                                                 gCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACg                                                                   <  1:652624/150‑1 (MQ=255)                                                                                 gCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACg                                                                   >  2:1088579/1‑150 (MQ=255)                                                                                       tgttgCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCAtt                                                             <  1:1776814/150‑1 (MQ=255)                                                                                        gttgCTCACTCGGCGGATTCAGGGCGGCAAACAGCTCGGGGGAGACGTCATAAGAAGTGTTCAGCAGTTCCGTTGCCTGTTCATGGCCGGCGAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCGATTg                                                            <  2:1652571/150‑1 (MQ=255)                                                                                            cTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCa                                                        >  1:208667/1‑150 (MQ=255)                                                                                              cAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCAcc                                                      >  2:72510/1‑150 (MQ=255)                                                                                              cAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCAcc                                                      >  2:1924813/1‑150 (MQ=255)                                                                                                aaCGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCt                                                    <  2:766894/150‑1 (MQ=255)                                                                                                     cggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATaa                                               <  1:1088579/150‑1 (MQ=255)                                                                                                     cggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATaa                                               >  2:728273/1‑150 (MQ=255)                                                                                                      ggATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATaaa                                              <  2:1174485/150‑1 (MQ=255)                                                                                                           cAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTaa                                         >  2:1746359/1‑150 (MQ=255)                                                                                                                   cAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTc                                 >  1:775196/1‑150 (MQ=255)                                                                                                                   cAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTc                                 >  2:1793061/1‑150 (MQ=255)                                                                                                                    aaaCAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCGGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCt                                <  2:1342239/150‑1 (MQ=255)                                                                                                                    aaaCAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCt                                <  2:1798715/150‑1 (MQ=255)                                                                                                                          cTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGc                          >  1:167824/1‑150 (MQ=255)                                                                                                                            cgcgGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCca                        >  1:1242789/1‑150 (MQ=255)                                                                                                                              cgGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCcatc                      <  2:1843739/150‑1 (MQ=255)                                                                                                                                       tCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTtata             <  1:1038314/150‑1 (MQ=255)                                                                                                                                        caacggAACTGCTGAACAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTCATAAATTTAATGACTTTCTAATAGCCATCATATTTATaa            >  1:1653688/6‑150 (MQ=255)                                                                                                                                         aacggAACTGCTGAACAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAAc           >  1:782079/5‑150 (MQ=255)                                                                                                                                         aTAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAAc           <  2:167824/150‑1 (MQ=255)                                                                                                                                          tAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACg          >  1:162645/1‑150 (MQ=255)                                                                                                                                          acggAACTGCTGAACAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACg          >  1:113986/4‑150 (MQ=255)                                                                                                                                           aaGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACgg         >  1:1975431/1‑150 (MQ=255)                                                                                                                                                cTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTg    >  1:559045/1‑150 (MQ=255)                                                                                                                                                 tGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGa   <  1:60193/150‑1 (MQ=255)                                                                                                                                                  gTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGaa  >  1:1057409/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CCTGTAAAATCGCCAGCGCCTGTTTTGATGACCCGGCATGAGATTGTTTTATCGTCAGTTTGTCGCCGCCGTTATCTTTTGCCCATTGTTGCTCAAACGGCGGATTCAGGGCGGCAAACAGCTCGCGGGAGACGTCATAAGAACTGTTCAGCAGTTCCGTTGCCTGTACATGGCCCGCCAGCAGCAGAGAAGCGACCAGCGCGAGTGAGTTCTTTTTCAGTAAGTTAACGGCCATTGCGCACCCTTATAAATTTAATGACTTTCTAATAGCCATCATATTTATAACGGACGTGAA  >  NC_000913/2543294‑2543588

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