Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,601,160	0	G	A	43.9%	63.8 / 31.7	57	F309F (TTC→TTT)	ftsX	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (27/5);  new base A (0/25);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 2.32e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATC  >  NC_000913/3601011‑3601307                                                                                                                                                      |                                                                                                                                                     agTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCg                                                                                                                                                      >  1:736429/1‑150 (MQ=255) agTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCg                                                                                                                                                      >  1:382174/1‑150 (MQ=255)      aaCACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑cc                                                                                                                                               >  1:658203/1‑151 (MQ=255)       acacGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCt                                                                                                                                              <  1:577387/151‑1 (MQ=255)         acGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCtgt                                                                                                                                            >  1:816376/1‑151 (MQ=255)          cGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCtgtg                                                                                                                                           >  2:792645/1‑151 (MQ=255)               ttattaTTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACt                                                                                                                                      >  2:852761/1‑151 (MQ=255)                  ttattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTc                                                                                                                                    <  2:560218/150‑1 (MQ=255)                    attCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGc                                                                                                                                 >  1:33579/1‑151 (MQ=255)                    attCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGc                                                                                                                                 <  2:429896/151‑1 (MQ=255)                         ggCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAAcc                                                                                                                             <  2:1036408/150‑1 (MQ=255)                         ggCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAAcc                                                                                                                             >  1:584479/1‑150 (MQ=255)                                aGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATg                                                                                                                      >  1:343292/1‑150 (MQ=255)                                    gCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAAt                                                                                                                 >  1:840113/1‑151 (MQ=255)                                        aaaTGTTGTCCCGGGGCAACCCCCGCTCCCCCCCAGCCAATCATCGAGCATACCAGCAGAAATAGCAGCCATTCAACAAATGATAACCCATTAATATAAACCTTCGTCCCAAAAC‑CCTGTCCCCCTCCCGCACCCGCCGATGACAATCgc                                                                                                              <  1:250380/150‑1 (MQ=255)                                         aaTGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCgc                                                                                                              >  1:233939/1‑149 (MQ=255)                                             ttggccggTGGCAACCCCCCCGCCCCCCCAACCAATCATCGACCACCCCAGCAGCAATACCAGGCATTCATCGAATAAAACGCAATTAATAAAAACCTTCGTACCAAAAA‑CCGGTCCCCCTCCCGAAACCGCCGATGACAATCgcagca                                                                                                          <  2:655928/142‑1 (MQ=255)                                                tACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCAcc                                                                                                        >  2:645955/1‑148 (MQ=255)                                                tACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCAcc                                                                                                        >  2:226167/1‑148 (MQ=255)                                                    gTGGAAACCCACGCGCCCCCCCACCAAATCACCGAGCATCCCAGAAGAAATAGCAGGCATTCATCAAATGAAAAGCCATTGATAAAAAACTTCGTCCCAAAAC‑CCTGTGCCACTCCCACAACCGCCGATGACAATCGCAGCACCAg                                                                                                      <  1:674552/146‑1 (MQ=255)                                                         aaGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGAc                                                                                            >  1:1007179/1‑151 (MQ=255)                                                            ccACGCTCCCCCCCAGCAAATCACCGACCATCCCAGAAAAAAAAGCAGGCATTCATCAAATGAAAAGCCATTAATATAAACCTTCGTCCCAAAAC‑CCTGTCCCACTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaa                                                                                          <  1:901563/150‑1 (MQ=255)                                                                 cccccccccacccaATCACCGACCAACCCAGCAGCAATAGCAGGCATTAACCAAATGATAAGCAATTGATATCAACCTCCTTCCCAAACC‑CCTGTCCCCTTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAt                                                                                       <  1:487611/142‑1 (MQ=255)                                                                 cTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTa                                                                                     >  1:68932/1‑150 (MQ=255)                                                                    ccccccAGCCAATCATCGACCATCCCAGCAACAAAGGCAGGCATTCACCGAATAAAAAGCCATTAAAATCAACTTTCTTCCCAAAAA‑CCTGTGCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAAt                                                                                   <  1:249740/146‑1 (MQ=255)                                                                       cccAGCAAATCAACGAGCATACCAGAAAGAAGAACGGGCATTCATAGAAGGATAAGCATTTGATACAAATCTCCGTCCCAAAAG‑CCGGGCCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGaca                                                                               <  1:613177/150‑1 (MQ=255)                                                                       cccAACAAATCACCGAGCATCCCAGAAGAAATAGCAGGCTTTCATCGAATGATAAGCATTTAATATAAAACTTCGTTCCAAAAC‑CCGGTGCCCTTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGac                                                                                <  1:305910/149‑1 (MQ=255)                                                                       ccaaacaaaaaatcGAGCATACCAGCAACAAGCGCGGGCATTCATCAAATAAAAAGCATTTGATATAAAACTTCGTCCCAAAAC‑CCTGTCCCTCCTCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaa                                                                              <  2:313526/140‑1 (MQ=255)                                                                              aaTCATCAACCATCCCAGCAGAAATAGCGGGCATTCATCAAATAATAAGCTATTAATAAAAATCTTCGTCCCAAAAC‑CCTGTGCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg                                                                         <  1:183656/149‑1 (MQ=255)                                                                                 catcGAGCATACCAGCAGAAAAAGCAGGCATTCATAAAATAAAAAGCCATTAATAAAAACCTCCGTCCCAAAAC‑CCTGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgcg                                                                       <  2:879077/148‑1 (MQ=255)                                                                                 caccGACCACACCAGCAACAATAGCAGGTTTTCATCAAATGATAAGCCTTTAATATAAATCTTCGTCCCAAAAC‑CCTCGCCCCCTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATg                                                                         <  1:849206/143‑1 (MQ=255)                                                                                         atcccAGCAGAAATAGCAGGCATTCATCGAATAAAAAGCTATTGATATCAACCTTCGTCCCAAAAC‑CCTGTGCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgc                                                                        <  1:45768/136‑1 (MQ=255)                                                                                            ccAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCTTTGATATAAACCTTCGTCCCAAAAC‑CCTGTCCCTCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATccc                                                           <  1:588243/149‑1 (MQ=255)                                                                                              agcagaAATAGCAGGCTTTCATGAAATAAAAACCAATTTATATAAACCTTCTTCCCAAAAC‑CCTGTGCCACTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaaaa                                                               <  1:734121/143‑1 (MQ=255)                                                                                                    aaaagcagcCATTCAAAGAATAATAAGCTATTAATATCAACCTTC‑TTCCAAAAACCCTGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTgc                                                   <  1:895661/146‑1 (MQ=255)                                                                                                    aaTAGCGGGCATTCATCAAATAATAAGCCATTAATATAAACCTTCGTCCCAAAAC‑CCTGTCCCCCTCCCAAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTg                                                    <  1:51117/148‑1 (MQ=255)                                                                                                      tAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGcc                                                >  1:690981/1‑150 (MQ=255)                                                                                                        gCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGAAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGccacc                                             >  2:633354/1‑151 (MQ=255)                                                                                                           ggaaTTCATCAAATAAAAACCCATTGAAATCAACCTTCTTTCAAAAAC‑CCTGTCCCCCTCCCGAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGcca                                               <  2:829643/143‑1 (MQ=255)                                                                                                              aTTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATAc                                         >  1:791923/1‑149 (MQ=255)                                                                                                               ttCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAgg                                      >  1:391741/1‑151 (MQ=255)                                                                                                                  aTCAAATGAAAACCCTTTGATATAAACCTTCTTCCCAAAAC‑CCTGTCCCACTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGa                                     <  1:824871/149‑1 (MQ=255)                                                                                                                          aTAACCCATTTATATAAACCTTCGTCCCAAAAC‑CCTGTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAg                           <  1:500607/151‑1 (MQ=255)                                                                                                                            aaGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAgg                          >  1:800312/1‑150 (MQ=255)                                                                                                                              gCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATg                       >  2:1066848/1‑151 (MQ=255)                                                                                                                               ctaTTAATATCAAACTTCGTCCCGAAAC‑CCGGCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                        <  1:513101/147‑1 (MQ=255)                                                                                                                               ccATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATg                       >  1:909162/1‑150 (MQ=255)                                                                                                                                  ttGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc                   >  2:499557/1‑151 (MQ=255)                                                                                                                                   taatatAAAACTTC‑TTCCAAAAACCCTGTCCCCCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATc                   <  1:429904/148‑1 (MQ=255)                                                                                                                                   tGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATcc                  >  1:288160/1‑151 (MQ=255)                                                                                                                                    gATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATcc                  >  1:972970/1‑150 (MQ=255)                                                                                                                                     atatCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCa                 >  1:601232/1‑150 (MQ=255)                                                                                                                                      tatCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGA‑AATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATc               >  2:315017/1‑150 (MQ=255)                                                                                                                                        acAACCTTCGTACCAAAAC‑CCTGCCCCCCTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCt              <  1:858177/149‑1 (MQ=255)                                                                                                                                         cAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGt            >  2:332493/1‑151 (MQ=255)                                                                                                                                          aaaCTGCGTCCCAAAAC‑CCTGTCCCCTTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTc           <  2:825544/151‑1 (MQ=255)                                                                                                                                           accTTCGTCCCAAAAC‑CCTGTGCCACTCCCCCAACCGCCGATGACAATCGCAGCACCTGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCg          <  1:676544/149‑1 (MQ=255)                                                                                                                                                  ttCCAAAAC‑CCTGTGCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCa        <  1:11515/146‑1 (MQ=255)                                                                                                                                                    ccGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATc  >  2:582461/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGTATAACACGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CCTGTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATC  >  NC_000913/3601011‑3601307

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