Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,206,435	0	C	T	44.1%	56.9 / 36.9	59	intergenic (‑179/‑28)	ttdR/ttdA	transcriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (1/32);  new base T (26/0);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 6.82e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAG  >  NC_000913/3206291‑3206579                                                                                                                                                 |                                                                                                                                                  aTTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCaaa                                                                                                                                               >  2:118073/1‑147 (MQ=255)    tCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACaaa                                                                                                                                           <  2:972123/150‑1 (MQ=255)      gCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCAGACaaaaa                                                                                                                                         >  1:957732/1‑148 (MQ=255)       cAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCaaaaaa                                                                                                                                            >  2:292536/1‑141 (MQ=255)          aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaaa                                                                                                                                    <  2:48871/151‑1 (MQ=255)          aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaaa                                                                                                                                    <  1:1138120/151‑1 (MQ=255)               tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATAaaaaaa                                                                                                                                   >  1:83690/1‑147 (MQ=255)               tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATgg                                                                                                                               >  1:217614/1‑151 (MQ=255)                gCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGa                                                                                                                              <  2:875638/151‑1 (MQ=255)                    ttCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAAAAATAAAAAAACTgggg                                                                                                                             >  2:1059340/1‑146 (MQ=255)                     tCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATgg                                                                                                                               >  1:1046476/1‑145 (MQ=255)                      cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGc                                                                                                                         <  1:289817/150‑1 (MQ=255)                       ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCa                                                                                                                       <  2:1110350/151‑1 (MQ=255)                         tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATTGAggttgcat                                                                                                                      >  1:212018/1‑143 (MQ=255)                          gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATAAAAAAAATGGAGGTGGCa                                                                                                                       >  1:621013/1‑148 (MQ=255)                          gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatga                                                                                                                    <  1:684190/151‑1 (MQ=255)                             ctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatgat                                                                                                                   <  1:343087/149‑1 (MQ=255)                             ctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatgat                                                                                                                   <  2:70742/149‑1 (MQ=255)                                  gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                                                                                            <  2:509778/151‑1 (MQ=255)                                          cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATTGAGGTGCCCTTATTAGCGGAAGtaata                                                                                                     >  1:1102430/1‑150 (MQ=255)                                           attcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaata                                                                                                     <  2:61489/149‑1 (MQ=255)                                             tcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATAAAAAAAAAGGGGGTTGCAAGATGAGTGAAAAtaaaaa                                                                                                    >  2:293616/1‑148 (MQ=255)                                                 aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAATTTTTTCAAACAAATAAAAAAAATGGGGGTGCCATGATGAGCGAAAATAATAAG‑Caa                                                                                               >  1:1019430/1‑148 (MQ=255)                                                  ttttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAAAAATACAAAAACTGGAGTTGCCATGTTTTGCGAAAATAATAAA‑Caaaa                                                                                             >  2:88284/1‑147 (MQ=255)                                                    ttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGc                                                                                          <  2:996463/150‑1 (MQ=255)                                                     tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCa                                                                                         <  1:129727/150‑1 (MQ=255)                                                      cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAAAAAATACAAAAAATGGAGTTGGCATGATTAGCGAAAGTAATAAa‑caaca                                                                                             >  2:252956/1‑145 (MQ=255)                                                        cgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGa                                                                                     <  1:891949/151‑1 (MQ=255)                                                                cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTTTATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGt                                                                               <  1:868808/149‑1 (MQ=255)                                                                 aCGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGGCATGATGGGCGGAAGTAAAAAA‑AAAAAGGCAGGGAAAAAAtt                                                                              >  1:430205/1‑149 (MQ=255)                                                                    aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACag                                                                         <  1:67667/151‑1 (MQ=255)                                                                       tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTCCACAAAAACAAAAACTGGAGGTGCCATGATGAGCGAAAAAAAAAAA‑CAACAGGCAGTGGATAAGTTGGCa                                                                          >  1:25576/1‑147 (MQ=255)                                                                           aaTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCTATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTg                                                                   <  2:662448/150‑1 (MQ=255)                                                                               aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCt                                                              <  1:640848/151‑1 (MQ=255)                                                                                cATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTa                                                             <  1:729013/151‑1 (MQ=255)                                                                                 atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCACACAAAAACAAAAAAAGGAGGTGTCATGATGAGGGAAAGTAATAAA‑CAAAAGGCAGTGAGAAAAATGACAGAGAATGt                                                                  >  2:450156/1‑145 (MQ=255)                                                                                      cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTa                                                       <  2:982585/151‑1 (MQ=255)                                                                                        cGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTAcc                                                     <  1:854894/151‑1 (MQ=255)                                                                                          gTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACTGGGGTTGCCATGATGAGCGGAAGTAATAAG‑CAACAAGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGc                                                   >  1:112470/1‑151 (MQ=255)                                                                                          gTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTAAACAAATAAAAAAAAATGGGTTGCCATGATGAGGGAAAATAATAAA‑CAAAAAGGAGAGAATAAAGTGACAGAGAATTTTGCTAAATTTaac                                                     >  2:100658/1‑147 (MQ=255)                                                                                             gCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAACCAAATAAAAAAACCAGAGTTGCCCTTATGAACCAAAGGAAATAA‑AAAAAAGCAGTTAATAAATTGACAGAGAATGTCGCCAAAtttttctccat                                                >  1:886958/1‑143 (MQ=255)                                                                                               aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTAAACAAAAAAAAAAACTGGGGGTGGCATGATGAGCGGAAGTAATAAA‑AAACAAGCAGTGGAAAAATTGACAGAGATTTTAGCTAACTTTACCGCCaag                                               >  2:172512/1‑148 (MQ=255)                                                                                               aaaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAATTTTTTCACACAAATAAAAAAAATGGAGTTGCCATGATGTGCGAAAGTAATAAA‑CAACAAGCAGTGAATAAGTTGACAGAGATTTTCGCTAAATTTTCCGccc                                                 >  2:444903/1‑147 (MQ=255)                                                                                                aaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCCGACAAAAACAAAAACTGGAGGTGCCCTGATGAGCGAAAAGAAAAAA‑AAACAAGCAGTGAAAAAGATGACAGAGAATGTCGCTAAATTTAACGCCATGAt                                             >  1:945890/1‑151 (MQ=255)                                                                                                       taAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCTGACAAATAAAAAAAATGGAGGTGCCATGATGAGCGAAAGAAATAAA‑AAAAAGGCAGAGGAAAAAGTGACAGAGATTGTCGCTAACTTTACCGCCCTGATTTTTAcc                                      >  1:190582/1‑151 (MQ=255)                                                                                                        aaaaCCACATTTTTAGTGGTAGTTGGTGGCGGTGAATTTTTCCAGACAAAAACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTAcc                                      <  1:925187/150‑1 (MQ=255)                                                                                                            ccACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAAt                                 <  1:738590/151‑1 (MQ=255)                                                                                                             cacaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAAAGGAGGTGGCCTGATGAGGGGAAAAAAAAAG‑AAAAAGGCAGTGAAAAAATTGACAAAGATTGTTGCTAACTTTAACCCCATGATTTTTATCCGaa                                  >  1:1022120/1‑149 (MQ=255)                                                                                                                aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAACCGGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCt                             <  1:956506/151‑1 (MQ=255)                                                                                                                  ttttAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATGGAGTTGCCATGATGAGCGGAAGTAATAAA‑CAAAAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTTTACCCGAATGcc                              >  1:755699/1‑148 (MQ=255)                                                                                                                  ttttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtg                            <  1:599543/150‑1 (MQ=255)                                                                                                                   tttAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACAGGGGGTGCCATGATGAGCGGAAGAAATAAA‑CAAAAGGCAGGGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCCAGATTTTTTCCCCAATGCCtga                           >  1:554599/1‑150 (MQ=255)                                                                                                                       gTGGTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatg                         <  2:550845/148‑1 (MQ=255)                                                                                                                        tGGTAGTTTGTGGCGTTGAATTTTCCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAC‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACg                      <  2:1110414/150‑1 (MQ=255)                                                                                                                        tGGTAGTTTGTGGCGGTAAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑AAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGAc                       <  1:411692/149‑1 (MQ=255)                                                                                                                         ggTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgt                     <  1:795140/150‑1 (MQ=255)                                                                                                                                      gtggAATTTTCCCAGACAAAGAAAAAAACTGGAGTGCCAATGAGGAGCGAAA‑TAAAAAGCCACCAGTATGTGAATAACTCGACAGAGTTTTTCGCTACCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTa        <  2:1118127/147‑1 (MQ=255)                                                                                                                                           aTTTTTCCAGACAAATACAAAAACTGGAGTTGCCGTGAAGACAGAAAATAATAAA‑CAACAGGCAGTGAATAAGTTGGCAGAGATTGTCGCTAACTTTACCGCCATGATTTCTCCCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAg  <  1:130561/151‑1 (MQ=255)                                                                                                                                            tttttCAAAAAAAAAACAAAAACTGGATTTGCCATAATGAGCGAAAGAAATAAC‑AAACAGCAAGTAAATAATTGGAAAGAGTTTTTCCCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAg  <  2:1102286/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  ATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAG‑CAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAG  >  NC_000913/3206291‑3206579

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