Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	47.1%	70.8 / 18.9	51	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (5/22);  new base G (24/0);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 8.69e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366378‑366674                                                                                                                                                       |                                                                                                                                                      aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGggtgg                                                                                                                                                     <  2:673956/151‑1 (MQ=255)   gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGtt                                                                                                                                                   <  1:781343/151‑1 (MQ=255)    ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGGGCCCGGGGGGgtt                                                                                                                                                  >  1:560157/1‑151 (MQ=255)                aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTGAc                                                                                                                                        <  2:792667/149‑1 (MQ=255)                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGGGGGTTTTTTTTTCACCAGTg                                                                                                                                   >  1:835584/1‑151 (MQ=255)                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGAA‑TTGGGGCCCGGGGGGGTTTTTTTTTTTAccgag                                                                                                                                    >  2:724822/1‑146 (MQ=255)                     gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                                 <  2:756629/151‑1 (MQ=255)                        agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GCC‑TTGGGGGCCCGGGGGGGTTTTTTTTTTACCacagag                                                                                                                                 >  2:822410/1‑143 (MQ=255)                        agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                                <  2:370959/149‑1 (MQ=255)                          cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTG‑TTTGGGGCCAGGGGGGGTTTTTTTTTTTCCCGGGAGA‑Cggg                                                                                                                           >  1:291750/1‑151 (MQ=255)                            aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTAACAGTGAGA‑Cggg                                                                                                                           >  1:517661/1‑149 (MQ=255)                                      aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGAt                                                                                                                <  1:619227/150‑1 (MQ=255)                                       aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGGGGCCGGGGGGGTTTTCTTTTTTCCCGTGAGA‑AGGGcca                                                                                                                        >  2:586867/1‑139 (MQ=255)                                             ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGC‑TTGGGCGCCCCGGGGGTTTTTTTTTTTTCCCGTGagg‑caggccagacccggat                                                                                                               >  2:647403/1‑133 (MQ=255)                                              tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTc                                                                                                        <  1:609781/150‑1 (MQ=255)                                              tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTc                                                                                                        <  1:646488/150‑1 (MQ=255)                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCa                                                                                                       <  1:808898/150‑1 (MQ=255)                                                  cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTT‑GCA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGAGAGA‑CGGGCAACAGCTGATTGTCCCTTAcccc                                                                                                   >  1:909167/1‑149 (MQ=255)                                                    cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGcc                                                                                                  <  1:322715/150‑1 (MQ=255)                                                        gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGc                                                                                              <  2:549575/150‑1 (MQ=255)                                                        gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGc                                                                                              <  2:355049/150‑1 (MQ=255)                                                           cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTg                                                                                          <  1:16100/151‑1 (MQ=255)                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGT‑TTGGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGA‑CGGGGGACAGCTGATTTCCCTTCCCCGCCTGGccc                                                                                            >  2:618332/1‑148 (MQ=255)                                                               ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGGGGCAGGGGGGGTTTTTTTTTTCCCAGAGGGA‑CGGGGAACAACTGATTTGCCTTTACCCCCCGGCCCCgag                                                                                        >  1:794644/1‑149 (MQ=255)                                                               ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGA‑TTGGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGA‑CGGGCAACAGCCGGTTGGCCTTCACCGCCTGGccc                                                                                            >  2:944098/1‑145 (MQ=255)                                                                   aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGA‑TTTGGCGGCCGGGGGGTTTTTTTTTTCACCCGCGAGG‑AGGGGAAAAAGTTATTGCCCTTTCCCGCCCGGccc                                                                                            >  2:91654/1‑141 (MQ=255)                                                                    gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GGA‑GTTGGGGGCACGGGGGGTTTTCTTTTTACCAGGGGGA‑TGGGGAAAACCCGGTTGGCCCTCCCACCCTGGccccggcggcgat                                                                                  >  2:534790/1‑140 (MQ=255)                                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTgcagca                                                                            >  2:677468/1‑150 (MQ=255)                                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCCGGGTGGGTTTTTTTTTTACCAAGGAGA‑CGGGCAACAGCGGATTGCCCCTCACCCCCCGGCCCCGggagag                                                                                    >  2:854988/1‑142 (MQ=255)                                                                            aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                         <  2:819270/151‑1 (MQ=255)                                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTCCCAGTGAGA‑CGGGCAACAGACCATTGCCCCTCACCCCCCGGCCCCGAGAGAGTTGCAGGAAACgg                                                                       >  1:673141/1‑150 (MQ=255)                                                                               tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GGA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCCGGGAGA‑CGGGGAACAGCTGATTGGCCCTCCCCCCCCGGCCCCGCGAGAGGTGGAGCAAGCgg                                                                       >  2:45577/1‑150 (MQ=255)                                                                                   gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTATTTGGGGCCCAGGGGGGTTTTT‑TTTTTTCCAATGAGA‑CCGGCAAAAGCTGATTGCCCTTCACCGCCCGGCCCTGGGGGGGTTGCAGCAAACGGTcc                                                                    >  1:500362/1‑149 (MQ=255)                                                                                    tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGG‑GTA‑CTGTGCCCCCGGGGGGGTTTTTTTTTTTCCATAGAGA‑CGGGCCACAGCCGGTAGCCCTTTTCCCCCCGGCCCCTGGGGAGGTGCAGCAAGAGGGcc                                                                    >  2:852692/1‑148 (MQ=255)                                                                                      ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                               <  2:425551/151‑1 (MQ=255)                                                                                        aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGA‑GGA‑TTGGGGGCCAGGGGGGGTTTTTTTTTTACCAGGGGGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTTGGCCTGGGAGAGTTGCAGCAAACGGTCCCCGCCgg                                                             >  1:87582/1‑151 (MQ=255)                                                                                             cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTG‑GTA‑TTTGGGGCCACGGTGGGTTTTTTTTTTACCCGTGTGA‑AGGGCAACAGCTGATTTCCCTTCACCTCCCTGCCCCGAGGGAGTTGCAGCGAGAGGGcccccccgggtt                                                          >  1:776237/1‑139 (MQ=255)                                                                                               ttAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTC‑GTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCCGCAAGCGGGCCACCCCGGGTTGccc                                                      >  1:434935/1‑151 (MQ=255)                                                                                                    gAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagc                                                  <  2:663782/150‑1 (MQ=255)                                                                                                      aTCGGCCAACGCGCGGGGAGAGGCGGTTTTT‑GTA‑ATGGGGGCCAGGGGGGTTTTTTTTTTTCCCAAGGGGAGCGGGC‑CGAGCCCATTTGCCTTTACCGCCGGGGCCCGAGGGAGTTGCCGGAAGCGGTGCACGCTGGTTTTcccc                                                     >  2:224615/1‑145 (MQ=255)                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGa                                            <  2:355097/151‑1 (MQ=255)                                                                                                           ccAACGCGCGGGGAGAGGCGGTTTTC‑GCG‑TTGGGGGCCAGGGGGGGTTTTCTTTTCACCACCGGGA‑TGGGCAACACCGGATTTCCCTTCCCCCCCTGGCCCCGAGGGAGATGCAGCAAACGGTCAACGCCGGGTTGGCCCAGCAGGCGaa                                           >  1:183886/1‑150 (MQ=255)                                                                                                            cAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCCGGGGGGGTTTTTTTTTCACCAGGGGGA‑AGGGCAACAGCTGATTGGCCCTCCCCGCCCGGCCCTGGGGGAGGTGCAGCAAACGGGCCCCGCCGGGTTGGCCCCGCAgggggaa                                          >  2:589224/1‑145 (MQ=255)                                                                                                                   cgGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                                   <  2:322673/150‑1 (MQ=255)                                                                                                                      ggAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                               <  2:764778/151‑1 (MQ=255)                                                                                                                                   gC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                  <  1:224077/151‑1 (MQ=255)                                                                                                                                          ttGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              >  1:127420/1‑150 (MQ=255)                                                                                                                                            gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat            <  1:132003/150‑1 (MQ=255)                                                                                                                                              gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              >  1:741525/1‑146 (MQ=255)                                                                                                                                               cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata         >  2:826718/1‑150 (MQ=255)                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    <  2:813569/150‑1 (MQ=255)                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTAAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    <  1:348892/150‑1 (MQ=255)                                                                                                                                                     ggtggtTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg  <  2:460801/151‑1 (MQ=255)                                                                                                                                                     ggtggtTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg  <  2:417138/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGC‑GTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGA‑CGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366378‑366674

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