Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,421,374 C→T S17L (TCA→TTA)  yhdY → putative ABC transporter membrane subunit YhdY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,421,3740CT94.3% 104.8 / ‑0.9 37S17L (TCA→TTA) yhdYputative ABC transporter membrane subunit YhdY
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (12/21);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.60e-01

TCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGC  >  NC_000913/3421133‑3421619
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
aatgacctctatgagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCcgcgc                                                                                                                                                                                                                                              <  1:253065‑M2/17‑1 (MQ=255)
                 ggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGaa                                                                                                                                                                                                                             <  1:133512‑M2/34‑1 (MQ=255)
                   tactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGaaa                                                                                                                                                                                                                            >  1:9310‑M2/216‑250 (MQ=255)
                       ggctgaagtagaacagcccctgttggacattgtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacgggatgaactaagtcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgtcttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAAt                                                                                                                                                                                                                        <  2:365501‑M2/39‑1 (MQ=255)
                       ggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATc                                                                                                                                                                                                                       <  1:15369‑M2/40‑1 (MQ=255)
                        gctgaagtagaacagcccctgttggacatggtgaggcaatacacccggggtagccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgatgaaaaaggcgctgctcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCt                                                                                                                                                                                                                      <  1:45044‑M2/41‑1 (MQ=255)
                           gaagtagaacagcccctggtggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGGGCGAAAAAATCTGtt                                                                                                                                                                                                                   >  2:383157‑M2/208‑251 (MQ=255)
                              gtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTtctc                                                                                                                                                                                                                >  1:387847‑M2/205‑251 (MQ=255)
                                          ctgttgtacagggatatgcaatacacccgtggttaccagtcccgtgctgcgctgatgatgggcatcatccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgtcttttttataggggtcacacaggaCATACCTTAAGCCGCGCGATGGTGTGGGTGCCAAAAAATCTGTTCTCCAGCTTGAGc                                                                                                                                                                                                      <  2:252776‑M2/57‑1 (MQ=255)
                                                                      gtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATgg                                                                                                                                                                        <  2:387847‑M2/87‑1 (MQ=255)
                                                                        ggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGt                                                                                                                                                                       <  1:252432‑M2/88‑1 (MQ=255)
                                                                           aaccagacccgtgctgcgcttatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggtgaagcgccttttttataggtgtctcaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGAt                                                                                                                                                                   <  1:189348‑M2/92‑1 (MQ=255)
                                                                                gacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                                                                                                                                              >  1:55381‑M2/155‑251 (MQ=255)
                                                                                                  gatgggcatcaaccgtggcacgctgcgtaaaaaattgaaaaaatacggcatgcaccaattccggctagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggcCATAACTTAAGCCGAGCGATGGTGAGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGc                                                                                                                                            >  2:170946‑M2/137‑251 (MQ=255)
                                                                                                        catcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcatggtagctaaatgcttgattaaaaaggcgctactcggcattgggtagcgcctgttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACt                                                                                                                                      <  2:257177‑M2/121‑1 (MQ=255)
                                                                                                                gtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCAttt                                                                                                                              >  1:217632‑M2/123‑251 (MQ=255)
                                                                                                                           cataaaaaattgaaaaaaaacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTGCTCCAGCTGGAGCAGTAGCCTGCTGACTAGTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTaa                                                                                                                   <  2:373390‑M2/140‑1 (MQ=255)
                                                                                                                                             atgaATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAAcgcg                                                                                                 >  2:404162/1‑251 (MQ=255)
                                                                                                                                                ggcaagaactaattcatgtaagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCCTAACTTAAGCAGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGAAGGCTACCTGGGTTGGCTCAACGCGTGc                                                                                              >  2:367891‑M2/91‑251 (MQ=255)
                                                                                                                                                          aattcaggttagctaaatgcttgattaaaaaggcgctactcagcatggggaagagccttttttataggtgtcaaaaaggaCATAACTAAAGCCGCGCGATGGGGTGGGTGCGACAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGAATATTGGCTGCATATGGTTAATGCGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAGCTGGGTTGGCTAAACGCGTGCCGACTGcac                                                                                     >  1:200484‑M2/81‑250 (MQ=255)
                                                                                                                                                           attcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGcaca                                                                                    >  1:155472‑M2/80‑250 (MQ=255)
                                                                                                                                                                   tagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCgg                                                                            <  1:22347‑M2/179‑1 (MQ=255)
                                                                                                                                                                    cgcgATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTg                                                                          <  1:404162/251‑1 (MQ=255)
                                                                                                                                                                        aaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTg                                                                      >  1:465789‑M2/67‑251 (MQ=255)
                                                                                                                                                                                gattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTc                                                               <  1:326454‑M2/192‑1 (MQ=255)
                                                                                                                                                                                    aaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCAt                                                          >  1:79385‑M2/55‑251 (MQ=255)
                                                                                                                                                                                            gcgactcggcatggggaagcgccttttttatagggggcacaaaggaCATACCTGAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGAAATAGCCTGCTGATTATTGGCTGCATATGGTTGATGTGGTAACTGATCCCCCCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCCCAAAAGCCGGTGCCTGTTGGGTCTGCATCCtcgaac                                                  <  2:430749‑M2/205‑1 (MQ=255)
                                                                                                                                                                                               actcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCacgaacga                                                <  2:311618‑M2/207‑1 (MQ=255)
                                                                                                                                                                                               actcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCacgaacga                                                <  1:338783‑M2/207‑1 (MQ=255)
                                                                                                                                                                                                  cggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTg                                            <  2:377917‑M2/211‑1 (MQ=255)
                                                                                                                                                                                                  cggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCGCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTg                                            <  2:41908‑M2/211‑1 (MQ=255)
                                                                                                                                                                                                      atggggaggcgccttttttataggtgtctcanaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACTCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCa                                        <  2:236677‑M2/215‑1 (MQ=255)
                                                                                                                                                                                                      atggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTc                                         >  1:206373‑M2/37‑250 (MQ=255)
                                                                                                                                                                                                              gcgctgtttttataggtggcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCGGACTATTGGCTGCAGATGGTTGATGTGGTAACTGATCCCACCGTTGCTGAACTGGGCTTTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgt                                <  1:215177‑M2/223‑1 (MQ=255)
                                                                                                                                                                                                                                         aCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCt     <  2:5595‑M2/250‑1 (MQ=255)
                                                                                                                                                                                                                                            tAACTTAAGCCGCGCGTTGGGGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGCACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCCACGCTGGc  >  2:445301/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaaCTTAAGCCGCGCGATGGTGTGGGTGCGANAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGc  <  2:206373/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
TCTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGC  >  NC_000913/3421133‑3421619

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