Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	380,012	2	.	G	69.8%	104.6 / 47.7	43	intergenic (‑131/+57)	frmR/yaiO	regulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (9/4);  new base G (16/14);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCA  >  NC_000913/379865‑380154                                                                                                                                                         |                                                                                                                                                    tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaa‑‑a                                                                                                                                                       <  2:1198640/148‑1 (MQ=255) tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑ggg                                                                                                                                                 >  2:221028/1‑151 (MQ=255) tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gg                                                                                                                                                  <  1:1059876/150‑1 (MQ=255)  ccGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑ggg                                                                                                                                                 <  2:913777/150‑1 (MQ=255)  ccGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑ggg                                                                                                                                                 <  2:374757/150‑1 (MQ=255)   cGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑g                                                                                                                                                   <  2:347828/147‑1 (MQ=255)            gggCGTTTCGCACCCCTTCACGTGGATGTAAGATTCTAGAGGGGGGTATTCTATAT‑GTCAATGCATCCCCCCCTATAGTACAATACACGGGAGAATAAAATTTATCCGGTGAATGCGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gggggg‑‑ggggg                                                                                                                                       <  1:1050290/148‑2 (MQ=255)             ggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gggggg‑‑ggggggct                                                                                                                                    >  2:723633/1‑146 (MQ=255)                 tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCgg                                                                                                                              >  2:469283/1‑151 (MQ=255)                 tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gggggg‑‑ggggggctaatcg                                                                                                                               >  2:972538/1‑142 (MQ=255)                 tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gggggg‑‑ggggggctaatcg                                                                                                                               >  1:577816/1‑142 (MQ=255)                              cctcTTCATTCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGgcc                                                                                                                   >  2:877786/3‑151 (MQ=255)                               atcTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGgccgc                                                                                                                 >  1:856937/1‑151 (MQ=255)                               atcTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAAAGTGGTCGGAAAACAAAGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAAAGGGAGGGGAAGgccccc                                                                                                                >  2:821124/1‑151 (MQ=255)                                      atgaattaacctaaagggggTATTCT‑TATAGTTAAATCCTACCCCCCCTAAGTAAAAAACCCGGGGAAAAAAAATTTTTCCGGGGGAAGGGGGCGGAAAAAAAAAGGGAA‑‑AAG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                          <  1:336024/136‑1 (MQ=255)                                         aaaatccTAAAGGGGGGTTTTTTAATT‑TTAAAACCATCCCCCCCTATAGTTAAAAACCGTGGAAAAAAAAATTTATCCGGGGAATGTGGTCGGAAAAAAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                          <  2:20392/142‑1 (MQ=255)                                          atatACTATAGGGGGGTATTATAAAT‑GTCAATGCATCCCCCCCTATAGTATAAAACACAAGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGAAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgg                                                                                                     <  2:321147/151‑1 (MQ=255)                                            ataCTATAGGGGGGTATTCTATAT‑GTCAAGGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATGGGCAGGGAAGGACGCCCCGGATAgcgg                                                                                                     >  1:823374/1‑149 (MQ=255)                                                 ataGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc                                                                                                >  2:870411/1‑150 (MQ=255)                                                 ataGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGAAGGGAAGGCCGCCCCGGATAGCGggcggc                                                                                                >  2:688760/1‑150 (MQ=255)                                                    ggggtgTAGTCTTTAT‑TTAAAGGCTTCCCCCCTTTTAATACAAAAACCAGGAGAAAAAAATTTATCCGGCAAAGGGGCTCGGAAAACAAAGGGAAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc                                                                                                   <  1:1033070/144‑1 (MQ=255)                                                    ggggggTTTTCTATAT‑TTCAATGCATACCCCCCTATAGTATAAAACACTGGAAAAAAAAATTTATCCGGTGAATGTGGGCGGAAAAAAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAg                                                                                              <  2:50050/149‑1 (MQ=255)                                                    ggggggTTTTCAAATT‑GTCAATGAATACCCCCCTATTATATAAAAACGTGGAAAAAAAAATTTTTCCGGGGAATGGGGGCGGAAAAAAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                                             <  2:646925/150‑1 (MQ=255)                                                     gggggTATTCTATAG‑GGCAATACATCCCCCCCTATAATAAATAACACTGGAGAAAAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc                                                                                                   <  1:233016/143‑1 (MQ=255)                                                           aTTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑GGGGGG‑‑GGGGGCCTATTCGCGAGGGAAGGCCGCCCCGGATGGCGGGCGGCAGAGGGAAt                                                                                       >  2:1104775/1‑149 (MQ=255)                                                             tCTATAT‑GTCAATGCAGACCCCCCTATAGTAAATAAAACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAAGGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAg                                                                                    <  1:326440/150‑1 (MQ=255)                                                             tCTATAT‑GACAATGCAACCCCCCCTAAATTATAAAACAGGGAAAAAAAAAATTTAGCGGGGGAATGTGGGGGGAAAAAAAAGGGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATc                                                                                      <  2:725654/148‑1 (MQ=255)                                                                 tat‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                                >  2:536718/1‑150 (MQ=255)                                                                  at‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑GGGGGG‑‑GGGGGCCTAACCGGCGGGGAAGGCCGCCCCGGGTACGGGGCCGAAGAAGGAGTCAGAttt                                                                                >  1:1084852/1‑149 (MQ=255)                                                                     gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                             >  2:1158366/1‑150 (MQ=255)                                                                            cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGAc                                                                    >  2:900501/1‑151 (MQ=255)                                                                            cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGT‑GTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg                                                                      >  2:359745/1‑149 (MQ=255)                                                                             aaaCCCCCTTATATTATTAAACACTGGAGAAAAAAATTTATCCGGGGAATGGGGTGGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCa                                                                       <  1:869114/146‑1 (MQ=255)                                                                                          gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt                                                        >  2:365614/1‑150 (MQ=255)                                                                                             tataACACGGGAGAATAAAATTTATCCGGTTAATGTGGTCGGAAAACAAGGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt                                                        <  1:972000/147‑1 (MQ=255)                                                                                                            tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGGGGGCGGCAGAAGGAATCAGAATTTCCAGGACAGACGGGCTGCAAGTTGCAGACGGTTAAAatcat                                      >  2:1087632/1‑150 (MQ=255)                                                                                                             aaaaTTTACCCGGTAAATGGGGTCGAAAAACAAAGAAAAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcatc                                     <  2:1163006/150‑1 (MQ=255)                                                                                                              aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCAGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCgg                                   >  1:968885/1‑151 (MQ=255)                                                                                                               aaTTTATCCGGTGAATGTGGTCGGAAAACAACGTGGAA‑‑A‑‑‑GGGGGG‑‑GGAC‑‑ATAATCGGCCGGGAAGGCGTGCATGGATAGCGGCTGGCCGAGGGGCAGAAAATTTCGAGTTCTATCGGGGTGAAAGTTGCGGCGCGTTAAAATCATCGGTTg                                >  1:654595/1‑151 (MQ=255)                                                                                                                   tATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTg                                >  2:780337/1‑149 (MQ=255)                                                                                                                     tCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATAGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg                             >  1:889049/1‑150 (MQ=255)                                                                                                                      ccGGTGAATGGGGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgt                            <  1:59199/150‑1 (MQ=255)                                                                                                                         gTGAGTGTGGTGGAAAAAAAAAGAGGAAGGG‑‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtgt                          <  2:141403/149‑1 (MQ=255)                                                                                                                          tGAATGGGGTCGAAAAAAAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCg                        <  2:751247/150‑1 (MQ=255)                                                                                                                               gtgGTCGGAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAcc                    >  2:443560/1‑149 (MQ=255)                                                                                                                                     ggAAAACAAAGAGGCA‑‑AGG‑GGGGGG‑‑GGGG‑‑CAAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGCATTTCCAGGTAAGACGGGCTGCAAGTGTCAGACCGTTAAAATCTTCGGTTGGCGTGTCGTACCACACttt             >  1:650826/1‑150 (MQ=255)                                                                                                                                      gAAAACAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAc           >  1:97617/1‑151 (MQ=255)                                                                                                                                       aaaaCAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTa            <  1:916267/149‑1 (MQ=255)                                                                                                                                         aaCAAAGAGGAA‑‑AGG‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGGGGGCGGCAGAAGGAATAAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt         >  2:796670/1‑150 (MQ=255)                                                                                                                                            aaaGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCg     >  1:276459/1‑150 (MQ=255)                                                                                                                                             aaGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc   <  2:763111/151‑1 (MQ=255)                                                                                                                                             aaGAGGAA‑‑A‑‑‑GGGGGGGCGGGG‑‑CTAAACGGCAGGGAAGGCCCCCCCGGATAGCGGGCGGGAGAAGGAACCAGAATTTCCAGGTCAGCCGGGCCGCAAGTTGCAGACCGTTAAAATCATCGGGTGGGGTGCCGCACCCCCCTTTACCTGcc      >  1:32678/1‑149 (MQ=255)                                                                                                                                               gagGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCa  <  1:469269/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                    TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATAT‑GTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCA  >  NC_000913/379865‑380154

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