Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,443,277	0	A	C	41.4%	93.1 / 22.8	58	V274G (GTA→GGA)	secY	preprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (30/4);  new base C (0/24);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.60e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443130‑3443427                                                                                                                                                     |                                                                                                                                                       gACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAccc                                                                                                                                                     >  1:630465/1‑151 (MQ=255)  aCGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTAcccc                                                                                                                                                    >  1:107289/1‑151 (MQ=255)      aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                                <  1:296224/151‑1 (MQ=255)          agTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCatat                                                                                                                                            <  2:52994/151‑1 (MQ=255)           gTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATAtt                                                                                                                                           >  1:427239/1‑151 (MQ=255)                acaTAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACtt                                                                                                                                      >  2:601188/1‑151 (MQ=255)                         ggTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGt                                                                                                                              >  2:1113688/1‑150 (MQ=255)                         ggTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGt                                                                                                                              >  1:791414/1‑150 (MQ=255)                               ccggccgcaaaaTAACGCAAATTTTTTGCCAGCCATTCCCACCCAG‑TCCCGCCCCCGACCCAGGCGGCGATGCTCGCCGGAAACAAAATAATACTGGAGGGGAAAATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTaa                                                                                                                            <  2:729917/137‑1 (MQ=255)                                 aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                     >  2:449054/1‑151 (MQ=255)                                  ggCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc                                                                                                                    >  1:1003555/1‑151 (MQ=255)                                  ggCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc                                                                                                                    >  1:1056636/1‑151 (MQ=255)                                   gCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                                   >  2:787405/1‑151 (MQ=255)                                      aaaaaTACAGCAAATTTGTTGCCAGCCAGTCCCAACCAG‑TCCCCCCCCCGAACCATGACGCGATGGTCGCCGGGAACAAAATAATACTGGAAGGGAAGATGGCGGGGATCACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                                  <  2:116507/147‑1 (MQ=255)                                        aaaTCCAGCAAAATTGTTGTCACCCAGTTCCAACCAG‑TCCCCCCCCCGAACCATGACGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                 <  1:96611/148‑1 (MQ=255)                                        aaaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg                                                                                                              >  2:861328/1‑151 (MQ=255)                                         aaTACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                             >  1:683294/1‑151 (MQ=255)                                         aaTACACCGAAATTGTGGCCGGCCAGTTCAACCCAG‑TCCCGCCCCCGAACATGGCCGGGGGGGCCGCGGGAAAAAAAAAAATTCGGGAAGGGAAGATGGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                             <  2:878487/151‑1 (MQ=255)                                           tACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                             >  1:654454/1‑149 (MQ=255)                                            aCAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                          >  2:1075135/1‑151 (MQ=255)                                                 gAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                        >  2:1013614/1‑148 (MQ=255)                                                  aaaTTGTTGCCAGCCATTCCCAACCAT‑TCCCCCCCCCGACCCATGACGGGTGGGTCGCCGGAAAAAAAAAAATACGGGAAGCGAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                        <  1:1075746/147‑1 (MQ=255)                                                       gttgtCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                               >  1:958009/1‑151 (MQ=255)                                                         tttccggcctttCCAACCGGTTCCCCCCCCCGAACCAGGCGGCGGTGGCCGGGGCAAAAAGATTAATACGGAAAGCGATGATGGCGGGGCTCACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATCGACacgac                                                                                               <  2:275721/140‑1 (MQ=255)                                                                 attcccAACCAT‑CCCCGCCCCCGACCCAGGCCGCGAGGGTGGCGGGAAACAAAAAAAGACTGGAAGCGAAGTGGGCGGGGCTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                         <  2:43197/143‑1 (MQ=255)                                                                     ccAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                                 >  1:807316/1‑151 (MQ=255)                                                                       aaCCAG‑TCCCGCCCCCGAACCATGACGCGATGGTCGCGGGAAAAAGAATAATAGGGGAAGCGAAGATTGCTGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACg                                                                                  <  2:395197/148‑1 (MQ=255)                                                                         ccAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                             >  2:362490/1‑151 (MQ=255)                                                                         ccAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcg                                                                             >  1:208600/1‑151 (MQ=255)                                                                           aG‑TCCCGCCCCCGACCCATGACGGGATGGTCGCGGGAAACAGAAAAATACTGAAGGCGAAGATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg                                                                                     <  1:1113996/141‑1 (MQ=255)                                                                            gggcccccccccGACCCAGGAGGGGTGGGTCGCGGGAAAAAAAAAAATACTGGAAGCGAAGTTGGCGGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgc                                                                            <  1:309886/143‑1 (MQ=255)                                                                             gccccccccccGAACCATGACGCGATGGTCGCGGGGAACAAAATAATTCTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                          <  1:769607/145‑1 (MQ=255)                                                                                   cccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                                     >  1:741329/1‑150 (MQ=255)                                                                                      aagAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTAACGTTTCGCGTAGTTTAc                                                                   >  2:863978/3‑149 (MQ=255)                                                                                        gAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                                >  2:532034/1‑150 (MQ=255)                                                                                              ggccgcgATGGTCGCGGGAAAAAAAAAAATACTGGAAGGAAAGATTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                            <  1:59089/145‑1 (MQ=255)                                                                                                 cgggATGGTCGCCGGAAAAAAAATAATATTGAAAGGAAATATTGCCGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                        <  2:1131193/146‑1 (MQ=255)                                                                                                 cgcgAGGGTCGCCGGAAAAAAAATAATACGGGAAGCGAAGATGGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                        <  1:689920/149‑1 (MQ=255)                                                                                                       ggccGCCGGAAAAAGAATAATCCTGAAAGCAAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                    <  1:150956/144‑1 (MQ=255)                                                                                                           gcgggaaaCAAAATAATACGGAAGGGGAAGTTGGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                               <  2:949753/142‑1 (MQ=255)                                                                                                           gCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGc                                             >  2:495087/1‑150 (MQ=255)                                                                                                              ggaaaaaGAATAATACGGGAAGCGAAGTTTGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTc                                           <  2:185725/143‑1 (MQ=255)                                                                                                                  aaaaaaaaagaCTGAAAACGAAGTTTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                        <  2:211094/138‑1 (MQ=255)                                                                                                                  aaaGAATAATCCGGGAAGCGAAGATCGCCGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACaaa                                     <  1:121661/149‑1 (MQ=255)                                                                                                                     gAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                  >  1:1033513/1‑151 (MQ=255)                                                                                                                     aaataatGCTGGAAGCGAAGATGGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                  <  2:1084776/150‑1 (MQ=255)                                                                                                                      aataatACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATac                                  >  2:100033/1‑150 (MQ=255)                                                                                                                         aataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                               >  2:337898/1‑150 (MQ=255)                                                                                                                               ggaagggaagATGGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                             <  1:653010/146‑1 (MQ=255)                                                                                                                                   gcaaaGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc                     <  1:470400/147‑1 (MQ=255)                                                                                                                                     gAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACt                  >  1:396459/1‑151 (MQ=255)                                                                                                                                      aataTTGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACg                       <  2:1156991/142‑1 (MQ=255)                                                                                                                                           ggCCGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa              <  2:264622/148‑1 (MQ=255)                                                                                                                                            gcgggaaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaaa             <  2:22220/143‑1 (MQ=255)                                                                                                                                             ccGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaata           <  1:288737/150‑1 (MQ=255)                                                                                                                                               gggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact         >  2:651444/1‑150 (MQ=255)                                                                                                                                               gggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact         >  2:725135/1‑150 (MQ=255)                                                                                                                                                ggATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatacta        >  2:510243/1‑150 (MQ=255)                                                                                                                                                     aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTg  >  1:443327/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAACCAG‑TACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACTG  >  NC_000913/3443130‑3443427

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