Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	52.6%	74.6 / 103.8	57	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (25/2);  new base C (0/30);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 3.53e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-11

GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGT  >  NC_000913/1376356‑1376644                                                                                                                                                           |                                                                                                                                            gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAacac                                                                                                                                            >  2:122233/1‑151 (MQ=255)      ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCg                                                                                                                                       <  1:757300/151‑1 (MQ=255)           gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                  >  1:131623/1‑151 (MQ=255)              gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                               >  1:830331/1‑151 (MQ=255)                ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                             >  1:1049829/1‑151 (MQ=255)                ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                             <  2:1122148/151‑1 (MQ=255)                ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCccg                                                                                                                             >  1:378107/1‑151 (MQ=255)                             cgtGGACGGCGACCGAAAAATGGTGGGGTATCCCCGGCGCGTACT‑GAACGGGGCGGCCGCGCGTCCCGGAA‑CG‑GTGGTGTATCTCAACCCGTTAAC‑CCGC‑TATAAGGTCCATATAACCAACCCCCTCGCCGATGCCCGCCGTTa                                                                                                                       <  2:517772/142‑1 (MQ=255)                                ggccGGCAACCGAAAAATGGGGGGTGATCCCCTGCGCGTACTGAAACGGGTC‑CCCGCGCTTCCCGGAC‑CC‑GTGGTGTATCTCGAACCGTAAAC‑CCGC‑TATCAGGATAATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGaa                                                                                                             <  1:325583/148‑1 (MQ=255)                                    ggCGCCCGCAAAAGGTGGAGTGTTCCCCTGCGCGTACTGAAACGGGTC‑GCCGCGCGTCCCGGAC‑CG‑GTGGTTTATCCCGACCCTTTAAC‑CCCT‑TATCAGGTTCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                            <  1:664866/148‑1 (MQ=255)                                    ggCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                         >  2:380539/1‑151 (MQ=255)                                        aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                     >  2:150111/1‑151 (MQ=255)                                         ccgaaaaggggtgggtgTTCCCCTGCGCTTACGGGACCGGGCC‑GCCGCGGGTCCGGGAC‑CG‑GGGGTGTATCTCGACCCTTAAAC‑CCGC‑TATCAGGATAATATAATCAACCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                       <  2:516839/139‑1 (MQ=255)                                         ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                    >  2:962762/1‑151 (MQ=255)                                          cGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                    >  1:760801/1‑150 (MQ=255)                                                aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                             >  1:14974/1‑151 (MQ=255)                                                  ggTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTCCCGAAC‑CG‑GTGTTGTATCTCGACCGGTAAAC‑CCGC‑TACCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                             <  1:1130433/149‑1 (MQ=255)                                                  ggTGAGTGATCCCCTGCGCTTACTGGAACAGGTC‑GCCGCGCGTCCCGGAC‑CC‑GGGGTGTATCTCGACCCGTTAAC‑CCGC‑TATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                             <  1:255488/149‑1 (MQ=255)                                                     gagtgaTTCCCTGCGCGTACGGGAACAGGTC‑GCCGCGCGT‑CCGGAACCC‑GTGGTGTATCTCGACCCGTAAAC‑CCGC‑TATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                          <  1:1184414/149‑1 (MQ=255)                                                       gtgaTTCCCGGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAC‑CC‑GTGGTGTATCTCAACCCGTAAAC‑CCGC‑TATAAGAATAATATGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                      <  1:1100953/151‑1 (MQ=255)                                                          aTTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                   >  1:682864/1‑151 (MQ=255)                                                             cccTGCGCGTACGGGAACAGGTC‑GCCGCGCGTCCCGAAC‑CC‑GTGGTGTACCTCGACCCGTAAAC‑CCGT‑TATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                  <  1:725951/149‑1 (MQ=255)                                                             cccGGCGCGTACTGGAACAGGTC‑CCCGCGCGTCCCGGAA‑CC‑GTGGTGTATCTCGAACCTTTAAC‑CCGC‑TATCAGAATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                     <  1:382368/146‑1 (MQ=255)                                                                   gcgTACTGAAACGGGTC‑CCCGCGCGTCCGGGAA‑CC‑GGGGTGTATCTCAACCCTTTAAC‑CCGC‑TATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                           <  1:1192339/150‑1 (MQ=255)                                                                     gTACTGGAACAGGTC‑GCCGCGCGTACGGGAC‑CG‑GTGGTGTATCTCG‑ACCTTTAAACCCTC‑TATCAGGATCATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                         <  2:527689/150‑1 (MQ=255)                                                                      tacggaaaCGGGTC‑CCCGCGCGTCCCGGAC‑CG‑GTGGTGTATCCCGACCCTTAAAC‑CCGC‑TATCAGAATAATATAATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGa                                                                             <  2:1021041/139‑1 (MQ=255)                                                                      tACTGGAACGGGTC‑GCCGGGCGTCCCGGAC‑CC‑GGGTTGTATCCCGAACCTTAAAC‑CCGT‑TAACAGGATCATATGAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                                        <  2:334376/150‑1 (MQ=255)                                                                         tGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                     >  1:803511/1‑150 (MQ=255)                                                                         ggaaacggtgc‑gccgcGCGTACCGGAA‑CG‑GGGGTGTATCCCGAACC‑TTAAA‑CCCCTTATCAGTATAATAGAATCAACCCCCTCGCCTTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                       <  2:1034406/138‑1 (MQ=255)                                                                               aGGTC‑GCCGCGGGT‑CCGGAA‑CCGGGGGTGTATCTCGACCGGTAAAC‑CCGC‑TATCTGGATTATATGATCAACCCCCCCGCAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                  <  1:539729/147‑1 (MQ=255)                                                                               aGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                               >  2:313162/1‑150 (MQ=255)                                                                                      ccgcGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                         >  2:379444/1‑150 (MQ=255)                                                                                      ccgcGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                         >  2:955535/1‑150 (MQ=255)                                                                                       cgcgcgTCCCGAAC‑CG‑GTGGTTTCTCCCGACCCTTAAAC‑CCGC‑TATCAGGATCATATAAACACCCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                        <  2:559577/150‑1 (MQ=255)                                                                                         cgcgTCCCGGAC‑CC‑GTGTTGTTCCACAACCCTTTAAA‑CCTC‑TATCAGGATCATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                      <  2:648424/150‑1 (MQ=255)                                                                                          gcgccccGGAC‑CG‑GTGGTGTTCCCCGCCCCTTAAAC‑CCGT‑TATCAGAATAATATAACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                      <  2:57167/144‑1 (MQ=255)                                                                                             tcccGGAC‑CC‑GGTGTGTATCTCGACCCTTAAAC‑CCGT‑TATCGGGATCATATGACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAagaaga                                                  <  1:305346/148‑1 (MQ=255)                                                                                              gtaccaaaccggtggtgTATCAATAACGGTAAAG‑CCTA‑TACCAGAATTATAAAATCATCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                        <  2:961666/134‑1 (MQ=255)                                                                                                 ggAA‑CG‑GTGGTGTATCTCGAACCGTAAAC‑CCGC‑TATCAGGATCATATGAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                 <  1:181755/147‑1 (MQ=255)                                                                                                       cggtgtttttcctcGACCCGTAAAC‑CCGC‑TATCAGGATCATATGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                              <  2:496795/135‑1 (MQ=255)                                                                                                       cggtggtgTTTCTCAACCCGTAAAC‑CCGC‑TATCAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                           <  1:222722/147‑1 (MQ=255)                                                                                                           ggtgTATCTCGAACCGTAAAC‑CCGC‑TATCAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                      <  1:530141/150‑1 (MQ=255)                                                                                                            ttgtACCTCGACCCTTTAAC‑CCGT‑TATAAGAATCATATAACAAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                           <  2:50649/143‑1 (MQ=255)                                                                                                             gttttcccGACCCGTAAAC‑CCGC‑AATCAGGATAATATGAACAACCCCCGCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                   <  2:619648/144‑1 (MQ=255)                                                                                                                 tctcGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                               >  2:427410/1‑151 (MQ=255)                                                                                                                        cgtttAAC‑CCGT‑TATCAGGATAATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                               <  1:482267/141‑1 (MQ=255)                                                                                                                        ccGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                         >  2:1201032/1‑150 (MQ=255)                                                                                                                           ttAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                     >  1:628880/1‑151 (MQ=255)                                                                                                                            tAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                    >  2:1006399/1‑151 (MQ=255)                                                                                                                             aaa‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                   >  1:704339/1‑151 (MQ=255)                                                                                                                              aaccccc‑TATCAGGATAATATGACCACCCCCCTCGCGTATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                   <  2:647375/145‑1 (MQ=255)                                                                                                                               a‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                  >  1:1107990/1‑150 (MQ=255)                                                                                                                                         cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg         >  2:338614/1‑151 (MQ=255)                                                                                                                                          aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg        >  1:112313/1‑151 (MQ=255)                                                                                                                                          aGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg        >  1:247651/1‑151 (MQ=255)                                                                                                                                           ggATTATAAGAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt       <  2:871648/151‑1 (MQ=255)                                                                                                                                             aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc      >  2:1060224/1‑150 (MQ=255)                                                                                                                                                atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt  >  1:741303/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                            GCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTC‑GCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTTAAA‑CCGC‑TATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGT  >  NC_000913/1376356‑1376644

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