Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	11,042	A→G	7.6%	V41V (GTA→GTG)	hutH →	histidine ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	11,042	0	A	G	7.6%	164.6 / 3.7	64	V41V (GTA→GTG)	hutH	histidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base G (2/3);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

CAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGTGAATCTGATACGCTCACATGCCTGTGGAC  >  CP000730/10904‑11178                                                                                                                                           |                                                                                                                                         cgaggagatAACATATATGAATTTATAGTTAGATGGTGAAACACTAAAAAGTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAAAAAAAAAAAAAAAAAGTAGAGCGGTAGt                                                                                                                                        <  2:71847/133‑1 (MQ=255)   aGGGGCTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTg                                                                                                                                      <  2:346304/141‑1 (MQ=255)    gggggTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGa                                                                                                                                     <  2:234998/141‑1 (MQ=255)    gggggTAAAATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAAAAAAAAAAAAAAAGTAGAGCGGTAGTTGa                                                                                                                                     <  2:234838/141‑1 (MQ=255)       ggTAACATATATGAGTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATGATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACg                                                                                                                                  >  2:529217/1‑141 (MQ=255)       ggTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACg                                                                                                                                  >  2:143819/1‑141 (MQ=255)        gTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGt                                                                                                                                 >  1:586901/1‑141 (MQ=255)           aCATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGtatt                                                                                                                              >  1:106856/1‑141 (MQ=255)           aCATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGGAGTTGAACGtatt                                                                                                                              >  2:225040/1‑141 (MQ=255)           aCATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGGTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGtatt                                                                                                                              >  1:106418/1‑141 (MQ=255)             atatatGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGtattat                                                                                                                            <  1:619392/141‑1 (MQ=255)                  tGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGaaa                                                                                                                       >  2:687134/1‑141 (MQ=255)                    aCTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAATTTAGGGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAAt                                                                                                                     >  1:178766/1‑141 (MQ=255)                     atttATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACTATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGCAGTTGAACGTATTATTGAAAATg                                                                                                                    <  2:553271/140‑1 (MQ=255)                           atTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGAAGTTGAACGTATTATTGAAAATGAGGaaa                                                                                                              >  2:311087/1‑141 (MQ=255)                               aGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGt                                                                                                          >  2:395318/1‑141 (MQ=255)                                 aTGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAATGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGAGGTAGTTGAACGTATTATCGATAATGAGGAAACGGttt                                                                                                        <  1:432533/141‑1 (MQ=255)                                 aTGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGttt                                                                                                        <  2:736495/141‑1 (MQ=255)                                  tGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTa                                                                                                       <  2:142491/141‑1 (MQ=255)                                     tgagtagcTAACAATTGCGGCTATTAAACCATTGTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACgg                                                                                                    <  1:363191/134‑1 (MQ=255)                                         acacTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTGGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTAtt                                                                                                >  1:725013/1‑141 (MQ=255)                                         acacTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTGGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTAtt                                                                                                >  1:724620/1‑141 (MQ=255)                                          cacCAACAATTAAGTAAATTAAATCATTTTTACAACAACAATGAAAGATTGAAATTATTGATGATGCGTTAGAACGAGTCAAAAAAAGTAGAGCGGTAGTTAAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATta                                                                                               <  1:131479/141‑1 (MQ=255)                                             tAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATtacta                                                                                            <  2:210755/141‑1 (MQ=255)                                             tAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATtacta                                                                                            >  1:304101/1‑141 (MQ=255)                                              aaCAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATtactac                                                                                           >  2:566115/1‑141 (MQ=255)                                                cAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTCTtactacgg                                                                                         >  1:316113/1‑139 (MQ=255)                                                aaaTTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGAGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGGATTATTGAAAATGAGGAAACGGATTACGGTATTACTACAg                                                                                         >  2:103354/2‑141 (MQ=255)                                                 aagtGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAgg                                                                                        <  1:473558/138‑1 (MQ=255)                                                 aaTTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAgg                                                                                        <  1:548468/141‑1 (MQ=255)                                                  aTTGAGGATATTAAATGCTTTTTACAACAACAATCAAAGATTGAAATTATTTATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGGGGAAACGGTTTACGGTATTACTACAGGt                                                                                       <  2:196244/141‑1 (MQ=255)                                                  aTTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGt                                                                                       >  1:212262/1‑141 (MQ=255)                                                       ggATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGATAGAACGAGTCAAAAAAAGTAGAGCGGAAGTGGAACGTATTATTGAAAATTAAGAAACGGTTTACGGTCTTACGACAGGATATgg                                                                                  >  2:775049/1‑141 (MQ=255)                                                             tAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACTTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTAtt                                                                            <  1:381982/141‑1 (MQ=255)                                                               aaTCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTa                                                                          <  1:225040/141‑1 (MQ=255)                                                                aTCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTAAAGGTTTTGGGTTATTTAg                                                                         >  2:693594/1‑141 (MQ=255)                                                                  cATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTg                                                                       <  2:334244/141‑1 (MQ=255)                                                                    tgtttACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTAAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGAt                                                                     <  2:686563/139‑1 (MQ=255)                                                                          caacgacaaTCAAAGATTGAAATTATTGATAATGCGTTAGAACGTGTCAAAAACAGTAGAGCGGGAGTTGAACGTATTATTGAAAATGAGGGAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACga                                                               >  2:471255/1‑140 (MQ=255)                                                                            acaacaaTCAATGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGGGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGtat                                                             <  2:263507/141‑1 (MQ=255)                                                                               acaaTCAAAGTTTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGa                                                          <  2:333716/141‑1 (MQ=255)                                                                               acaaTCAAAGGTTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACTTATAGa                                                          <  2:427793/141‑1 (MQ=255)                                                                                 aaccAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGAcc                                                        <  1:218375/138‑1 (MQ=255)                                                                                 aaTCAATGATTGAAATTATTGATGATGCGTTAGAACGCGTAAAAAAAAGTAGCGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGAcc                                                        <  1:72478/141‑1 (MQ=255)                                                                                  aTCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTGGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGAccc                                                       <  2:725013/141‑1 (MQ=255)                                                                                  aTCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTGGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGAccc                                                       <  2:724620/141‑1 (MQ=255)                                                                                  aTCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGCAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGAccc                                                       >  1:263603/1‑141 (MQ=255)                                                                                  aTCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGAccc                                                       <  1:371613/141‑1 (MQ=255)                                                                                   tCAAAGATTGAAATTATTGCTGGTGAGTTCGAAGGTACCATACAGCGTCGAGCTGTAGTTGATCTTACTGTTGAAAATGATGAGGCTGTTGACGGTAGTGCTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCg                                                      <  2:448648/141‑1 (MQ=255)                                                                                    cAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACGCGa                                                     >  1:361215/1‑141 (MQ=255)                                                                                      aaGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACg                                                   <  1:142933/141‑1 (MQ=255)                                                                                       aGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTCCGTATAGACCCGACGc                                                  >  2:407562/1‑141 (MQ=255)                                                                                               attattGATTATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATAATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAAt                                          <  2:479064/141‑1 (MQ=255)                                                                                                   ttGATGATGCGTTAGAACGTGTCACACCCAGCAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAAt                                      >  2:130488/1‑141 (MQ=255)                                                                                                      atgatgCGTTATAACGGGTCAAAAAAAGTAGAGCGGTAGTTTAACTTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTAc                                   <  2:305205/141‑1 (MQ=255)                                                                                                           gCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGTg                              >  1:68185/1‑141 (MQ=255)                                                                                                                   aCGTGTCAAAAAAAGTAGAGCGGTGGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAGTGAACTACAAGTGAATCTGAt                      <  1:395318/141‑1 (MQ=255)                                                                                                                       gtaaaaaaaaGTAGGGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGTGAATCTGATACGc                  <  2:734901/138‑1 (MQ=255)                                                                                                                          aaaaaaaGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGTGAATCTGCAACGCCca               >  2:388705/1‑141 (MQ=255)                                                                                                                           aaaaaaGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAGGTGAATCTGATACGCTcac              <  1:423477/141‑1 (MQ=255)                                                                                                                                 gTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACAGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGTGAATCTGATACGCTCACATGCCt        >  2:496521/1‑141 (MQ=255)                                                                                                                                  gagagCGGTAGTTGAACGTATTAGTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGTGAATCTGATACGCTCACATGCCtg       <  2:57158/140‑1 (MQ=255)                                                                                                                                    gagCGGTAGTTGAACGTATTGTTGAAAATGAGGAAACGGTTTACGGTATTACGACAGGATTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGAGAATCTGATACGCTCACATGCCtgtg     >  1:471933/1‑141 (MQ=255)                                                                                                                                       cGGTAGTTGAACGTATCATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGTGAATCTGATACGCTCACATGCCTGGGGAc  >  1:384265/1‑141 (MQ=255)                                                                                                                                           |                                                                                                                                         CAAGGGGGTAACATATATGACTTTATATTTAGATGGTGAAACACTAACAATTGAGGATATTAAATCATTTTTACAACAACAATCAAAGATTGAAATTATTGATGATGCGTTAGAACGTGTCAAAAAAAGTAGAGCGGTAGTTGAACGTATTATTGAAAATGAGGAAACGGTTTACGGTATTACTACAGGTTTTGGGTTATTTAGTGATGTACGTATAGACCCGACGCAATATAATGAATTACAAGTGAATCTGATACGCTCACATGCCTGTGGAC  >  CP000730/10904‑11178

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