Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,580,634	0	C	T	38.6%	86.0 / 18.2	59	E277K (GAA→AAA)	hsdS	specificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base C (6/29);  new base T (22/0);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.08e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATT  >  NC_000913/4580488‑4580784                                                                                                                                                    |                                                                                                                                                        tAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCa                                                                                                                                                        >  2:158951/1‑148 (MQ=255)  aaTCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGATTA‑TGTTACttt                                                                                                                                                    >  2:136955/1‑151 (MQ=255)   atctatttctaataatTTTATCATCATATAACAAATTTTTATTTTATAATTTTTTCAATAAACTACAAACACCATCAAATTCTAAACTTCCTTTATATCTAATAAATAAAACATCTCCATCTTGTAATTTGTGTGGGTTTA‑GTTCAGtt                                                                                                                                                     <  1:350077/139‑1 (MQ=255)       aGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACttt                                                                                                                                                    >  2:55430/1‑146 (MQ=255)         cTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑TTTCACttttttaa                                                                                                                                               >  2:668442/1‑144 (MQ=255)          tCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACTTTTTTAACa                                                                                                                                             >  2:722041/1‑150 (MQ=255)           cGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACa                                                                                                                                             <  2:924880/149‑1 (MQ=255)             aaTAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTTCttttttgatctttt                                                                                                                                         >  1:267792/1‑140 (MQ=255)             aaTAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACttttttaa                                                                                                                                               >  1:62087/1‑140 (MQ=255)              aTAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTGTCAATAACCCACAAACACCAACAAATTCTATACTTCCGTTTTAGCAAGTAAAGAAAAGATCTCCATCTTATAATTTGTGGGGGGTTA‑GTTCACTTTCTGAACA‑tt                                                                                                                                          <  1:597788/148‑1 (MQ=255)                      aTCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGa                                                                                                                               <  2:437369/151‑1 (MQ=255)                           gATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAatat                                                                                                                           <  1:33324/150‑1 (MQ=255)                             tataGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAAC‑TCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑TTTTACTTTTTTAACA‑ATTTAGAAAACGCAAATCGtt                                                                                                                       >  2:357499/1‑151 (MQ=255)                                aGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGtttt                                                                                                                     <  1:313294/151‑1 (MQ=255)                                 gCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑TTTTACTTTTTTAAAA‑ATCTAAAAAACCAAAATCGtttt                                                                                                                     >  2:1036373/1‑150 (MQ=255)                                  cAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGGTTACATTAACTTT‑TGAACC‑TTTTAGAAAcctgaattacct                                                                                                                       >  1:131811/1‑137 (MQ=255)                                        ttGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACTTTTTTAAAA‑TTCTAGAAACCGAAAATAGTTTTTATCt                                                                                                                >  2:694403/1‑148 (MQ=255)                                           aTGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTCACTTTTTTGAAA‑ATCTAGAAAACGAATATCGTTTTTAttt                                                                                                                >  2:991718/1‑143 (MQ=255)                                             gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACTTTTTTAACA‑TTCTAAAAAACGAATATCGTTTTTATTTACATgggc                                                                                                        >  2:49787/1‑149 (MQ=255)                                               tgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCa                                                                                                       <  2:272192/150‑1 (MQ=255)                                                 tAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑TTTTACTTTTTTTACA‑TTTTTTAAAACCAATATAGttttttt                                                                                                                  >  1:558978/1‑134 (MQ=255)                                                  aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTTCTTTCTTAAAA‑TTCTAGAAAACCAATATCGTTTTGATCTACATTGCCAGc                                                                                                     >  2:1065590/1‑149 (MQ=255)                                                  aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCAc                                                                                                   <  2:606966/151‑1 (MQ=255)                                                     tttttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTACTTTGTGGCGGCTTT‑CTTTACCCTCTCAAAA‑TTCTCAAAACCCAACATCTTTTTTATCTAcctcccctcc                                                                                                     >  1:337700/1‑137 (MQ=255)                                                       tttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTAc                                                                                               <  1:735254/150‑1 (MQ=255)                                                            aTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GGTTACTTTCTTAAAA‑TTCTAGGAACCGAATATTGTTTTGTTCTTCAAGGCCCGCACCGACAGaa                                                                                           >  2:248264/1‑149 (MQ=255)                                                                 ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTCCTTTTTTAACA‑ATTTTGAAAACGAATATCGTTTTTATCTTCATAGCCACCACCTACAGAACTAATgg                                                                                    >  2:812973/1‑150 (MQ=255)                                                                 ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATg                                                                                     <  1:191220/150‑1 (MQ=255)                                                                  cacaAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACTTTTTTAACA‑TTCTAGAAAACGAATATAGTTTTTATCTACATGGCCCGCACGTACAGAAATAATGCg                                                                                   >  2:262160/1‑151 (MQ=255)                                                                           caacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTgg                                                                          <  2:290625/151‑1 (MQ=255)                                                                                            ccGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACttt                                                         <  1:112873/151‑1 (MQ=255)                                                                                             cGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACTTTTTTAACA‑TTCTAGAAAACGAATATCGTTTTTATTTACATTGCCCGCACGTACAGAACTAAAGCGTAGGATTTGGTGACCAACAACACctt                                                         >  2:995788/1‑150 (MQ=255)                                                                                              gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACTTTTTGGACA‑TTTTAGAAACCGAATATCGTTTTTATCTACATGGCCAGCACGGACAAAAATAAAGCGTAGTATTTGATGACCAACACCACtt                                                          >  2:73070/1‑148 (MQ=255)                                                                                                 ataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACCTTTTTAACA‑TTCTTGAAAACGAATATTGTTTTTATCTTCATGGCCAGCACGGACAGAACTAATTCGGAGTATTGGATGACCAACACCACCTTTAttt                                                    >  2:673705/1‑151 (MQ=255)                                                                                                  taGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACTTTCTTAACA‑TTTTAGAAAACGAAAATCGTTTTTATCTACATGGCCAGCACGTACAGAAATAATGCGTAGTATTTGGTGACCAAAACCACTTTCAttt                                                    >  2:271757/1‑150 (MQ=255)                                                                                                   aGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GGTTAATTTTTTTAAA‑TTTTAAAAAACTAATATCGTTTTTTTTTATATCGTTAGCACGTAAAAAAATAATGGGAAGGTTTGGATTGCCAACACCACtcttttt                                                     >  1:34808/1‑142 (MQ=255)                                                                                                   aGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTgg                                                  <  2:511915/151‑1 (MQ=255)                                                                                                   aGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGCTTA‑TTTCCCTTTTTTAACA‑TTCTATAAACCCAAAATCGTTTTTATCTACATTGTCCGCACCGACAGAAATAAAGCGTAGTATTTGGGGAACAACACCCCttttaattt                                                   >  2:765505/1‑144 (MQ=255)                                                                                                      gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTTACTTTTTGAACA‑TTCTAGAAAACGAATATCGTTTTTATTTACATGGGCAGCACGTACAGAAATAATGCGTAGTATTTGATGACCAAAACCACTTTTATTTTGCtt                                               >  2:520074/1‑151 (MQ=255)                                                                                                      gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTT‑GTTGCCTTTTTGAACA‑TTCTAGAAAACGAACATCGTTTTTTTTTACATGGGCACCCCGTGCAAAACCAATGCGTAGTAATCGATAACCCACACCCCatttattt                                                    >  2:819720/1‑146 (MQ=255)                                                                                                      gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCtt                                               <  2:131856/151‑1 (MQ=255)                                                                                                      gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCtt                                               <  1:866472/151‑1 (MQ=255)                                                                                                       aGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCttt                                              <  1:680730/151‑1 (MQ=255)                                                                                                          aaataaaAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTtga                                            <  2:1124915/150‑1 (MQ=255)                                                                                                           aataaaAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACGCCACTTTCATTTGGCTTtgatg                                          <  1:695597/151‑1 (MQ=255)                                                                                                           aataaaAGATCTCCATCTTGTAAGTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTATAGAAACCGAATATCGTTTTGATCTACATTGCCAGCACTTACAGAACTAATGCGTAGTATTGGATGACCAACACCTCTGTCATTTGGCGTtgatg                                          <  1:136909/151‑1 (MQ=255)                                                                                                              aaaaGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGaaa                                       <  1:812565/151‑1 (MQ=255)                                                                                                                     ctcCATCTTGTAATTTGTGGCGGTTTT‑TTTTACTTTCTGAACA‑TTCTTGAAACCGCATATCGTTTTTATTTACATTGGCAGCACCCACAGAAATAATACGGAGTATTTGGTGAACCACACCCCTTTCATTTTGCtttttttgaaaaaccat                                >  2:184944/1‑137 (MQ=255)                                                                                                                     ctcCATCTTGTAATTTGTGGCGGTTTA‑GTTCAACTTTTTAACACTTCT‑CATAAAAAATATTGGTTTGATTTACATTGCCAATACATTCAGAAATAAAACGTAATTTTGGTTGACCACCACCACTTTCTTTTGTCTTTGTTTAAAAAACAtt                                >  1:124648/1‑151 (MQ=255)                                                                                                                       ccATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAGCTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTa                               <  1:927916/150‑1 (MQ=255)                                                                                                                       ccATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCa                                  <  2:923472/147‑1 (MQ=255)                                                                                                                             tGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTaa                          <  2:885298/149‑1 (MQ=255)                                                                                                                                    tgtgGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTa                  <  2:239032/150‑1 (MQ=255)                                                                                                                                      tgGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCGCTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAag                <  1:19370/150‑1 (MQ=255)                                                                                                                                         cGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagatag            <  1:450053/151‑1 (MQ=255)                                                                                                                                          ggTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagataga           <  2:1051073/151‑1 (MQ=255)                                                                                                                                          ggTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat              >  2:646834/1‑148 (MQ=255)                                                                                                                                              tA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACAAGGCCAGCGCGTACAGAACTAAGGCGTAGTATGGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTGAAGATAGAtt         <  1:997562/149‑1 (MQ=255)                                                                                                                                                 gTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCaaa     <  2:155885/151‑1 (MQ=255)                                                                                                                                                    cACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAAtt  <  1:508136/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        TAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTA‑GTTCACTTTCTGAACA‑TTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAAATT  >  NC_000913/4580488‑4580784

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