Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	53.4%	86.7 / 119.4	73	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/39);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 1.10e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-10

ATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984296‑2984591                                                                                                                                                     |                                                                                                                                                       atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCt                                                                                                                                                    <  2:460816/150‑1 (MQ=255)   gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTgg                                                                                                                                                  <  2:287692/150‑1 (MQ=255)     cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGGCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGTAATGCAGAGCTGTTCGTCACATT‑CAATGGTCA‑ATG‑CCACACCCTGGcc                                                                                                                                                <  1:1114824/150‑1 (MQ=255)          gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCccgcc                                                                                                                                          >  2:600069/1‑151 (MQ=255)           aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCccgccg                                                                                                                                         <  2:1270729/151‑1 (MQ=255)            aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGa                                                                                                                                        <  1:1446644/151‑1 (MQ=255)                 aaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATAcg                                                                                                                                    >  1:824385/1‑149 (MQ=255)                  aaaTGAGGGAGCCAAGAAAATGTATGTCGCCTGCCCAGCCATCATTATTCAAAAGATAAATGACGGCAGATTTTTATCGGAAAGGAGAAAATGCAGAGCTTTTCTCCACGTT‑CAATGGTAA‑ATG‑CCACCCCCTGGCCCCCGCCGATACaaa                                                                                                                                  <  2:1263256/151‑1 (MQ=255)                           aGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGcc                                                                                                                          >  1:1249507/1‑150 (MQ=255)                                aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                    >  2:244817/1‑151 (MQ=255)                                  aaaaTGTATGTCGATTTCCCAGCCATCATTATACAACAAATTAATGACGGCATATATTTAGCGGTAAGGAGGAAATGCAGAGTTATTCGTCACTTT‑AATTGGCCA‑ACG‑CCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  1:213369/151‑1 (MQ=255)                                   aaaTTTATGTCGCTTTGCCAGCCATAATTATTCAACAGATTAAGGAGGGAAGATTTTTATCGGAAGGGAGAAAATGCAGAGTTTTTCGTCACGTT‑AAATGGTAA‑ATC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                                 <  1:416426/151‑1 (MQ=255)                                    aaTGTAGGTCGATTTGCCAGCCATCATTATACAAAAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGAAGAGCTTTTCGTCACGTT‑CAAGGGTAA‑ATC‑CCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                                <  1:1209632/151‑1 (MQ=255)                                             cgtcTTGCCAGCCATCATTATTCAACAGATAAAGGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGTTTTTCGCCACTTT‑CAAGGGCCA‑ATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                       <  2:384300/148‑1 (MQ=255)                                             cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                       >  1:669991/1‑151 (MQ=255)                                               aCTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                                     >  1:351981/1‑151 (MQ=255)                                                    cccgcaataatttttCAAAAAATAAAGGACGGAAGATATTTATGGGAAAGGAGAAAATGCAGAGCTTTTCGTCCGTTT‑CAATGGTAA‑ATC‑CCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                                 <  2:1439106/141‑1 (MQ=255)                                                    ccagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAc                                                                                                >  2:890968/1‑151 (MQ=255)                                                     cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               >  1:1234295/1‑151 (MQ=255)                                                      agccaTCATTATTCAACAGATAAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTT‑AAATGTTCA‑ATG‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               <  1:509553/150‑1 (MQ=255)                                                        ccaTCATTATTCAACAGATTAATGGCGGCAAATTTTTATCGGTAGGGAGGAAATGAAGAGCTTTTCGTCACGTT‑AAATGGTCA‑ATG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                             <  1:503650/150‑1 (MQ=255)                                                        ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                            >  1:1422935/1‑151 (MQ=255)                                                        ccaTCATTATTCAAAAGATTAATGACGGCAGTTATTTATGGGTAAGGAGGAAATGAAGAGCTTTTCGTCACGTT‑AAAGGGTCA‑ATC‑CCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                   <  1:722676/144‑1 (MQ=255)                                                        caatcatTATCAAAGAAATTAATACCGGCAGATTTTTATCGTTAAGGAGTAAAGGTAGTCTTATTGGTTCATTT‑CAATGGTAA‑CTC‑CCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                                 <  2:1425950/144‑1 (MQ=255)                                                         catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                           >  2:243623/1‑151 (MQ=255)                                                          atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          >  1:1069813/1‑151 (MQ=255)                                                          atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          >  1:53038/1‑151 (MQ=255)                                                           tcatTATTCAACAGATTAATGACGGAAAATATTTTTGGGTAAGGAGAAAATGCAGAGTTTTTCGTCCCTTT‑AAATGGTAA‑AGG‑CCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                                 <  1:461957/143‑1 (MQ=255)                                                             attattCAACAAATTAATGCCGGAAGATATTTATCGGTAGGGAGGAAATGCAGAGTTTTTCGTCCCGTT‑AAATGGTCA‑ATG‑CCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          <  2:1078621/148‑1 (MQ=255)                                                              ttatgaaaCAGATTAAGGCCGGCAAATATTTATCGGTAAGGAAGAAAAGGAGATTTTTTGGGC‑CGTTGAAAGGGTCA‑TGC‑CCACACCGTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               <  1:1189024/136‑1 (MQ=255)                                                                    aaaaGATAAATGAGGGAAAATATTAATGGGAAAGAGGAAAATGCAATCCTATTCGCCGCTTT‑CACTGTTAA‑ATC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                  <  2:568028/146‑1 (MQ=255)                                                                    aaCAGATTAATGCCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑AAAGGGTAA‑ATG‑CCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                      <  1:449253/145‑1 (MQ=255)                                                                    aaCAGATTAATGACGGCATATATTTATCGGTAAGGAGAAAATGAAGATCTTTTCGTCACTTT‑CAATGGTCA‑ATG‑CCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                   <  1:1424433/148‑1 (MQ=255)                                                                     aCAGATTAATGACGGCAGATTTTTATCGGTAAGGAGGAAATGCAGAGTTATTCTTCACGTT‑CAATGGTCA‑ACG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCagag                                                                               <  2:465176/151‑1 (MQ=255)                                                                       aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                              >  2:213703/1‑150 (MQ=255)                                                                          taaaggACGGAAGATATTTATGGGTAAGGAGAAAATGCAGAGTTTTTCGTCACTTT‑CAATGGTCA‑ATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                             <  2:336318/143‑1 (MQ=255)                                                                          taaaTGAGGGAAGATATTTATGGGTAAGGAGGAAAGGAAGATTTTTTCGTCACTTT‑AAATGGTAA‑ATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                              <  2:657311/145‑1 (MQ=255)                                                                            aaTGACGGCAGATTTTTATGGGAAGGGAGAAAATGAAGAGCTTTTCGTCACGTT‑CAATGGTCA‑AGC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                             <  2:816233/146‑1 (MQ=255)                                                                            aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                          >  2:559813/1‑149 (MQ=255)                                                                                aCGGCAGATTTTTATGGTAAAGGATAAAATGCAGCCTTTTTCGTCA‑GGT‑AAAAGGTCA‑AAGCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                     <  2:799750/150‑1 (MQ=255)                                                                                 cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACAAGAGAAACCAGGATTc                                                                   >  2:745803/1‑151 (MQ=255)                                                                                   gaagttttttATCGTAAAGGAGAAAATGCAGATTTTTTCTCCACTTT‑AAATGGTCA‑TGG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                                  <  2:991280/143‑1 (MQ=255)                                                                                    cAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                >  2:292553/1‑151 (MQ=255)                                                                                        tatTTATGGGTAGGGAGAAAATGCAGAGCTTTTCGTCACTTT‑CAATGGTAA‑ATG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                                 <  1:213799/146‑1 (MQ=255)                                                                                         atTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                           >  2:167668/1‑151 (MQ=255)                                                                                          tttACCGGAAAGGAAAAAATGCAGATTTTTTCTTCACTTT‑CATTGATCA‑AGC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa                                                              <  2:1212163/147‑1 (MQ=255)                                                                                              tgggTAGGGAGAAAAGGCAGATTTTTTCGTCACGTA‑CAATGGTAA‑ATC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa                                                        <  2:389182/147‑1 (MQ=255)                                                                                              acGGGAAGGAGAAAAGGCAGAGTTTTTCGTCACTTA‑AAATGGTAA‑ATC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                      <  2:515460/150‑1 (MQ=255)                                                                                                ggTAAGGAGGAAAGGCAGATTTTTTCGTCACTTC‑AAATGGTCA‑CTC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa                                                        <  2:812592/147‑1 (MQ=255)                                                                                                   aaGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                                    >  1:136444/1‑148 (MQ=255)                                                                                                        ggaAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                            >  1:1464700/1‑151 (MQ=255)                                                                                                          aaaTGCAGAGCTATTCGTCCTGTC‑CAATGGTCA‑ATG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                            <  1:548254/149‑1 (MQ=255)                                                                                                           aaTGCAGAGTTTTTCGTCCCGTA‑CAATGGTAA‑CGC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                                         <  2:1032492/151‑1 (MQ=255)                                                                                                            aTGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                        >  2:1125155/1‑151 (MQ=255)                                                                                                              gCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  2:716306/1‑151 (MQ=255)                                                                                                                 gagtttttCGTCCCGTT‑AAATGGTCA‑ATC‑CCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                                   <  1:282414/145‑1 (MQ=255)                                                                                                                 gagCTTTTCGTCCCGTT‑AAATGGTCA‑ATG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                                   <  2:810368/151‑1 (MQ=255)                                                                                                                    cttttCGCCACGTC‑CATTGCCAA‑ATG‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGTGCCTATCGGGTGACCGAGTg                                  <  2:1069712/146‑1 (MQ=255)                                                                                                                    cTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCACCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                >  1:715316/1‑151 (MQ=255)                                                                                                                       ttCGTCACGTT‑CAGTGCTGG‑CAC‑CCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                       <  1:130636/141‑1 (MQ=255)                                                                                                                        tCGTCACTTT‑AAATGGTAA‑ATC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             <  1:386086/150‑1 (MQ=255)                                                                                                                            caagTT‑CAATGGTTAGATC‑CCCCCCCCTGCCCACCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             <  2:1249408/144‑1 (MQ=255)                                                                                                                            cACGTT‑AAAGGGTAA‑ATG‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                            <  1:1414073/147‑1 (MQ=255)                                                                                                                            cACGTA‑AAATGGTCA‑AGC‑CCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             <  1:1020333/146‑1 (MQ=255)                                                                                                                                tt‑CATTGGTCA‑ATC‑CCCCCCCCTTGCCCCCGCCGATACGAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                                   <  2:824303/136‑1 (MQ=255)                                                                                                                                 t‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                   >  2:1003612/1‑151 (MQ=255)                                                                                                                                  taattGGTCA‑ATC‑CCCCCCCCCGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcca                       <  1:1226654/143‑1 (MQ=255)                                                                                                                                           a‑aTG‑CCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           <  1:256382/150‑1 (MQ=255)                                                                                                                                           a‑aTG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           >  1:1054547/1‑150 (MQ=255)                                                                                                                                           a‑aTG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           >  2:260443/1‑150 (MQ=255)                                                                                                                                             aTG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCattgacattgac         >  2:903675/1‑151 (MQ=255)                                                                                                                                              tG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc        >  1:318672/1‑151 (MQ=255)                                                                                                                                                 ccACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg      >  1:414481/1‑151 (MQ=255)                                                                                                                                                     acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  2:71243/1‑151 (MQ=255)                                                                                                                                                     acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  1:1055865/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       ATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTT‑CAATGGTCA‑ATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984296‑2984591

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