Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,740,735 A→G L139P (CTG→CCG)  aroF ← 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,740,7350AG100.0% 196.8 / NA 56L139P (CTG→CCG) aroF3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tyrosine‑repressible
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/33);  total (23/33)

TTAAAACCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCAG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCG  >  NC_000913/2740591‑2740874
                                                                                                                                                |                                                                                                                                              
ttAAAACCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑cccc                                                                                                                                        >  1:67085/1‑150 (MQ=255)
ttAAAACCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑ccc                                                                                                                                         <  2:156615/150‑1 (MQ=255)
  aaaaCCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCata                                                                                                                                      <  2:416992/151‑1 (MQ=255)
   aaaCCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCatat                                                                                                                                     <  2:103745/151‑1 (MQ=255)
     accaaccGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTc                                                                                                                                   <  1:180993/151‑1 (MQ=255)
       caaccGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCAc                                                                                                                                 <  2:311745/151‑1 (MQ=255)
         accGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑TCCCTAttttcc                                                                                                                                >  2:486785/1‑146 (MQ=255)
          ccGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAg                                                                                                                              <  2:347226/151‑1 (MQ=255)
           cGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTTACCaa                                                                                                                              >  2:619543/1‑149 (MQ=255)
            ggCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCt                                                                                                                            <  2:729887/151‑1 (MQ=255)
              cATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑TCCCCATTTTCCAGc                                                                                                                             >  2:658376/1‑148 (MQ=255)
              cATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGGATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCa                                                                                                                          <  1:393427/151‑1 (MQ=255)
                          ggAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCAcc                                                                                                                                >  2:171689/1‑133 (MQ=255)
                          ggAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCa                                                                                                                               <  1:171656/134‑1 (MQ=255)
                           gAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTAc                                                                                                             <  1:631162/151‑1 (MQ=255)
                              gCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCTATTTACCAACCCAAGCAGCAAttttagc                                                                                                           >  2:81875/1‑146 (MQ=255)
                               ccATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTTACCAGCTCCAGCAGCAAattttg                                                                                                            >  1:302925/1‑145 (MQ=255)
                               ccATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTAcgcgc                                                                                                         <  2:701934/151‑1 (MQ=255)
                                cATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTAcgcgcg                                                                                                        <  2:313666/151‑1 (MQ=255)
                                  tttCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCTATTTACCAGCTCAAGCAGCAATTTAcgcgc                                                                                                         >  2:38855/1‑148 (MQ=255)
                                      aCGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGc                                                                                                  <  1:299155/151‑1 (MQ=255)
                                      aCGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑ACCACATTTTCCCGCTCAAACCGCAATTTTCACGCGGTCt                                                                                                    >  1:57594/1‑149 (MQ=255)
                                       cGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑TCCATTTTCACCAGCTCCAACAACAAATTACGCGCGctct                                                                                                    >  1:759708/1‑148 (MQ=255)
                                       cGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCa                                                                                                 <  2:7515/151‑1 (MQ=255)
                                                tGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑TACCCCAATTAACC‑GCCCCAACAGCAAATTACGAGCGagcgagcggccc                                                                                            >  1:438666/1‑136 (MQ=255)
                                                 gCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCtt                                                                                        <  1:276122/150‑1 (MQ=255)
                                                  cGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCCATTCACCAACTCCAACAGCAAATTATGCGGGCGCTGCCGCCCgg                                                                                           >  1:231918/1‑146 (MQ=255)
                                                       ccGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACAt                                                                                  <  2:288996/150‑1 (MQ=255)
                                                        cGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCTATTTACCCGCCCAAACAGCAATTTACGCGCGGTCTGCAGCCCCGCctct                                                                                      >  1:170575/1‑145 (MQ=255)
                                                        cGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCTATTCACCACCTCAAACCGCAATTTATCCGCGGTCTGCAGCCCCGCctcttcca                                                                                  >  2:381770/1‑145 (MQ=255)
                                                        cGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACa                                                                                   <  1:381722/148‑1 (MQ=255)
                                                                gAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGcc                                                                        <  1:602034/151‑1 (MQ=255)
                                                                    aCCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCAATTTACCCGCCCCAGCAGCAAATTATGCGCGGTCTGCAGCCCGGCCTCTACATCCAAAAAGccaccc                                                                    >  2:74364/1‑151 (MQ=255)
                                                                      cAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATTTTCACCAGCTCCAACAGCAATTTACGCCCGCTCTGCCGCCCGGCCTCTTCATCAAAACAGccatcc                                                                    >  1:75179/1‑149 (MQ=255)
                                                                       aaTTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCata                                                                 <  2:333347/151‑1 (MQ=255)
                                                                           gCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATggg                                                             <  1:74333/151‑1 (MQ=255)
                                                                              gACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATgggg                                                            <  2:45804/149‑1 (MQ=255)
                                                                              gACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGAt                                                          <  1:744050/151‑1 (MQ=255)
                                                                                     tAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGGAT‑CCATTATTACCAAGCCCAAGAAG‑AAATTATGCGCGGGCTGGCGCCCGGCTTCTTCATTAAAAAAGCCATCCATATGGGGATC‑Gctaa                                                   >  2:407454/1‑150 (MQ=255)
                                                                                       aaCAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTTCCCAGCTCAAGCAGCAAATTTCGCGCGATCTGCAGGCCGGCCTCTACATCAAAAAAGACATCCCTAATGGGGGC‑Gtttatt                                                 >  1:53858/1‑150 (MQ=255)
                                                                                         cAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑TCCCTATTTACCCACTCAAACAGCAACTTTCGCGCGATCTGCCGCCCGGCCTCTTCATCCAAAGAGACCTCCCTATAGGGGTG‑Gttaattta                                               >  1:583850/1‑148 (MQ=255)
                                                                                           gATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCTATTTACCAGCTCAAACAGCAATTTATGCGCGGTCTGCCGCCCGGCCTCTTCCACAAAAAAGACCTTCATATGGGGGTC‑Gtttat                                                  >  1:500760/1‑145 (MQ=255)
                                                                                            aTCGCCCAGGTATTGCGGGCTATTCGGATCTAACTCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCt                                           <  1:172574/151‑1 (MQ=255)
                                                                                              cGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCttt                                         <  1:441859/151‑1 (MQ=255)
                                                                                                  cAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGC‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGTTCTGCAGCCCGGTTTCTACATCAAAATAGCCATCCATATGGGGATC‑GTTAATTAACCCGTTCCa                                      <  2:498170/150‑1 (MQ=255)
                                                                                                         tGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACa                               <  2:170608/150‑1 (MQ=255)
                                                                                                                 aTTCGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACgg                      <  1:272252/151‑1 (MQ=255)
                                                                                                                    cGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACgggg                    <  2:449140/150‑1 (MQ=255)
                                                                                                                    cGGATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGt                   <  2:17012/151‑1 (MQ=255)
                                                                                                                     ggATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCCATTTACCACCCCCAGCAGCAATTTTCACGCGGCCTTTCGGCCCGCTTTCCCATTAAAAAAAACACCTATATCGGGGTCGGTTAATTATCCTTTTCCAG‑CGACAgggggaacggg                    >  1:749830/1‑142 (MQ=255)
                                                                                                                     ggATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGtt                  <  1:560693/151‑1 (MQ=255)
                                                                                                                     ggATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGtt                  <  1:640288/151‑1 (MQ=255)
                                                                                                                      gATCTAACGCTTCCGTCGCCAGTGGCGG‑T‑CCCCTATTCACCAGCTCAAACACCAATTTTCGCGCGATCTGCAGCCCCGCCTCTACATCCAAAAAAGCAACCCTATTGGGGTC‑GTTTATTAAACCTTTCCCGCCGACAGTGGGACggggggt                 >  1:631538/1‑148 (MQ=255)
                                                                                                                           aaCGCTTCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTc              <  2:523109/149‑1 (MQ=255)
                                                                                                                                 tCCGTCGCCAGTGGCGG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATa        <  1:157857/149‑1 (MQ=255)
                                                                                                                                      cGCCAGTGGCGG‑T‑CCCATATTTACCAGCTCCAGCAGCAAATTTCGCGCGGACTGCCGCCCGGCTTCTACATCAAAAAAACCAACCAAATGGGGGTC‑GTTAATTTAACCTTTTCCGCCGCCAGTGGGACGGGGGGTTTCAAAATAAAcgcg  >  1:523508/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TTAAAACCAACCGGCATGGAAAGCCCGGAGGCCATTTCACGGTGAGTTTGCGATTCCGTTGTACGAGCACCAATTGCTGACCAGCTAAACAGATCGCCCAGGTATTGCGGGCTATTCGGATCTAACGCTTCCGTCGCCAGTGGCAG‑T‑CCCATATTCACCAGCTCAAGCAGCAATTTACGCGCGATCTGCAGCCCGGCTTCTACATCAAAAGAGCCATCCATATGGGGATC‑GTTAATTAACCCTTTCCAGCCGACAGTGGTACGGGGTTTTTCAAAATAGACGCG  >  NC_000913/2740591‑2740874

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