Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	55.4%	50.0 / 99.1	65	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (3/26);  new base G (36/0);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 3.64e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-11

CTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941567‑941859                                                                                                                                                   |                                                                                                                                                     ctactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGaat                                                                                                                                                  >  2:148985/1‑149 (MQ=255)    ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                 >  1:1138708/1‑149 (MQ=255)    ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAcct                                                                                                                                                 >  2:1170911/1‑147 (MQ=255)      cctccGCGCAAATTGCGGAAGGTTTGAACTATACCCACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACt                                                                                                                                              <  2:579282/150‑1 (MQ=255)         ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACt                                                                                                                                              >  2:93633/1‑147 (MQ=255)          cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTAtttt                                                                                                                                         >  1:979285/1‑151 (MQ=255)           gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATAttt                                                                                                                                        >  1:105102/1‑150 (MQ=255)               aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTaa                                                                                                                                      >  2:333689/1‑149 (MQ=255)                aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAg                                                                                                                                   <  1:382387/151‑1 (MQ=255)                    gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATCTTTAACAgg                                                                                                                                  >  2:806938/1‑148 (MQ=255)                        aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCa                                                                                                                          <  1:723987/151‑1 (MQ=255)                             ttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGaaa                                                                                                                      <  2:1017650/150‑1 (MQ=255)                               gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGC‑ACGCCCGaaa                                                                                                                      >  2:584272/1‑148 (MQ=255)                                   tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTAATATCTTTAAAAAGC‑ACGCCAGAAAACTaa                                                                                                                 >  1:806854/1‑149 (MQ=255)                                    ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATg                                                                                                               <  2:738457/150‑1 (MQ=255)                                     taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTAATATTTTGAACAAGCTACCCCAGAAA‑CTAATGc                                                                                                              >  1:31435/1‑150 (MQ=255)                                       ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCg                                                                                                           <  2:875507/151‑1 (MQ=255)                                          aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGc                                                                                                          <  2:182990/149‑1 (MQ=255)                                          aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCa                                                                                                         <  2:4667/150‑1 (MQ=255)                                            ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTTAACAAGC‑ACGCCAGAAAACTAATGCCCGCCTg                                                                                                       >  2:952246/1‑150 (MQ=255)                                               ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTAATATTTTTAAAAGGC‑ACGCCCGAAAAATAATGCCCGCATg                                                                                                       >  1:718969/1‑147 (MQ=255)                                                 cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAAGC‑ACGCCAGAAATCTAAAGCCCCCATGatcctc                                                                                                 >  1:104340/1‑151 (MQ=255)                                                     gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                             <  2:454168/151‑1 (MQ=255)                                                     gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                             <  1:821738/151‑1 (MQ=255)                                                     gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGC‑ACGCCAGAAAACTAATGCCCGCATGatcatcaac                                                                                              >  1:904762/1‑148 (MQ=255)                                                      cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGatc                                                                                                    >  1:883199/1‑143 (MQ=255)                                                                   cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑Acc                                                                              >  2:1128773/1‑151 (MQ=255)                                                                    ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCCTATTATTTTTAACAAGC‑ACCCCCCAAATTTAATGCCCCCCACTTCATTATCCATCCGCGCTTT‑Acc                                                                              >  1:5128/1‑150 (MQ=255)                                                                     cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGa                                                                            <  1:930065/151‑1 (MQ=255)                                                                           atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCgg                                                                      <  2:402850/151‑1 (MQ=255)                                                                            tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGC‑ACGCCCGAAATCTAATAACCACAGGATCATCATCGATCCGCGCCAT‑ACCGACAcc                                                                        >  2:136496/1‑148 (MQ=255)                                                                            tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCgg                                                                      >  2:426359/1‑150 (MQ=255)                                                                               cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGAGCCGCGCTAT‑ACCGACACCAGTGcc                                                                  <  1:1024306/151‑1 (MQ=255)                                                                                gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑AACCACACCCGTGcc                                                                  >  2:158376/1‑150 (MQ=255)                                                                                 aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTTTTTTTAAAAAGC‑ACGCCCGAAATCTAATACCCCCATGTTCTTCCTC‑ATCCCCGCCCTAACCGCCACCGGGGCggg                                                                >  2:1126701/1‑151 (MQ=255)                                                                                   aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGgc                                                              <  1:1112208/151‑1 (MQ=255)                                                                                      aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAAAAGGC‑ACGCCAGAAATCTAAAGCCCGCATGATCATCATCCATCCCCGCTAT‑ACCCACACCCGTGGCggg                                                               >  2:258214/1‑147 (MQ=255)                                                                                                 aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTAATATCTTTAAAAAGC‑ACGCCAGAAAACTAAAGCCCCCCTGAACATCCTCGATCCGCGCTAT‑ACCGCCACCCGGGCCGGGGGGGAAGATGAgggg                                                >  1:237200/1‑151 (MQ=255)                                                                                                    gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                              <  1:162724/150‑1 (MQ=255)                                                                                                    gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCATACCACCTTGCACAACG‑ACTCCACAAAATCAATGCCCCCCTGATCACCATCGTATCGGGCTAT‑TCCCACATCCCGGCCGGGGGCCAAGAAGGGTg                                                 >  1:1000500/1‑147 (MQ=255)                                                                                                      tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGC‑CCGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGGGCCGGGGGCGGAGATGAGTGGGTccc                                           >  2:689883/1‑151 (MQ=255)                                                                                                       ggttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATcc                                            <  2:177804/147‑1 (MQ=255)                                                                                                       gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTAATATCTTTAACAGGC‑CCGCCAGAAATCTAATGCCCGCATGATCATCATTGATCCGCGCCAT‑ACCGACACCCGTGCCGGGCGCGAAAATGAGGGGATccc                                           >  1:235936/1‑150 (MQ=255)                                                                                                         ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                         <  2:970936/150‑1 (MQ=255)                                                                                                          tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGC‑ACGCCAGAAAACTAATGCCCCCCTGATCATCATCGATCCGCGCTAT‑ACCGACAACGGGGCCGGGGGCGAAGATGAGGGGGTCCCTAtt                                       >  2:548357/1‑151 (MQ=255)                                                                                                           ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                       <  1:497316/150‑1 (MQ=255)                                                                                                              aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATTTTGAACAGGC‑ACGCCAGAAATTTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGGGCGGAGATGAGTGGATCCCTATTCGTc                                   >  2:193752/1‑151 (MQ=255)                                                                                                               ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGC‑ACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGGGCCGGGCGCGAAGATGAGGGGATCCCTATTCATcc                                  >  2:329456/1‑151 (MQ=255)                                                                                                                 aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCCACACCGGGGCCGGGCGCGGAGATGAGTGGATCCCTATTCGTCCgg                                >  2:520583/1‑151 (MQ=255)                                                                                                                   cccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAACAGGC‑ACGCCAGAAAACTAAAGCCCGCCTGATGATCCTCGATACGCGCTAT‑ACCCACACCGGGGCCGGGGGCGGAGATGAGGGGATCCCCATTCGTCCggg                               >  1:1139079/1‑150 (MQ=255)                                                                                                                       ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                            <  1:690944/149‑1 (MQ=255)                                                                                                                       ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                           <  2:915201/150‑1 (MQ=255)                                                                                                                          gAAACGCGAATGAGTGGCGGAGGGGGGGCTTATTCTCTTGAAAAGGC‑ACGGCAGAAATCTAAAGCCCGCCTGACCATCCTCGCTTCGCGCTCT‑AACGACACCGGGGCCGGGCGGGGAGAAGAGTGGATACCCTTTCTTCCCGGGACaaacg                       >  1:1038240/1‑147 (MQ=255)                                                                                                                            aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                     <  1:964207/151‑1 (MQ=255)                                                                                                                            aaCGCGAATGAGTGGCGGATGGGGGGCTTAATTTCTTTAAAAAGG‑ACCCCAGAAATCTTAAGCCCGCATGATCATCATCGATCCCCGCTTT‑ACCGCCACCTGGGCCGGGCGCGAAAAATTGTGGAGCCCTATTTTTCCCGGAACAGa                         >  1:498835/1‑147 (MQ=255)                                                                                                                             aCGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATGcc                     >  2:9132/1‑150 (MQ=255)                                                                                                                                cgAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCAT‑ACCGAAACCGGGGCCGGGCGCGAAGATGAGGGGGTCCCTATTCGTCCGGGAACAGAAGCCGccc                 >  2:970321/1‑149 (MQ=255)                                                                                                                                 gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                <  2:545220/151‑1 (MQ=255)                                                                                                                                 gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                <  1:397470/151‑1 (MQ=255)                                                                                                                                 gAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGGGGGTCCCCATTCGTCCCGGGAAAGAAGCCGc                   >  1:199787/1‑148 (MQ=255)                                                                                                                                      aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa           <  2:928531/151‑1 (MQ=255)                                                                                                                                      aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa           <  2:201938/151‑1 (MQ=255)                                                                                                                                       gTGGCGGTGGGGGGAATTACTATCTTGAAAAGGC‑ACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCCGGGCCGGGCGCGGAGATGAGTGGAACCCTATTCGTCCGGGGACAGATGCCGCCCTgggt            >  1:945931/1‑147 (MQ=255)                                                                                                                                            ggTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc     >  1:39657/1‑151 (MQ=255)                                                                                                                                               ggggTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg  <  2:998706/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGC‑ACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTAT‑ACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941567‑941859

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