Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,748,198	G→A	6.2%	L356L (CTG→CTA)	ydhS →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,748,198	0	G	A	6.2%	209.3 / 3.4	64	L356L (CTG→CTA)	ydhS	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base A (2/2);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01

ATGTCGATTCCTTTGCTCAAGCCTGGTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGCCTTACGTGAAGCGGGA  >  NC_000913/1748056‑1748345                                                                                                                                               |                                                                                                                                                    aTGTCGATTCCTTTGCTCAAGCCTGGTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGaa                                                                                                                                               <  2:140376/149‑1 (MQ=255)    tCGATTCCTTTGCTCAAGCCTGGTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGcc                                                                                                                                            >  4:251034/1‑149 (MQ=255)         tCCTTTGCTCAAGCCTGGTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACa                                                                                                                                       <  3:224549/149‑1 (MQ=255)                  cAAGCCTGGTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGt                                                                                                                              >  3:462758/1‑149 (MQ=255)                    aGCCTGGTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTc                                                                                                                            <  3:324457/149‑1 (MQ=255)                    aGCCTGGTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTc                                                                                                                            <  3:465212/149‑1 (MQ=255)                          gTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATc                                                                                                                      >  3:625054/1‑149 (MQ=255)                              gCCGAGCGCAAACAACGCGACCCATTTGACTGAGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCtgaa                                                                                                                  <  2:255734/149‑1 (MQ=255)                               ccGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCtgaat                                                                                                                 <  3:628673/149‑1 (MQ=255)                               ccGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCtgaat                                                                                                                 <  1:433740/149‑1 (MQ=255)                                  aGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCtgaatgaa                                                                                                              <  2:603034/149‑1 (MQ=255)                                   gcgcAAACAACGCGACCCATTTGACTTGGCAGAAAAAAATCTCCAGGAAGGCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGTGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAac                                                                                                             <  2:160378/149‑1 (MQ=255)                                   gcgcAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAac                                                                                                             <  3:268478/149‑1 (MQ=255)                                       aaaCAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGa                                                                                                         >  1:648705/1‑149 (MQ=255)                                       aaaCAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGa                                                                                                         >  1:672449/1‑149 (MQ=255)                                       aaaCAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGGCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATGAGACTGAACACGa                                                                                                         >  2:540990/1‑149 (MQ=255)                                        aacaacGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGAc                                                                                                        <  2:186517/149‑1 (MQ=255)                                         acaacGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTACATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGAcc                                                                                                       >  1:582193/1‑149 (MQ=255)                                          caacGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCa                                                                                                      >  4:369099/1‑149 (MQ=255)                                            acGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAacac                                                                                                           <  1:677985/142‑1 (MQ=255)                                            acGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAacac                                                                                                           >  2:677985/1‑142 (MQ=255)                                                   ccATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGtt                                                                                             <  4:478327/149‑1 (MQ=255)                                                      tttGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTACATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGAcc                                                                                                       >  1:491341/1‑136 (MQ=255)                                                      tttGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTACATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGAcc                                                                                                       <  2:491341/136‑1 (MQ=255)                                                                 aGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTg                                                                               >  3:665987/1‑149 (MQ=255)                                                                    aaaaaaTCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGccc                                                                            <  1:656485/149‑1 (MQ=255)                                                                       aaaTCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCggg                                                                         <  4:403019/149‑1 (MQ=255)                                                                        aaTCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGt                                                                        <  4:549903/149‑1 (MQ=255)                                                                            tcCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTc                                                                    >  3:77070/1‑149 (MQ=255)                                                                               aGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATc                                                                 <  3:41936/149‑1 (MQ=255)                                                                                  aaGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGAt                                                              <  1:365431/149‑1 (MQ=255)                                                                                    gTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATaa                                                            <  2:470154/149‑1 (MQ=255)                                                                                     tCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAAt                                                           <  1:644615/149‑1 (MQ=255)                                                                                        aaCGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTAt                                                        >  3:331306/1‑149 (MQ=255)                                                                                          cGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGc                                                      <  4:462758/149‑1 (MQ=255)                                                                                          cGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTACATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGc                                                      <  2:582193/149‑1 (MQ=255)                                                                                             aaTAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGc                                                   >  2:246457/1‑149 (MQ=255)                                                                                             aaTAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGc                                                   >  3:273859/1‑149 (MQ=255)                                                                                                 aaCGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATc                                               >  3:370047/1‑149 (MQ=255)                                                                                                       aaaaaCATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCa                                         <  3:500876/149‑1 (MQ=255)                                                                                                       aaaaaCATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCa                                         <  3:245529/149‑1 (MQ=255)                                                                                                           aCATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGt                                     <  4:209880/149‑1 (MQ=255)                                                                                                              tACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTa                                  >  2:668054/1‑149 (MQ=255)                                                                                                              tACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTa                                  >  2:480274/1‑149 (MQ=255)                                                                                                                cTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTAtt                                >  1:559244/1‑149 (MQ=255)                                                                                                                cTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTAtt                                <  3:438188/149‑1 (MQ=255)                                                                                                                      ccTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGa                                                               >  3:313102/1‑112 (MQ=255)                                                                                                                      ccTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGa                                                               <  4:313102/112‑1 (MQ=255)                                                                                                                       cTGGCGTCATGTCATTCTGCTCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGCACACTACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGcc                         <  1:540990/149‑1 (MQ=255)                                                                                                                        tGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCt                        >  1:37218/1‑149 (MQ=255)                                                                                                                        tGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCt                        >  2:364820/1‑149 (MQ=255)                                                                                                                        tGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCt                        >  3:394669/1‑149 (MQ=255)                                                                                                                         ggCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTa                       >  4:292204/1‑149 (MQ=255)                                                                                                                              tATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGc                  <  4:665987/149‑1 (MQ=255)                                                                                                                              tATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGc                  <  1:274602/149‑1 (MQ=255)                                                                                                                              tATGTCAGTCTGCACCTGCATGAAGCCGTACAGGAAATTGTTCCACGTCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCAGCAATCCCTACAGAGTCTATTCGTCGCCTACTTGc                  <  2:228268/149‑1 (MQ=255)                                                                                                                                 gTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTa                                  <  4:291551/130‑1 (MQ=255)                                                                                                                                 gTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTa                                  >  3:291551/1‑130 (MQ=255)                                                                                                                                  tCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGCCTTa              >  1:262449/1‑149 (MQ=255)                                                                                                                                  tCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTGTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTGGCCttc              >  4:487453/1‑148 (MQ=255)                                                                                                                                   cATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGCCTTAc             <  3:433813/149‑1 (MQ=255)                                                                                                                                   cATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGCCTTAc             <  3:123243/149‑1 (MQ=255)                                                                                                                                             cTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGCCTTACGTGAAGCggg   >  4:570936/1‑149 (MQ=255)                                                                                                                                              tGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGCCTTACGTGAAGCGGGa  >  1:369991/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                    ATGTCGATTCCTTTGCTCAAGCCTGGTTTGCCGAGCGCAAACAACGCGACCCATTTGACTGGGCAGAAAAAAATCTCCAGGAAGTCGAACGCAATAAACGAGAAAAACATACTGTTCCCTGGCGTTATGTCATTCTGCGCCTGCATGAAGCCGTACAGGAAATTGTTCCACATCTGAATGAACACGACCATAAACGGTTCAGTAAAGGCCTTGCCCGGGTTTTCATCGATAATTATGCGGCAATCCCTTCAGAGTCTATTCGTCGCCTACTTGCCTTACGTGAAGCGGGA  >  NC_000913/1748056‑1748345

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