Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	298,960	0	A	C	39.7%	115.6 / 25.6	63	Y659D (TAT→GAT)	paoC	PaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (36/2);  new base C (0/25);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.44e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACAC  >  NC_000913/298810‑299105                                                                                                                                                          |                                                                                                                                                  cGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATa                                                                                                                                                   <  2:406895/151‑1 (MQ=255)    tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATAccc                                                                                                                                                >  2:703034/1‑151 (MQ=255)        cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGcc                                                                                                                                            >  2:822046/1‑151 (MQ=255)        cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGcc                                                                                                                                            <  1:278492/151‑1 (MQ=255)          ccGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCa                                                                                                                                           <  1:1081852/150‑1 (MQ=255)                     aTAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTa                                                                                                                                >  1:108058/1‑150 (MQ=255)                      tAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAAc                                                                                                                              >  1:713480/1‑151 (MQ=255)                       aGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACg                                                                                                                             >  1:1092258/1‑151 (MQ=255)                       aGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACg                                                                                                                             >  1:630059/1‑151 (MQ=255)                         ccGCGCCCACGCCGCACAGCCCCAGCCCCCCGACCCTTTTGCCCTTCATGGGGGAGGATTTGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTTGGGG‑ATACCCGAAGGACCCGG‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACg                                                                                                                             <  1:541486/149‑1 (MQ=255)                               tCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGcc                                                                                                                     >  2:115395/1‑151 (MQ=255)                                cccgccgcACGGGCCCAGCCCACCGACCCCTTTGCCCTTCATCGGGGAGGATAGGGGGTGGGTATCACCCAGAAAAATCCCCTCCTTTTTTGGG‑ATACCCGCAGGACCCGG‑CACCTC‑ATCCCCCGCCATATCGTGATTAACg                                                                                                                             <  2:708731/140‑1 (MQ=255)                                cccgccgcACAGGCCCAGCCCACCGCCCCCTTGGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCACCTCCTTTTTTGGG‑ATATCCGCATGACCCGG‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTa                                                                                                                        <  2:767188/145‑1 (MQ=255)                                  cgccgcCCGGGCCCAGCTCCCGGACCCTTTTGGCATTCCACGGGGAGGTTAGGGGGTCGGTTTCATCCAGAAAATTCCCCTCCGGTTTGGGA‑ATACCCGCATGAGCCGG‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTa                                                                                                                        <  2:743693/145‑1 (MQ=255)                                  ccccgcACGGGCCCACCTCACCGACACTTTTGCCCTTCAGCGGGGAGGTTAGGGGGTCGTTATCATCCAGGAAAACCCCCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                                                                                                  <  2:1074141/149‑1 (MQ=255)                                   gccgcACGGCCCCAGCCCCCCGACCCTTTTGGCCTTCAGGGGGGAGGATAGGGGGTGGATATCATCCAGAAAAACCCCCCCCTGTTTTGGG‑ATTCCCGCAGGACCCGC‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                    <  1:587039/148‑1 (MQ=255)                                    ccgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                                >  1:566645/1‑151 (MQ=255)                                     cgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg                                                                                                               >  1:1006068/1‑151 (MQ=255)                                       cacaGCCCCAGCCCCCCGACCCTTTTGGCTTTCATGGGGGAGGATAGGGGGTCGTTATCACCCAGAAAAATCCCCTCCTGTTTTGGA‑ATATCCGCATGAACCGC‑CCCCCC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                                <  1:290574/148‑1 (MQ=255)                                        acaGGCCCAGCCCCCCGACCCCTTGGGCCTTCATGGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAATCCCCCCCTGTTTTGGG‑ATATCCGCATGACCCGC‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg                                                                                                              <  2:604412/149‑1 (MQ=255)                                                aGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg                                                                                                    >  1:160672/1‑151 (MQ=255)                                                  cTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc                                                                                                  >  1:36513/1‑151 (MQ=255)                                                   tCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc                                                                                                  >  2:982013/1‑150 (MQ=255)                                                         acacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcct                                                                                           >  1:682803/1‑151 (MQ=255)                                                            ccTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCtcctcc                                                                                         >  1:634653/1‑150 (MQ=255)                                                                ggcccTTCAGGGGGGAGAATAGGGGGTCGTTATCACCAAGAAAAACCCCCCCCTTTTTGGGA‑ATACCCGCAGGACCCGC‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa                                                                                     <  1:189746/147‑1 (MQ=255)                                                                   ctttCATGGGGGAGGATAGGGGTTCGGTTTCACCAAGAAAAATCACCTCCTGTTTGGGT‑ATACCCGCATGACCCGC‑CCCCCC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGt                                                                                   <  1:828493/147‑1 (MQ=255)                                                                          cGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc                                                                          >  2:933216/1‑151 (MQ=255)                                                                           ggggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc                                                                          >  1:886074/1‑150 (MQ=255)                                                                           ggggAGAATAGGGGGTGGTTATCACCCAGAAAAATCCCCCCCTTTTTGGGG‑ATACCCGCATAACCCGC‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc                                                                          <  2:544158/150‑1 (MQ=255)                                                                            ggaaggtTAGGGGGTCGTTATCACCCAGAAAAACCCCCCCCTGTTTTGGG‑ATACCCGAATGACCCGC‑CCCTCC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                        <  2:556698/148‑1 (MQ=255)                                                                                ggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                                    >  2:799603/1‑151 (MQ=255)                                                                                   tatGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                  >  2:526620/1‑150 (MQ=255)                                                                                            gttATCACCCA‑GAAAAACCCCCCCTTTTTTGGGAATACCCGCA‑GAACCCG‑CCCCCCAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                         <  2:356000/148‑1 (MQ=255)                                                                                            ggTATCACCCAGAAAACCCCCCCCCTTTTTGGGG‑ATACCCGAATAACCCGC‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAtt                                                           <  1:803056/148‑1 (MQ=255)                                                                                            ggTATCACCCAGAAAAACCCCCTCCTGTTTTGGG‑ATACCCGCAGGACCCGC‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                         <  2:367385/150‑1 (MQ=255)                                                                                               atcatcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                      >  1:886280/1‑150 (MQ=255)                                                                                                tcatcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCa                                                     >  1:164738/1‑150 (MQ=255)                                                                                                accaccgagAAAACCCCCCCCTTTTTTGGG‑AATCCCCCAAGGACCGGCCCCCTC‑ACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                      <  1:763527/142‑1 (MQ=255)                                                                                                 catcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                   >  1:868004/1‑151 (MQ=255)                                                                                                         aaaccccccccgttttttGGA‑ATCCCCGCATAACCCCC‑CCCCTC‑AGCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAcc                                               <  1:773896/133‑1 (MQ=255)                                                                                                          aaaTCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCt                                          >  2:243607/1‑151 (MQ=255)                                                                                                             tCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                       >  2:1040111/1‑151 (MQ=255)                                                                                                             tCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                       >  1:17977/1‑151 (MQ=255)                                                                                                              ccccccctGTTTGGGT‑ATACCCGCATGACCGGC‑CCCCCC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                         <  1:15345/143‑1 (MQ=255)                                                                                                                 ctcctGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                    >  1:110511/1‑150 (MQ=255)                                                                                                                   cctGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGt                                 >  1:977118/1‑151 (MQ=255)                                                                                                                    ctGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                                >  2:30546/1‑151 (MQ=255)                                                                                                                    ctGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                                >  1:949452/1‑151 (MQ=255)                                                                                                                    ctGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGt                                 >  2:70351/1‑150 (MQ=255)                                                                                                                     tGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCACCTCCTCCATCAGGGCCGCGCCCACGCCCATAGTCATTGCGCCAATGACCTAGCTGCGCGCAGt                                 >  1:139263/1‑149 (MQ=255)                                                                                                                     tGTTTGGGA‑ATACCCGAAGGACCCGC‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                                <  2:1041530/150‑1 (MQ=255)                                                                                                                                   ccGCATAACCCGC‑CCCCTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAg                     <  2:108076/148‑1 (MQ=255)                                                                                                                                    cGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg                 >  2:908456/1‑151 (MQ=255)                                                                                                                                            ccGG‑CGCCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg          >  1:425901/1‑150 (MQ=255)                                                                                                                                            ccGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc         >  1:423409/1‑151 (MQ=255)                                                                                                                                             ccgcccccTA‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc            <  2:616206/142‑1 (MQ=255)                                                                                                                                             ccgcccccTA‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt           <  1:522273/143‑1 (MQ=255)                                                                                                                                              gg‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc         >  1:1034160/1‑149 (MQ=255)                                                                                                                                              cgcccccTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgca        <  1:510221/146‑1 (MQ=255)                                                                                                                                               g‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagc      >  1:116175/1‑151 (MQ=255)                                                                                                                                                  cccTC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagca     <  1:170769/149‑1 (MQ=255)                                                                                                                                                     tC‑ATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac    <  1:468277/148‑1 (MQ=255)                                                                                                                                                     tC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac  >  1:432761/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                                  CGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CACCTC‑ATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACAC  >  NC_000913/298810‑299105

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