Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 216,568 C→T A146T (GCC→ACC)  yaeF ← putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913216,5680CT97.6% 285.5 / ‑2.0 84A146T (GCC→ACC) yaeFputative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (49/33);  total (50/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

ACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAAC  >  NC_000913/216475‑216664
                                                                                             |                                                                                                
aCGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAcgc                                                                                            >  2:1251807/1‑100 (MQ=255)
 cGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAAcgcg                                                                                           >  1:213451/1‑100 (MQ=255)
  ggATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCgg                                                                                          >  2:480500/1‑100 (MQ=255)
   gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                         >  1:665829/1‑100 (MQ=255)
   gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                         <  1:1320555/100‑1 (MQ=255)
    aTTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTg                                                                                        <  1:1038290/100‑1 (MQ=255)
    aTTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTg                                                                                        <  1:841946/100‑1 (MQ=255)
       cAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGata                                                                                     <  1:923455/100‑1 (MQ=255)
        aGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGatat                                                                                    >  2:1369303/1‑100 (MQ=255)
        aGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGatat                                                                                    >  2:241076/1‑100 (MQ=255)
        aGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGatat                                                                                    >  2:874086/1‑100 (MQ=255)
            aGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgt                                                                                >  1:501379/1‑100 (MQ=255)
            aGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgt                                                                                >  1:237757/1‑100 (MQ=255)
             gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                               <  1:1139325/100‑1 (MQ=255)
                caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCt                                                                            >  2:576786/1‑100 (MQ=255)
                caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTGTATTGGTAAACGCGGTGATATCTGTGGCt                                                                            >  2:281176/1‑100 (MQ=255)
                   cTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgt                                                                         <  2:1555662/100‑1 (MQ=255)
                     ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                       >  1:154871/1‑100 (MQ=255)
                     ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                       >  2:789114/1‑100 (MQ=255)
                     ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                       >  1:1394739/1‑100 (MQ=255)
                      gCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGc                                                                      >  2:544314/1‑100 (MQ=255)
                          gTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgggg                                                                  <  2:709510/100‑1 (MQ=255)
                            cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                                <  2:824181/100‑1 (MQ=255)
                            cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                                <  1:997108/100‑1 (MQ=255)
                            cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                                <  2:501378/100‑1 (MQ=255)
                            cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                                >  2:1187329/1‑100 (MQ=255)
                            cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                                >  2:64617/1‑100 (MQ=255)
                             aCCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTaa                                                               >  2:1136601/1‑100 (MQ=255)
                               cATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGa                                                             >  2:215615/1‑100 (MQ=255)
                                 tAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATc                                                           >  2:475560/1‑100 (MQ=255)
                                  aaaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATcc                                                          <  1:1446925/100‑1 (MQ=255)
                                   aacggAAGAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                         >  2:584792/4‑100 (MQ=255)
                                   aaGGGAATAAATTCGATAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                         >  2:496794/1‑100 (MQ=255)
                                   aaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                         >  1:1057796/1‑100 (MQ=255)
                                   aaGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                         >  2:125612/1‑100 (MQ=255)
                                    aGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGAGCCgg                                                        <  1:1369310/100‑1 (MQ=255)
                                     tggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCggg                                                       <  1:877786/99‑1 (MQ=255)
                                     gggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCggg                                                       <  1:634853/100‑1 (MQ=255)
                                      ggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGAGTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGa                                                      >  2:877781/1‑100 (MQ=255)
                                       gAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCg                                                         >  2:132121/1‑96 (MQ=255)
                                         aTAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTc                                                   >  1:1168633/1‑100 (MQ=255)
                                          tAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTc                                                   >  2:1388229/1‑99 (MQ=255)
                                          tAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCg                                                  <  1:1526691/100‑1 (MQ=255)
                                           tctttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGt                                                 <  2:705649/97‑1 (MQ=255)
                                           aaaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGt                                                 <  2:59807/100‑1 (MQ=255)
                                           aaaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGt                                                 >  1:118753/1‑100 (MQ=255)
                                            aaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTa                                                <  1:1214991/100‑1 (MQ=255)
                                             aTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTa                                                >  2:1033452/1‑99 (MQ=255)
                                              ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTaga                                              >  1:705651/1‑98 (MQ=255)
                                              ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                              <  1:1195558/100‑1 (MQ=255)
                                                cGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGAc                                            <  1:1334534/100‑1 (MQ=255)
                                                 gACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACg                                           <  2:806221/100‑1 (MQ=255)
                                                   cAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTCAGACGaa                                         >  1:650781/1‑100 (MQ=255)
                                                     aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaaa                                       <  1:549541/100‑1 (MQ=255)
                                                     aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaaa                                       >  1:937705/1‑100 (MQ=255)
                                                      tGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAg                                      >  1:542869/1‑100 (MQ=255)
                                                        ccGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCt                                    >  2:1518492/1‑100 (MQ=255)
                                                          gcgATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTa                                  >  2:983695/1‑100 (MQ=255)
                                                            gATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATc                                >  1:1461644/1‑100 (MQ=255)
                                                              tAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCAc                              >  2:984980/1‑100 (MQ=255)
                                                                   ataACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGc                         <  1:64617/100‑1 (MQ=255)
                                                                    taACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCt                        <  2:570368/100‑1 (MQ=255)
                                                                     aaCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCtt                       <  2:937700/100‑1 (MQ=255)
                                                                     aaCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCta                       >  1:570369/1‑99 (MQ=255)
                                                                       ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa                     >  1:1111894/1‑100 (MQ=255)
                                                                       ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa                     >  1:1228425/1‑100 (MQ=255)
                                                                       ccGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCAGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCa                     >  1:632071/1‑100 (MQ=255)
                                                                        cGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCAt                    <  1:873191/100‑1 (MQ=255)
                                                                              cTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttt              >  1:1515534/1‑100 (MQ=255)
                                                                               tttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttt             >  1:355703/1‑100 (MQ=255)
                                                                                ttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttttt            >  1:1579272/1‑100 (MQ=255)
                                                                                ttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttttt            <  1:894610/100‑1 (MQ=255)
                                                                                 tGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttttt           <  2:1536025/100‑1 (MQ=255)
                                                                                  gATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCGTATCACTATGCTTCATCGCTTTTTTa          <  2:1388019/100‑1 (MQ=255)
                                                                                   aTTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAGGCTTATCACTATGCTTCATCGCTTTTTTaa         >  2:1575217/1‑100 (MQ=255)
                                                                                   aTTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTaa         >  1:544853/1‑100 (MQ=255)
                                                                                   aTTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTaa         >  1:223517/1‑100 (MQ=255)
                                                                                     tttATTGGTATACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAATAGTTTATCACTATGCTTCATCGCTTTTTTAAgg       <  2:1038791/100‑1 (MQ=255)
                                                                                     tttATTGGTAAACACGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAg        <  2:1326453/99‑1 (MQ=255)
                                                                                      ttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAggg      <  1:997516/100‑1 (MQ=255)
                                                                                      ttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAggg      <  1:137366/100‑1 (MQ=255)
                                                                                         ttGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGaaa   >  1:507714/1‑100 (MQ=255)
                                                                                          tGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAAc  >  1:1039756/1‑100 (MQ=255)
                                                                                          tGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAAc  <  1:792664/100‑1 (MQ=255)
                                                                                             |                                                                                                
ACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAGGGAAAC  >  NC_000913/216475‑216664

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: