Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,748	0	C	A	41.4%	80.5 / 30.7	58	E55* (GAG→TAG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (30/4);  new base A (0/24);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.60e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGG  >  NC_000913/3177601‑3177893                                                                                                                                                     |                                                                                                                                                  tttCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCaa                                                                                                                                                 >  2:882620/1‑150 (MQ=255)  ttCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCaa                                                                                                                                                 <  1:436955/149‑1 (MQ=255)     aCCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt                                                                                                                                            >  2:327022/1‑151 (MQ=255)            tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                     >  2:995628/1‑151 (MQ=255)                  atcCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCAcc                                                                                                                                >  1:337624/1‑150 (MQ=255)                      cGGGAACCATATCCCGTAGCCAAGGGGGTTCGCTGGGTTCTTACGCCGCCACCCGAATCTGCTCAATCGGCACAAATCAAACACGCGCTTGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCCCCCTCAAAAATTTCCCGCCCCCCCTc                                                                                                                              >  1:611876/1‑148 (MQ=255)                         cAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACt                                                                                                                        >  1:916967/1‑151 (MQ=255)                                    aGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt                                                                                                             >  2:686243/1‑151 (MQ=255)                                     gTAGCCAGGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGtt                                                                                                            <  2:410614/151‑1 (MQ=255)                                        gCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGa                                                                                                          >  1:522333/1‑150 (MQ=255)                                         ccAGGGGGTTGCGCGGGTGCCGTCCGCGGCAACCCGACTCTGCTCAACCAGCACACCTACATCACCCACTGGGTCAAAGGGT‑AGCAAGCCTGCGGCGTGACCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGtt                                                                                                            <  1:282717/147‑1 (MQ=255)                                                 tttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCACCACGCACGGGGTAAAAGGGT‑AGCAAGCCTGGGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGAt                                                                                                 <  2:754207/150‑1 (MQ=255)                                                  ttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                >  2:95953/1‑150 (MQ=255)                                                  tccGCTGGTGTCAGACGCGGAAATCCGAATCTGCTCAATCAGCACAACTTCACCCCCCACTGGGTCAAGGGGT‑AGCAAGACTGCGGCTTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                <  1:149946/148‑1 (MQ=255)                                                   tCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                               <  1:464520/150‑1 (MQ=255)                                                       tGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa                                                                                           >  2:247130/1‑150 (MQ=255)                                                        ggTGTCGTACGCGGAAATCCGAATCTGCTCAACCAGCACAACTTCATCACGCACGGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                                                                                          <  1:874550/150‑1 (MQ=255)                                                              gtccgcgGCAATCCAAATCTCCACAATCAGCACAACTTCATCCCCCACTGGTTCAAGGGGT‑AGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa                                                                                   <  2:860756/148‑1 (MQ=255)                                                                 cgcgGCAATCCGAATCTGCTCAATCAGCACAACTGCACCACGCACTGGGTCAAAGGGT‑AGCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat                                                                                     <  1:304273/146‑1 (MQ=255)                                                                 cgcgGCAATCCGAATCTGCTCAAACAGCACAACTTCACCCCGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAg                                                                                 <  1:246128/150‑1 (MQ=255)                                                                      cAACCCGACCCTGCTCAATCACCACAACTTCATCCCGCA‑TTGGTTAAAGGTTCCAAAAGTCGGCGGCGTGCCCCCGGCAAAAAATCACCCGCCGCACCTCATGACGCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa                                                                                   <  1:572392/143‑1 (MQ=255)                                                                        aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg                                                                         >  2:665338/1‑151 (MQ=255)                                                                            gAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaa                                                                      >  1:615137/1‑150 (MQ=255)                                                                            gAATCTGCTCAACCAGCACACCTCCACCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg                                                                         <  1:338154/147‑1 (MQ=255)                                                                             aatttccTCAATCAGCACAACTTCATCACCCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaa                                                                    <  2:881828/145‑1 (MQ=255)                                                                                 tGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGc                                                                >  2:174954/1‑151 (MQ=255)                                                                                 tGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAAGATGACGAAATCTATAAAGATCTAATGAAAAAAAg                                                                 >  1:333317/1‑150 (MQ=255)                                                                                  gCTCAATCAGCACAACTTCACCACGCACTGGGTAAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAg                                                                 <  2:877062/149‑1 (MQ=255)                                                                                   cccAATCAGCACAACTTCACCCCGCACTGGGTAAAGGGGT‑AGCAAGACCGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcg                                                              <  2:610139/149‑1 (MQ=255)                                                                                   cccAACCAGCACATCTTCATCCCCCACTGGGTAAAAGGGT‑AGCAAGACTGCGGCGTGCCCCGGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaaa                                                                  <  1:644921/145‑1 (MQ=255)                                                                                   cTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcg                                                              >  1:306339/1‑151 (MQ=255)                                                                                    tCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgc                                                             >  1:222432/1‑151 (MQ=255)                                                                                        tCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                                          >  2:764941/1‑150 (MQ=255)                                                                                         cAGCACAACTTCATCCCGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                                          <  2:520835/149‑1 (MQ=255)                                                                                          aGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaa                                                        >  2:804468/1‑150 (MQ=255)                                                                                            cacaACTTCACCCCGCACGGGGTCAAAGGGT‑AGAAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaaa                                                       <  1:217564/149‑1 (MQ=255)                                                                                             acaACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtg                                                    >  1:554654/1‑151 (MQ=255)                                                                                              caACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                   >  1:626084/1‑151 (MQ=255)                                                                                                aCTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                 >  1:108012/1‑151 (MQ=255)                                                                                                      tcaCGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTg                                           >  2:259346/1‑151 (MQ=255)                                                                                                            aCTGGGTCAAAGGGT‑AGCAAGACTGCGGCTTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataa                                     <  2:782170/151‑1 (MQ=255)                                                                                                            aCTGGGTAAAGGGGT‑AGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataa                                     <  1:798552/151‑1 (MQ=255)                                                                                                               gggTCAAAGGGT‑AGCAAGACTGCGGGGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataa                                     <  2:380322/148‑1 (MQ=255)                                                                                                               gggTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAttt                                  >  1:298292/1‑151 (MQ=255)                                                                                                                    aaaGGGT‑AGCAAGCCTGCGGCGTGACCGCGCTCAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACa                             <  2:207438/151‑1 (MQ=255)                                                                                                                     aaGGGT‑AGAAGGCCTGCGGCGTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACAt                            <  2:559259/151‑1 (MQ=255)                                                                                                                       gggT‑AGCAAGATTGCGGAGTGACCGCATTAAAAAATTCACCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAACAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCg                          <  1:189379/151‑1 (MQ=255)                                                                                                                        ggT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGt                         >  2:131733/1‑151 (MQ=255)                                                                                                                          t‑AGAAAGACGGCGGCGGGCCCCCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtt                        <  1:851186/150‑1 (MQ=255)                                                                                                                             gCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGttttt                     >  1:971262/1‑151 (MQ=255)                                                                                                                               aaGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTATACATCGTTTTTGc                   >  2:874378/1‑151 (MQ=255)                                                                                                                               aaGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGttttt                     <  1:887210/149‑1 (MQ=255)                                                                                                                                aggcTGCGGCGTAACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGcc                  <  2:638133/148‑1 (MQ=255)                                                                                                                                    tGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAAt              >  1:372085/1‑151 (MQ=255)                                                                                                                                    tGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAAt              >  2:736385/1‑151 (MQ=255)                                                                                                                                     gcggcgTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAAt              >  2:381475/1‑150 (MQ=255)                                                                                                                                       ggcgTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTa           <  2:234303/151‑1 (MQ=255)                                                                                                                                                gcgcTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAgg  >  2:152261/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGT‑AGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGG  >  NC_000913/3177601‑3177893

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