Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	44.8%	82.1 / 15.0	59	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (24/8);  new base C (0/26);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 1.40e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCG  >  NC_000913/907460‑907751                                                                                                                                                |                                                                                                                                                      cgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCAc                                                                                                                                               <  2:420790/151‑1 (MQ=255)       catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCa                                                                                                                                          <  1:30934/150‑1 (MQ=255)       catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCa                                                                                                                                          <  1:840836/150‑1 (MQ=255)          ccGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGGGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAACCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCAt                                                                                                                                         <  2:187403/148‑1 (MQ=255)           cGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATaaa                                                                                                                                      >  2:838359/1‑150 (MQ=255)                  agATCGCGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGGGGGCGGTTCGGAAGAAATCACATTGAGAAACCACACCGAAATCCAGCGGGAGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCa                                                                                                                               >  2:356199/1‑150 (MQ=255)                         gtTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTg                                                                                                                       >  1:627440/1‑151 (MQ=255)                          ttccAGCCCCTGAACCACCGCCCGGGGGGTGGTGGGGGGTACGAAATAAATCACATTGAAAACCCACCCCAAAACCCAGGGGGGGGGGAGGGGGTGGGGGAT‑CACCCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTg                                                                                                                       <  2:435584/147‑1 (MQ=255)                                ccccTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATa                                                                                                                 >  1:507902/1‑150 (MQ=255)                                     gAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                             >  2:873476/1‑149 (MQ=255)                                       accAACCCCCTGGCGGTGGTGGGGGGTACGAAAGTAACCACATAAAAAACCCCCACCAAAATCGAGGGGGTGGGGGGGGGGATGGGAAT‑CAACGGGCGGCGCTTCCCTCTGTCCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                                              <  1:96956/144‑1 (MQ=255)                                          aaCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                       >  1:103671/1‑150 (MQ=255)                                            cGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                    >  1:400953/1‑151 (MQ=255)                                            cGCCCGGGGGGTGGGGGGGGGAAGGAAATAAATCACATAAAGAACCCCCCCCGAAACCCGGGGGGTGGGGGGGGGGAGGGGGAT‑AACCCGGCGGCCCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                    <  2:838007/151‑1 (MQ=255)                                                  gggggtgggGGGCGGTAGGAAAGTAATCACAGTGAAAAACCACCCCAAAACCCGGCGGGGGGGGGGGGGGAGGGGGAT‑CAACCAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                                               <  1:824243/147‑1 (MQ=255)                                                      gtggtgGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                          >  1:861732/1‑151 (MQ=255)                                                      gggggggggggTAGGAAAGTAATCACATTGAGAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGAT‑CACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                          <  2:788412/146‑1 (MQ=255)                                                       tggtggGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCggcggc                                                                                         >  1:680263/1‑151 (MQ=255)                                                        ggtggGGGGTACGAAAGTAATCACATTGAGAACCCCCACCGAAACCCGGGGGGGGGGGGGGGGGAGGGGAAT‑CACCCGGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCt                                                                                        <  1:692611/151‑1 (MQ=255)                                                        ggggggCGGTACGGAAGAAACCACATTGAAAACCCCCACCAAACTCCGGGGGGGGGGGGGGGGGTGGGGAAT‑CACCCAGCGGCGCTCACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCt                                                                                        <  2:832638/148‑1 (MQ=255)                                                         gtggGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCtt                                                                                       >  1:472710/1‑151 (MQ=255)                                                          ggggCGGTAGGGAAGAAATCACATTGAAAACCCACACCGAAACCCGGGGGGGGGGGGGGGGGAGGGGGAT‑CAACCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                                      <  1:586052/150‑1 (MQ=255)                                                                gTAGGGAAGAAATCAGATGAAAACACCACACCAAATTCCAGGGGGTGGGGAGGGGGAGGGGAGT‑CAACAGGCGGCGCTCACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGaa                                                                                   <  2:504113/148‑1 (MQ=255)                                                                 ggggaaaatAACCACATGGAAAACCCCCCCCAAAGCCCAGGGGGGGGGGGGGGGGATGGGGAT‑AACCGGGCGGCTCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATg                                                                               <  2:507846/143‑1 (MQ=255)                                                                  aCGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGc                                                                              >  1:850578/1‑151 (MQ=255)                                                                   gtaaaGAAATCACTATAAAAAACCACGCCTAAACCCAGGGGATGGGGGGGGGTGAGGGTAT‑CAACGTGTGTCTTTCTCCCCTCACCGCCATAAAATGTTCAAACAGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                                             <  2:19962/148‑1 (MQ=255)                                                                    ggAAGTAACCACATTGAGAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGATGGGGAT‑CACCCGGCGGCTCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATg                                                                               <  2:379639/148‑1 (MQ=255)                                                                      agtacaTCACATGGAGAACCCCCCCGGAAAACGGGGGGGGGGGGGGGGGGGAGGAGAC‑CACCCAGCGGCGCGTCCCTCACCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGa                                                                                    <  2:727626/136‑1 (MQ=255)                                                                      aaatAATCAAATTGAAAACCCCCCCCAAAACCCGGGGGGTGGGGAGGGGGAGGGGGCT‑CACCCGGCGGCTCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                          <  2:604218/148‑1 (MQ=255)                                                                      aaGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                          >  1:524776/1‑151 (MQ=255)                                                                       aGTAATCACATGGAAAAACCACCCCAAAATCCGGGGGGGGGGGGGGGGGAGGGGAAT‑CAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtgt                                                                         <  1:161270/151‑1 (MQ=255)                                                                        gTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGtt                                                                        >  2:620398/1‑151 (MQ=255)                                                                                  tGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                              >  2:479039/1‑151 (MQ=255)                                                                                  ggacaaaCCCCCCCGAACTCCGGCGGGGGGGGAGGGGGATGGGGAT‑CACCCAGGGGCTCTCCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                                                <  2:231539/145‑1 (MQ=255)                                                                                    aaaaaCCCCACCAAACCCCAGGGGGGGGGGGGGGGGTGGGGGAT‑CACCCGGGGCCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCt                                                                    <  2:223706/141‑1 (MQ=255)                                                                                    aaaaaCCCCACCAAAAACCGGGGGGGGGGGGGGGGGAGGGGGAT‑CAACCGGGGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                            <  2:452155/149‑1 (MQ=255)                                                                                            ccccaaaccccggggggTGGGGGGGGGGAGGGGGAT‑CAACCGGCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaa                                                     <  1:434102/136‑1 (MQ=255)                                                                                            cacCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGGGTTCGCTGGCCGGATGATCTa                                                      >  2:749596/1‑149 (MQ=255)                                                                                                 agaccgggggggTGGGGAGGGGGTGGGGAAA‑CAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaac                                                 <  1:558989/140‑1 (MQ=255)                                                                                                 aaccccAGGGGGGGGGGAGGGGGAGGGGGAC‑AACCCGGCGGCTCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                                <  1:228276/146‑1 (MQ=255)                                                                                                 aaaCCCAGCGGGGGGGGAGGGGGA‑GGGGATACAACCAGCGGCTCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaa                                                     <  1:604940/145‑1 (MQ=255)                                                                                                   aTCCAGGGGGTGGGGGGGGGGATGGGGAT‑CAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCaa                                             <  1:773267/151‑1 (MQ=255)                                                                                                     ccggggggTGGGGGGGGGGAGGGGGAT‑CAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAc                                            <  1:834910/145‑1 (MQ=255)                                                                                                     ccAGGGGGGGGGGAGGGGGATGGGGAT‑CAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                              <  1:218568/148‑1 (MQ=255)                                                                                                      cAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                           >  2:568981/1‑150 (MQ=255)                                                                                                            ggggggggggGTAGGGGGAC‑CACCGAGCGGGTTTTCCCTCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCGACCAg                                         <  2:95892/144‑1 (MQ=255)                                                                                                                     ggagtggggac‑CACCCGGGGGCCCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGt                           <  1:236150/147‑1 (MQ=255)                                                                                                                      gggtgggaat‑CACCCGGCGGCTTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                          <  1:146521/148‑1 (MQ=255)                                                                                                                       gatggggat‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                         >  1:604577/1‑151 (MQ=255)                                                                                                                         tggggat‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTc                       >  1:511369/1‑151 (MQ=255)                                                                                                                          gggtat‑CACCCGGCGCCGTTTCCCCCCCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                      <  2:623012/151‑1 (MQ=255)                                                                                                                           gggaaacacccggcggcgCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                      <  2:641203/139‑1 (MQ=255)                                                                                                                              at‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcg                  >  1:493582/1‑151 (MQ=255)                                                                                                                                 accccAGCGGCGCTACCCTCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcg                  <  1:93695/146‑1 (MQ=255)                                                                                                                                  aaCCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa                >  1:346285/1‑150 (MQ=255)                                                                                                                                   aCCAACGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAAt              >  2:789222/1‑151 (MQ=255)                                                                                                                                          gcgcTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtg       >  1:64600/1‑151 (MQ=255)                                                                                                                                          gcgcTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAACACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgtg       >  1:507936/1‑151 (MQ=255)                                                                                                                                               tACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCg  >  1:425558/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                      CGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCG  >  NC_000913/907460‑907751

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: ‑