Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	503,834	0	T	C	43.1%	74.6 / 23.5	59	D288G (GAT→GGT)	fsr	putative fosmidomycin efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (20/13);  new base C (0/25);  total (20/39)
Fisher's exact test for biased strand distribution p-value = 4.51e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAGAAATA  >  NC_000913/503685‑503979                                                                                                                                                        |                                                                                                                                                     ccAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑c                                                                                                                                                    >  1:803255/1‑151 (MQ=255)     aaTGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCt                                                                                                                                                <  1:731606/151‑1 (MQ=255)        ggCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGAAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACa                                                                                                                                            <  2:878719/151‑1 (MQ=255)         gCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAg                                                                                                                                           <  2:191125/151‑1 (MQ=255)         gCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAg                                                                                                                                           <  1:597300/151‑1 (MQ=255)          cAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAgg                                                                                                                                          >  2:927785/1‑151 (MQ=255)            gagaATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑cc                                                                                                                                       >  1:223269/1‑151 (MQ=255)             agaATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑ccc                                                                                                                                      <  1:1222931/151‑1 (MQ=255)                aTGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑Cccgcc                                                                                                                                   <  2:388527/151‑1 (MQ=255)                    cGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATc                                                                                                                               >  2:1052661/1‑151 (MQ=255)                    cGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATc                                                                                                                               >  2:445174/1‑151 (MQ=255)                              aTAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTg                                                                                                                     >  2:258752/1‑151 (MQ=255)                               tAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCt                                                                                                                      >  1:863060/1‑149 (MQ=255)                                aaaTCCAATAATCACCGTTAAACCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGcc                                                                                                                   <  2:525388/151‑1 (MQ=255)                                aaaTCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGc                                                                                                                    <  2:69355/150‑1 (MQ=255)                                aaaTCCAATAACCCCCGTTAAACCCCCCGTCCAGTGCAGGCTGGCGAAGGGTAAAATCAGCGAAAACGGCGCAACGCCGAGGAAAGAGCCCCAAATCACATATT‑CCCGCCCAA‑TTGTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGcc                                                                                                                   <  1:916192/151‑1 (MQ=255)                                    ccAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCa                                                                                                                >  2:810887/1‑150 (MQ=255)                                          atCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCa                                                                                                          >  1:998607/1‑150 (MQ=255)                                                   aaaaCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAgg                                                                                                >  2:1067422/1‑151 (MQ=255)                                                    aaaCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                               >  2:366930/1‑151 (MQ=255)                                                       cccccGCCCAGTGCGGGCGGGCGGAGGGTAAAATCGCCGAAACCGGCCAACCGCCGAGGATAGAGCCCAAAATCACATATT‑CCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCaaa                                                                                            <  2:1006615/151‑1 (MQ=255)                                                       cccccGCCCAGGGCGGGCTGG‑GAAAGGGAAAAACAAAGTAAACCGCCACACCCCGAGGAATGGGCCCCAAAACCCCTTTTACCCCCCCAA‑TTTTAC‑CCCCG‑AGAGG‑CCCGCCGATCCCCATACCTGCCGCAACGGCAAACAGGAAGGc                                                                                               <  1:650888/148‑1 (MQ=255)                                                         cccGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAAGCa                                                                                              >  2:1214850/1‑147 (MQ=255)                                                          ccGTCCGGTGCGGGCGGGCTTGGGGAAAAATCGGAGAAAGCGGCACAGCCCGGAGAATAGGGCCCAAAAACACATTTT‑TCCCCCAAA‑TTTTAC‑CCCCG‑ACGGG‑CCCGCCGAACACCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAg                                                                                         <  2:374843/151‑1 (MQ=255)                                                          ccGCCCAGGGCAGGCGGGGGTAGGGAAAAACCAGCGAAAACGGCGCAACCCCGAGGATAGACCCCAAAACCACATTTT‑CCCGCCCAT‑TTTTAC‑CCCCT‑ACAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAg                                                                                         <  2:755000/151‑1 (MQ=255)                                                              ccAGTGCGGGCGGGCGTGGGGTAAAACCAGCGTAAACGGCGCAACGCCGAGGATAGACCCCAAAACCACATATT‑CCCCCCAAA‑TTTTAC‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGAt                                                                                       <  2:604375/149‑1 (MQ=255)                                                                   gCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCAGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCt                                                                                 >  1:982953/1‑150 (MQ=255)                                                                      ggcggggggaGGGAAAAATCAGCTAAACCGGCCCACCGCGGAGGAGAGACCCCAAAATCACAT‑TT‑TCCCCCCAA‑TTTTTTCCCCCT‑ACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGc                                                                             <  2:803283/142‑1 (MQ=255)                                                                       gttGGCGTGGGAAAAAACCAGGGAAACCGGCCCACCCCCGAGAATAGACCCCAAAACCACATATC‑CCCCCCCAT‑TTTTCC‑CCCCA‑AAGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGa                                                                               <  1:796962/146‑1 (MQ=255)                                                                        ctgcagtggggTAAATACACCGTAACCGGCGCAGGGCCGGGAATGGCTCCCAAAGCCGCATTTT‑CCCGCCCAT‑TTTTCG‑CCCCA‑ACGGC‑CCCGCGGATCGCCATTCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGc                                                                                  <  2:816658/136‑1 (MQ=255)                                                                          ggCGTAGGGAAAAATCAGCGAAACCGGCGCAACCCCGAGGATAGACCCCAAAAACAAATATT‑CCCCCCAAT‑TTTTAC‑CCCTT‑ACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTc                                                                         <  2:761151/151‑1 (MQ=255)                                                                            cGTGGGGAAAAACCGGGGAAAGCGGCGCAACGCGGAGGAAAGAGCCCAAATACACAT‑TT‑CCCCCCCCA‑TTTTTC‑CCCCTAACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCt                                                                        <  2:1017430/150‑1 (MQ=255)                                                                                   aaaaagcggctaaaaCGGCCCACCGCCGAGGAAAGACCCCAAAACCACATATT‑CCCGCCAAT‑TTTAAC‑CCCCA‑AAAGC‑CCCGCCGATCGCCGTGCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGgataga                                                                 <  1:425459/138‑1 (MQ=255)                                                                                               aacgggcacaccccggAGGAAAGACCCCAAAATCACAT‑TT‑CCCCCCCAATTTTTAC‑CCCCA‑ACAGC‑CCCGCCGATCACCGTACCTGCCGTAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAtttt                                                    <  1:1248382/136‑1 (MQ=255)                                                                                               aaccggccaaccGCCGAGGATAGCCCCCAAAAACACATATT‑CCCCCCAAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCTTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATcc                                                           <  2:190532/133‑1 (MQ=255)                                                                                               aaccGGCGCACCGCGGAGAATAGACCCCAAAACCACATTTT‑CCCCCCAAA‑TTTTAC‑CCCCA‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAttt                                                     <  2:70852/147‑1 (MQ=255)                                                                                                   ggCGCACCGCCGAGGAAAGACCCCAAAATCCCATTTT‑CCCGCCCAT‑TTTTAC‑CCCCA‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTg                                                  <  1:345432/149‑1 (MQ=255)                                                                                                        acccccGAGGATAGACCCCAAAATCACATTTT‑CCCGCCAAA‑TTTTAC‑CCCCT‑AAGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCAt                                               <  2:77014/143‑1 (MQ=255)                                                                                                          cGCCGAGAAAAGCGCCCAAAACAAAATTTT‑CCCCCCAAA‑TTTTAC‑CCCCA‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGAt                                          <  1:712526/150‑1 (MQ=255)                                                                                                             cGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATaaa                                       >  1:154051/1‑150 (MQ=255)                                                                                                              gAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATaaa                                       >  2:358906/1‑149 (MQ=255)                                                                                                                  aTAGACCCCCAAACCACATATT‑CCCCCCAAA‑TTTTAC‑CCCCG‑ACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAg                                     <  2:748449/147‑1 (MQ=255)                                                                                                                   tAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGgtg                                  >  2:384783/1‑149 (MQ=255)                                                                                                                       gCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAAt                              >  2:1219549/1‑149 (MQ=255)                                                                                                                          ccAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctg                         >  2:484591/1‑151 (MQ=255)                                                                                                                            aaaCCACATTTT‑CCCCCCCAA‑TTTTAC‑CCCCA‑AGAGG‑CCCGCCGATCACCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctgc                        <  2:580719/150‑1 (MQ=255)                                                                                                                            aaaCCAAATATT‑CCCCCCAAT‑TTTTCC‑CCCCA‑ACAGG‑CCCGCCGATCACCTTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctg                         <  1:736982/149‑1 (MQ=255)                                                                                                                             aacccccaatT‑TCCCCCCAA‑TTTTTT‑CCCCA‑AGTGGCCCCCCCGATCCCCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctgc                        <  1:1153416/142‑1 (MQ=255)                                                                                                                                cacaTATT‑CCCCCCAAA‑TTTTAC‑CCCCT‑ACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGAt                    <  2:617200/150‑1 (MQ=255)                                                                                                                                  caTATT‑TCCGACCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCt                 >  1:709449/1‑151 (MQ=255)                                                                                                                                   cttcttcccccccAT‑TTTTCC‑CCCCT‑CCACC‑CCCCCCAATCCCCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATg                   <  1:299125/139‑1 (MQ=255)                                                                                                                                       t‑tCCCCCCAATTTTTCC‑CCCCA‑ACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCg               <  1:282819/149‑1 (MQ=255)                                                                                                                                       t‑tCCCCCCAATTTTTAC‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTc                <  2:1166503/148‑1 (MQ=255)                                                                                                                                         ccccccaaa‑TTTTAC‑CCCAT‑ACGCC‑CCCTCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCt                 <  1:434708/138‑1 (MQ=255)                                                                                                                                           ccccaaa‑TTTTCC‑CCCCC‑ACAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCatat         <  1:287716/146‑1 (MQ=255)                                                                                                                                              caaa‑TTTTCC‑CCCCA‑AAAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAGa      <  2:95432/149‑1 (MQ=255)                                                                                                                                               cAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTTGATAGATAATCCGAATTTTTGCATCAGTTAAAAGGAGTAATAGCTGCTGATGCTCGCCATATAGa      >  1:463040/1‑150 (MQ=255)                                                                                                                                                aa‑TTTTAC‑CCCCT‑ACGCG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAg       <  2:92231/148‑1 (MQ=255)                                                                                                                                                    tttAC‑CCCCT‑ACGGC‑CCCGCCGATCACCGAACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAGAAata  <  2:689177/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                     CCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATT‑TCCGCCCAA‑TTTTAT‑CCCCT‑ACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCCATATAGAAATA  >  NC_000913/503685‑503979

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