Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,512	A→G	35.5%	S194S (TCA→TCG)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,512	0	A	G	35.5%	51.7 / 57.1	65	S194S (TCA→TCG)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (18/24);  new base G (6/17);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

TACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGA  >  USA300TCH1516_ALE/111374‑111641                                                                                                                                           |                                                                                                                                  tACTGAGGAAGTCAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAg                                                                                                                                  <  2:396637/140‑1 (MQ=255)  aCTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAATGGTGCCGAGTTATTAAga                                                                                                                                 >  1:265997/1‑139 (MQ=255)  aCTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGAAATATCATATAATCCAGTGGTGCTGAGTTATTCGGc                                                                                                                                 <  2:459697/140‑1 (MQ=255)     gAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCa                                                                                                                              <  2:537279/140‑1 (MQ=255)      aGGAACTGAAAGAAAAAAAAAAAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCaa                                                                                                                             <  2:36342/140‑1 (MQ=21)           gTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATGTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCa                                                                                                                        <  1:70591/140‑1 (MQ=255)              aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                                                                     <  1:248645/140‑1 (MQ=255)               aagaagAAAACCTAAAAAAAGAAATCGAAAACGTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                                                    <  2:309270/140‑1 (MQ=18)               aagaagAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATAAATTa                                                                                                                    >  2:385277/1‑140 (MQ=255)                 gaataaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                                                  <  1:518097/140‑1 (MQ=255)                 gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAAc                                                                                                                  >  1:529735/1‑140 (MQ=255)                  aagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTTCTCTATATTTAAATATTAAAGATTTAATGAATTATAAAGATGGAGATATAGCATATAATCCAGAAGTCTCGAGTTATTCGGCTAAAGATCAATTAACt                                                                                                                 >  2:204320/1‑140 (MQ=14)                  aagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                 <  1:370919/140‑1 (MQ=255)                    gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                                                               <  1:503062/140‑1 (MQ=255)                    gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTaa                                                                                                               >  2:468747/1‑140 (MQ=255)                     aaaaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                                              <  2:326125/140‑1 (MQ=255)                       aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGATATATATCATATAATCAAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatga                                                                                                            >  1:43450/1‑140 (MQ=255)                        aaaTAAAAAAAGAAATCGAAGACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                                                           <  2:516046/138‑1 (MQ=255)                        aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTCTGCGTAAAAGGGAACTTTTGAAGATTTAATGCATTATAAAGATGGAGATATATCATATAATCCAGAGGTACCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                           >  1:265024/1‑140 (MQ=14)                         caaaaaaaaaGAAAACGAAAAATTTAAGGTTCTTGCGCAAAAGAGCAATTTCAAAGATTTAAGGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCAGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                          <  2:242801/137‑1 (MQ=18)                         cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgatg                                                                                                          <  1:60295/140‑1 (MQ=255)                         cATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTAAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgctg                                                                                                          >  2:369079/1‑140 (MQ=255)                          aaaaaaaaaGAAAGCGAAAAATTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGTTCTATCTTATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAtgatga                                                                                                         <  1:290023/138‑1 (MQ=21)                          aTAAAAAAAGAAATAGTAAACTCTAAGTGTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATATAGATGGAGATATATCATATAAGCCAGAGGTGCCGACTTATTCAGCTCAATATCAATTAACTAAtgatga                                                                                                         <  2:158909/140‑1 (MQ=18)                           tAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCCATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGATATTCGGCTAAATATCAATTAACTAATgatgat                                                                                                        >  1:99619/1‑140 (MQ=255)                              aaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                     <  1:387081/140‑1 (MQ=255)                                  ggAATTCTAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATAACAATTAACTAATGATGATTATAATg                                                                                                 >  2:142812/2‑140 (MQ=18)                                  acaaaTCTATAACTTTAAGTTTTTTGCACAATATGGCTATTTTAATCAGTTAATGATTTTTATAGTTGGAGATATATCATAGAATCCAGAGGTGCCGAGTTAGTCGGCTAAATATCAATGAACTAATGATGATTATAATg                                                                                                 <  1:195674/138‑1 (MQ=14)                                  aGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATg                                                                                                 <  1:167452/140‑1 (MQ=255)                                    aattCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTa                                                                                               <  2:183423/137‑1 (MQ=255)                                      aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTTAAGATTTAATGAATTATAAAGATTGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                             <  2:518185/140‑1 (MQ=255)                                      aTCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                                                             >  2:291753/1‑140 (MQ=255)                                       tCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaaaa                                                                                            <  1:506209/140‑1 (MQ=255)                                            aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTACTGAATTATAAAGATGGAGATATATCATATAATACAGAGGTTCCGAGTTATTCAGCCCAATATCAATTAACTAATAATGTTTATAATGTAAAACAAtt                                                                                       >  1:199829/1‑140 (MQ=255)                                              cTTTATGTTTTTTGCGCAATATTGGTATTTTTATGATTTTATGATTTATAAAGATGGAGATATATCATATAGTCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                                                     <  1:252393/140‑1 (MQ=14)                                              cTTTAAGTTTTTTGCCCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                                                                     <  1:29332/140‑1 (MQ=255)                                                 tAAGTTTTTTGCGCAATTTTGGAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgca                                                                                  <  2:495868/140‑3 (MQ=17)                                                 tAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgca                                                                                  <  2:388604/140‑3 (MQ=255)                                                          ttcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATg                                                                         <  1:447401/138‑1 (MQ=255)                                                          tGCGCAATTGAGCAATTTTTAAACTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATg                                                                         <  1:370207/140‑1 (MQ=14)                                                          tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAATACAATTACGCAAAAGATATg                                                                         <  1:334040/140‑1 (MQ=17)                                                          tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATg                                                                         <  1:278094/140‑1 (MQ=255)                                                          ggcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATg                                                                         <  2:60147/139‑1 (MQ=255)                                                                    aGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGTTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaac                                                               <  1:348031/140‑1 (MQ=255)                                                                              aaGATTTAATGAATCATAAGGATGGAGATATATCATATAATCCACAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg                                                     >  1:302381/1‑140 (MQ=255)                                                                               aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAAGAATTACGTAAAAGATATGACATACCAACCAATAAAgcgc                                                    >  2:160426/1‑140 (MQ=255)                                                                               aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcgc                                                    >  1:449891/1‑140 (MQ=255)                                                                                 aTTAAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg                                                  <  2:481761/140‑1 (MQ=255)                                                                                  tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGa                                                 <  2:155972/140‑1 (MQ=255)                                                                                  tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGa                                                 >  1:284125/1‑140 (MQ=255)                                                                                       tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCGATTAACTAATGAAGATTATAATGTAAAACAATTACGAAAAAGATATGACATAACAAGGAAAAAAACGCCGAAgaag                                            >  2:408614/1‑137 (MQ=17)                                                                                       tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                            <  1:444282/140‑1 (MQ=255)                                                                                           ttATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttgt                                        <  1:535356/140‑1 (MQ=255)                                                                                             ataAAGATTGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGttgttg                                      <  1:84710/140‑1 (MQ=255)                                                                                                   aTGGAGATATATCATATAATCCAGAGGTGCAGAGTTAGTCAGCTCAATATCAATTAACTAATGATGATTATAATTTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGt                                <  2:265997/140‑1 (MQ=255)                                                                                                       agaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                            >  2:317976/1‑140 (MQ=255)                                                                                                        gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                           >  1:78145/1‑140 (MQ=255)                                                                                                        gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGAGCCGAAGCTGTTGTTGAAAGGTAAAga                           >  1:78121/1‑139 (MQ=255)                                                                                                                atatAATCCATATGTGTCGAGTTATTCAGCTCAATATTTATTAACTTATGACGATTATAATGTAAAACAATTACGTAAAAGATAGTACATACCAACAAATAAAGCGCTGAAGCTGTTGTTGAAAGGTACAGGTAATTTaa                   <  2:75311/140‑1 (MQ=14)                                                                                                                  ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTGGAAAGGTACAGGGAATCTaaaa                 >  2:131860/1‑140 (MQ=255)                                                                                                                  ataATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGATGTTGTTGAAAGGTACAGGGAATTTaaaa                 >  2:181528/1‑140 (MQ=255)                                                                                                                    aaTTCAGAGGTGCCGAGTTATTCAGCTCAATATCCATTAATTAATGAGGATTATAATGTAAAATGATTACGTAATAGATATGACATTCCAACTAATAGATCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAgg               <  2:429571/140‑1 (MQ=14)                                                                                                                     aTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGt              >  2:13791/1‑140 (MQ=255)                                                                                                                        cAGAGGTGCCGAGTTATGCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtca           <  1:82478/140‑1 (MQ=255)                                                                                                                                 cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGGAAACCAACTACGTAAACGATATGACATACCAACAAATAAAGCGCCGAAGTTGTTGTTGAAAGGTGCAGGGAATTTGAAAGGACACTCAGTTGGa  >  2:323044/1‑140 (MQ=14)                                                                                                                                           |                                                                                                                                  TACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGA  >  USA300TCH1516_ALE/111374‑111641

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