Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,516	C→A	33.7%	Q196K (CAA→AAA)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,516	0	C	A	33.7%	62.3 / 51.1	65	Q196K (CAA→AAA)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (17/26);  new base A (6/16);  total (23/42)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

GAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAT  >  USA300TCH1516_ALE/111378‑111652                                                                                                                                           |                                                                                                                                         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)                                                                                                                                        gCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAt  <  2:449891/140‑1 (MQ=255)                                                                                                                                        gCTCAATATCAAATAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAt  <  1:473880/140‑1 (MQ=255)                                                                                                                                        gCTCAAGATCAATGAACTAATGATGATTATAATGTAAAACGATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAt  <  2:337230/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                         GAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAAT  >  USA300TCH1516_ALE/111378‑111652

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