Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,601,176	0	T	C	43.3%	73.3 / 33.5	60	E304G (GAA→GGA)	ftsX	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (31/3);  new base C (0/26);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 5.23e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATG  >  NC_000913/3601026‑3601325                                                                                                                                                            |                                                                                                                                                      tattatTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACtt                                                                                                                                                       >  1:672303/1‑151 (MQ=255)  attattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACtt                                                                                                                                                       <  1:752408/150‑1 (MQ=255)   ttattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTcc                                                                                                                                                     >  1:759724/1‑151 (MQ=255)    tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTcc                                                                                                                                                     <  2:225739/150‑1 (MQ=255)    tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCg                                                                                                                                                    >  2:636920/1‑151 (MQ=255)      ttCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGc                                                                                                                                                   >  1:467910/1‑150 (MQ=255)         aGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAAcc                                                                                                                                               <  1:826707/151‑1 (MQ=255)                    ggCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGTCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACa                                                                                                                                     >  1:786092/1‑150 (MQ=255)                     gCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAAt                                                                                                                                   >  1:538709/1‑151 (MQ=255)                     gCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAAt                                                                                                                                   >  1:623604/1‑151 (MQ=255)                       gTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCCATAGCAGGCATTCATCGAATGATACGCCATTG‑ATCTCAAACCTCGCTCCGCAAA‑CC‑TGT‑C‑CCACTTCCGCAACCCCCGCCGACaacc                                                                                                                                  >  1:309746/1‑148 (MQ=255)                       gTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCg                                                                                                                                 >  1:795281/1‑151 (MQ=255)                        tAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCgc                                                                                                                                >  2:673755/1‑151 (MQ=255)                         aaaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCgc                                                                                                                                >  1:855956/1‑150 (MQ=255)                                gtACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATAACAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCAcc                                                                                                                          >  2:948557/1‑149 (MQ=255)                                 tACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCAcc                                                                                                                          >  1:625295/1‑148 (MQ=255)                                  aCCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGa                                                                                                                       >  1:191977/1‑150 (MQ=255)                                          aaGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGAc                                                                                                              >  1:612764/1‑151 (MQ=255)                                             ccccccTGCCCCCCACCAAACAACCGAGCATACCAGCATAAAGAGCGGGCATTCATCGAATGAAAAGCCATTA‑ATACAAAACTTCTTCCCAAAAC‑CC‑TGT‑C‑CCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaaa                                                                                                           <  2:960305/146‑1 (MQ=255)                                             ccccGCTGCCCCCCAGCAAATCATCGAGCACCCCAGCAGAAATAGCAGGCATTCATCAAATAATAAGCCATTA‑ATATCAAACTTCTTCCCAAAAC‑CC‑TGT‑C‑CCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTg                                                                                                                <  1:826180/143‑1 (MQ=255)                                              cgcGCTCCCGCCCGGCAAACCACCGAGCATCCCACAAGCAATAGCAGTCCTTCA‑AGAATTAAAAGGCCTTGAATACCAAGCTTCGTGCCAAAAC‑CC‑TCT‑C‑CCATCTCCACAGCCGCCGATGACAATCGCAGCACCAGAATTTCTGACaa                                                                                                            <  2:960932/147‑1 (MQ=255)                                                cgctccccccaaccaaaTCATCGAGCATACAAGAAAAAATGGCAGGATTTCATCAAATGAAAACCCTTTA‑ATATCAACCTTC‑TTCCCAAAA‑CC‑CGTCC‑CCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAt                                                                                                         <  2:551357/135‑1 (MQ=255)                                                           acccAATCACCGAGCATACCAACAACAAGAGCAGCCATTCAACAAATAATAACCAATTA‑ATATCAAACTCCTTCCCAAAAC‑CC‑TGT‑C‑CCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaaca                                                                                              <  2:65145/148‑1 (MQ=255)                                                             caaaTCACCGAGCATACCAGAAAAAATAGCAGGCATTCATCAAATGAAAAGCCTTTG‑ATAAAAACCTTCGTCCCAAAAC‑CC‑TGT‑C‑CCCTTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg                                                                                           <  1:788031/149‑1 (MQ=255)                                                               aaTCATCGAACACACCACAAAAAATAGCAGGCTTTCATAAAATAAAACGCCATTG‑ATATAAACCTTCGTCCCAAAAC‑CC‑TTT‑C‑CCTTTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaaca                                                                                              <  2:387976/146‑1 (MQ=255)                                                                   atggAGCACCCCAGAAGAAATAGGAGGCTTTCAGAAAATGATAAGCTATTA‑ATA‑AAAACCTCCGTCCCAAAA‑CC‑TGG‑GCCCCCTCCCGAAGCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAg                                                                                     <  2:432061/148‑1 (MQ=255)                                                                        gCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAAt                                                                                >  2:51719/1‑151 (MQ=255)                                                                          aTACCAGCAGAAATAGCAGGCTTTCATCGAATAAAAACCCATTA‑ATATCAACCTTCGTCCCAAAAC‑CT‑TGT‑C‑CCACTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAAt                                                                                <  1:897125/149‑1 (MQ=255)                                                                          aTACCAGCAGAAATAGCAGACATTCTTCAAATAAAAGCCCATTA‑ATATAAAACTCCTTCCCAAAAC‑CC‑TGT‑C‑CCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATc                                                                               <  2:146152/150‑1 (MQ=255)                                                                          aTACCAGCAGAAAGAGAAGACATTCACCAAATAAAAACCCTTTA‑ATAAAAACCTTCCTCCCAAAAC‑CC‑TGT‑C‑CCACTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaaaa                                                                                 <  2:720868/148‑1 (MQ=255)                                                                             caagaaaaaaTAGCAGCCTTTCACCAAATGAAAACCAATTA‑ATACCAAACTTCG‑TCCAAAAA‑CCTTGT‑C‑CCCCTCCCACATCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAAt                                                                                <  2:785186/138‑1 (MQ=255)                                                                                gcagcaATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcagc                                                                          >  2:432213/1‑149 (MQ=255)                                                                                 cagcaATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGt                                                                       >  1:50222/1‑151 (MQ=255)                                                                                 cagcaATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGt                                                                       >  1:655318/1‑151 (MQ=255)                                                                                  agaaaaaGCAGGTTTTCATGAAATGAAAAGCAATTA‑ATATCAAACTTCGTCCCAAAAC‑CC‑TGT‑C‑CCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCcagcag                                                                        <  1:641837/143‑1 (MQ=255)                                                                                    caATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgc                                                                     >  2:169757/1‑150 (MQ=255)                                                                                     aaaaGCAGGCATTCATCGAATAAAAACCCATTA‑ATACCAAACTCCTTCCCAAAAC‑CC‑TGT‑C‑CCACTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgc                                                                     <  2:759097/146‑1 (MQ=255)                                                                                         gCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGccacc                                                               >  2:817576/1‑151 (MQ=255)                                                                                            gtaaTTCATCGAATGATAAGCCATTG‑ATACAAAACTTCGTCCCAAAAC‑CC‑TGT‑C‑CCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGccacca                                                              <  2:821015/146‑1 (MQ=255)                                                                                            ggCATCCAACAAATAAAAACCCATTA‑ATATCAAACTTCGTCCCAAAAC‑CC‑TGC‑C‑CCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCAt                                                             <  2:493766/150‑1 (MQ=255)                                                                                                   aTCAAATGAAAAGCCATTA‑ATACCAACCTTCGTCCCAAAAC‑CC‑TGT‑C‑CCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAAc                                                     <  1:496426/151‑1 (MQ=255)                                                                                                       aataatAAGCTTTTA‑ATAT‑AAACCTCCTTCCAAAAACCC‑TGT‑C‑CCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACggg                                                  <  1:688606/150‑1 (MQ=255)                                                                                                        aTGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcg                                                >  1:185412/1‑151 (MQ=255)                                                                                                          gATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcgc                                               >  1:369537/1‑150 (MQ=255)                                                                                                              taccctattataaCAAACCTTCTTCCAAAAA‑CC‑CGTCG‑CCCTTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAg                                             <  2:623541/141‑1 (MQ=255)                                                                                                              aGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                                          >  1:34918/1‑151 (MQ=255)                                                                                                               acctatttataaaaaaCTTCGTTCCAAAAC‑CC‑TGT‑C‑CCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                                          <  1:744742/138‑1 (MQ=255)                                                                                                                    ttaatatAAAACTCCTTCCCAAAAC‑CC‑TGT‑C‑CCACTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGa                                        <  1:766851/145‑1 (MQ=255)                                                                                                                        tatCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATc                                 >  2:562026/1‑151 (MQ=255)                                                                                                                            aaccTTCGTCCCAAAAC‑CG‑TGC‑C‑CCACCCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCAt                                  <  2:863166/143‑1 (MQ=255)                                                                                                                             accttcgcccaaaaac‑CC‑TGT‑G‑CCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCg                            <  2:65840/139‑1 (MQ=255)                                                                                                                             aaCTCCGTCCCAAAAC‑CC‑TGG‑C‑CCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATcc                                    <  1:720974/143‑1 (MQ=255)                                                                                                                                 tCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCAcc                        >  1:616270/1‑151 (MQ=255)                                                                                                                                    tcccaaaac‑CC‑TGT‑C‑CCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCa                          <  1:703470/141‑1 (MQ=255)                                                                                                                                        gAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt                 >  2:778841/1‑151 (MQ=255)                                                                                                                                        aaaac‑CC‑TGT‑G‑CCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAg                  <  1:337503/149‑1 (MQ=255)                                                                                                                                         aaaa‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAGTCAGtt                >  1:510591/1‑151 (MQ=255)                                                                                                                                          aaa‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt               >  1:68118/1‑151 (MQ=255)                                                                                                                                            a‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCt             >  1:878611/1‑151 (MQ=255)                                                                                                                                                 tgt‑g‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGt           >  1:796989/1‑150 (MQ=255)                                                                                                                                                       ccACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGt       <  2:79356/150‑1 (MQ=255)                                                                                                                                                        cACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGtt      >  1:241710/1‑150 (MQ=255)                                                                                                                                                           ttCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATg  >  1:457979/1‑151 (MQ=255)                                                                                                                                                            |                                                                                                                                                      TATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTG‑ATATCAAACTTCGTTCCGAAAA‑CC‑TGT‑G‑CCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACGTTAATG  >  NC_000913/3601026‑3601325

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