Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,431	0	A	C	41.8%	84.7 / 17.0	56	G34G (GGT→GGG)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (24/8);  new base C (0/23);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 4.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAT  >  NC_000913/2248287‑2248581                                                                                                                                                 |                                                                                                                                                        aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGt                                                                                                                                                    >  1:921367/1‑150 (MQ=255)     gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑g                                                                                                                                               >  1:1034616/1‑150 (MQ=255)         cGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATgg                                                                                                                                           >  1:226976/1‑150 (MQ=255)           gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAAc                                                                                                                                        <  2:580491/151‑1 (MQ=255)           gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAAc                                                                                                                                        <  2:716254/151‑1 (MQ=255)            aaaCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAAcc                                                                                                                                       <  1:410741/151‑1 (MQ=255)              aCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGc                                                                                                                                     <  2:660900/151‑1 (MQ=255)                 ggCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCa                                                                                                                                   >  1:935662/1‑150 (MQ=255)                 ggCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCa                                                                                                                                   >  1:629919/1‑150 (MQ=255)                              ccAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGt                                                                                                                     >  1:1038296/1‑151 (MQ=255)                                atttCCCCGAGCCTGGTCATCAGAAAGAAACCCACCGGCCCAACAAGAATATACGAGAGAAATCCCCCCCCGGCCCCTCCCTGAAAAATCTTTCCGCCAGAGCAAAAAAAAAGCCCTGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCa                                                                                                                   <  1:891001/149‑1 (MQ=255)                                  ttCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAgg                                                                                                                 >  2:871081/1‑151 (MQ=255)                                    cACCGAGACTGGTCATAAGAAAGTAACCCATCAGGCCAATCAGCATATACGAGAGAAATCCCCCCCCCGGACCTGCCTGAAAAATCGTGGCGCCAGAGGAAAAAAAAAGCCCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTg                                                                                                               <  1:826969/151‑1 (MQ=255)                                          gaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACgcgc                                                                                                         >  1:144046/1‑151 (MQ=255)                                              gttcatcaGGAAAAAACCCACCAGGCCAACCAGCATATACGAGAAAAATCCCCCGCCCGGACCTGCCTGAAAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                           <  2:951759/143‑1 (MQ=255)                                              gttcataggaaaGAAAACCATCAGGCAAACAAAAATAGCCGAGAGCAATCCCCCGCCGGCCCCTCCCTGAAAAATCGTGGCGCCAGAGACAAAAAAAAGCCCTGCCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                                     <  2:514941/141‑1 (MQ=255)                                                 catcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAt                                                                                                  >  1:668515/1‑151 (MQ=255)                                                  atcaGGAAGAAACCCACCGGGCAAAACAACATAAAGGAGACAAACCCCCCGCCCGGCCCTCCGTGAAAAACCGTGGCGCAAGAGAAAAAAAAAAGCCCTTTCCG‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTaa                                                                                                   <  2:40861/149‑1 (MQ=255)                                                     agaaaGAAAACCACCAGCCCAATCAGAATATACGAGAGCATTCCCCCGCCCGGCCCTGCCTGAAAAATCGTTGCGCAAGAGGCAAAAAAAAGCCCTGTACC‑TATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                                                                 <  2:391990/145‑1 (MQ=255)                                                        aaataaccctacggaccAATCAAAATATACAAGACAAATGCCCCACCCGGCCCTGCCTAAAAATTTGTTGCGCAAGATGCAA‑AAAAAGGCCCGTTCCGGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGc                                                                                           <  2:284285/136‑1 (MQ=255)                                                               ccATCGGGCAAAACACCAAATAGGAGAGCAACGCCCCCCCCGGCCCTGCTGGAGAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGcc                                                                                    <  2:935652/151‑1 (MQ=255)                                                                          accAGAATATCCGAGACAAACGCCCCCCCCGGCCCGGCATAAAAAATCGTTGCCCCAGAGACAAAAAAAAGCCCTGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                          <  2:835239/148‑1 (MQ=255)                                                                          aTCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtg                                                                         >  2:558990/1‑151 (MQ=255)                                                                          aTCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtg                                                                         >  2:221485/1‑151 (MQ=255)                                                                                 tataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGCGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                                  >  2:1002868/1‑151 (MQ=255)                                                                                    acaagacCAATCCCCCCCCCGGCCCTGCCTGAAAAACCGTGGCGCCAGGGCAAAAAAAAAGCCCTGTCCG‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                                  <  2:871017/145‑1 (MQ=255)                                                                                    aCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                                               >  2:340201/1‑151 (MQ=255)                                                                                         aactgccccccGCCGGGCCCAGCCGGAGAAGGGGTGGCGCCGGAGCAAAAAAAAAGACCTGTCCT‑ATTGGACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                                               <  2:904605/139‑1 (MQ=255)                                                                                            aacgCCCCGCCCGGCCCTGCCGGAGAATTCGTTGCGCAAGGGCCAAAAAAAAGCCCTGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                        <  1:528128/147‑1 (MQ=255)                                                                                            aaccccccGCCCGGCCCTGCCGGAAAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                        <  1:454186/146‑1 (MQ=255)                                                                                            aaTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                        >  1:40769/1‑150 (MQ=255)                                                                                            aaTCCCCCGCCCGGACCTGCTTGAAAAATCGTGGCGCCAGAGGAAAAAAAAAGCCCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGa                                                         <  1:626982/149‑1 (MQ=255)                                                                                                ccccccccGGCCCTGCCAGAAAAATCGTTGCGCCAGAGACAAAAAAAAGCCCTGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCa                                                    <  1:180473/150‑1 (MQ=255)                                                                                                ccccGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                                                   >  1:772918/1‑151 (MQ=255)                                                                                                  ccgccggAACCTGCATGAAAAATCGTGGCGCAAGAGGAAAAAAAAGGCCTTGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCa                                                    <  2:806273/148‑1 (MQ=255)                                                                                                        ggACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAAt                                           >  1:369532/1‑151 (MQ=255)                                                                                                         gACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                          >  1:843955/1‑151 (MQ=255)                                                                                                         gACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                          >  2:145991/1‑151 (MQ=255)                                                                                                           cctgcctgAAAAACCTTTGCGCCAGAGAAAAAAAAAACACCTGTACC‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                          <  1:86872/149‑1 (MQ=255)                                                                                                           ccgccctaaaaaaTCTTGGCGCCAGAGCAAAAAAAAAGCCCTGCACC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATg                                          <  2:9351/139‑1 (MQ=255)                                                                                                              gcctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc                                     >  1:882286/1‑151 (MQ=255)                                                                                                               cctgAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc                                     <  1:643909/150‑1 (MQ=255)                                                                                                                   aaaaaTCGTTGCGCAAGAGGAAAAAAAAAGCCCGTTACG‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                                    <  2:993172/145‑1 (MQ=255)                                                                                                                    gaAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                               >  1:663262/1‑151 (MQ=255)                                                                                                                     aactcgtcccgcgAGAGAAAAAAAAAAGCCCTGTCCG‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                              <  1:870262/142‑1 (MQ=255)                                                                                                                     aaacggTTGCGCCAGAGCAAAAAAAAAGCCCGGTACG‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                               <  2:959810/145‑1 (MQ=255)                                                                                                                     aaaTCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                               >  1:548667/1‑150 (MQ=255)                                                                                                                     aaaTCGTGGCCCCAGAGGCAAAAAAAAGCCCTGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                               <  2:946506/150‑1 (MQ=255)                                                                                                                      aaTCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGAGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                              >  2:676818/1‑150 (MQ=255)                                                                                                                          tgtGCGCGAGGGCCAAAAAAATGACCTGCACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaa                             <  2:194536/145‑1 (MQ=255)                                                                                                                              cgcAAGAGGAAAAAAAAAGCCCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaaa                          <  1:245874/146‑1 (MQ=255)                                                                                                                                    agcaaaaaaaaaCCCCTGTCCA‑GAGGAACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGa                <  2:1034602/143‑1 (MQ=255)                                                                                                                                        aaaaaaaaGCCCTCTCCG‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCTCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTa           <  2:416810/148‑1 (MQ=255)                                                                                                                                        aaaaaaaaGACCTGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAtt              <  1:424262/145‑1 (MQ=255)                                                                                                                                            aaaaGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGc         >  1:926419/1‑149 (MQ=255)                                                                                                                                              aaGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAc     >  2:39469/1‑151 (MQ=255)                                                                                                                                                 aCCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAt  >  2:105214/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                        AAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAT  >  NC_000913/2248287‑2248581

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