Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,934 T→C 6.3% Q6R (CAA→CGA)  ybiA ← DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913832,9340TC 6.3% 215.0 / 3.9 63Q6R (CAA→CGA) ybiADUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (2/2);  total (31/32)
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.40e-01

TTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTCTCATG  >  NC_000913/832787‑833072
                                                                                                                                                   |                                                                                                                                          
ttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCttt                                                                                                                                            <  1:779356/148‑1 (MQ=255)
ttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGGGTTGGATCCttt                                                                                                                                            <  2:236270/148‑1 (MQ=255)
         tGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGaa                                                                                                                                   <  2:986611/148‑1 (MQ=255)
                    ttCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTaaa                                                                                                                        <  1:566507/148‑1 (MQ=255)
                         aGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCa                                                                                                                   >  4:706899/1‑148 (MQ=255)
                                        tCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTa                                                                                                    <  3:706899/148‑1 (MQ=255)
                                             tCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCg                                                                                               >  4:280861/1‑148 (MQ=255)
                                              cAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGa                                                                                              <  3:280861/148‑1 (MQ=255)
                                                   ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGc                                                                                                >  2:244391/1‑141 (MQ=255)
                                                   ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGc                                                                                                <  1:244391/141‑1 (MQ=255)
                                                   ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGc                                                                                         >  3:585236/1‑148 (MQ=255)
                                                       tGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTg                                                                                                 >  1:984465/1‑136 (MQ=255)
                                                       tGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTg                                                                                                 <  2:984465/136‑1 (MQ=255)
                                                             aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                                                                                                   <  4:50442/128‑1 (MQ=255)
                                                             aGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                                                                                                   >  3:50442/1‑128 (MQ=255)
                                                               gcggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGt                                                                                                             >  2:627445/1‑116 (MQ=255)
                                                               gcggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGt                                                                                                             <  1:627445/116‑1 (MQ=255)
                                                                       gTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAAcc                                                                     <  4:585236/148‑1 (MQ=255)
                                                                           gAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGc                                                                                                                    <  2:145858/97‑1 (MQ=255)
                                                                           gAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGc                                                                                                                    >  1:145858/1‑97 (MQ=255)
                                                                             aaaaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATACTTTGTGCTCGAACgg                                                                                                                                <  2:210151/83‑1 (MQ=255)
                                                                             aaaaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATACTTTGTGCTCGAACgg                                                                                                                                >  1:210151/1‑83 (MQ=255)
                                                                                  cACCGTAGTCATCGCTGGTGCTGTAGAACTTTATGATGCTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAActa                                                                    <  4:377582/138‑3 (MQ=255)
                                                                                  cACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCa                                                                    >  3:377582/1‑138 (MQ=255)
                                                                                                  ggTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCg                                                                                                                     <  3:121441/73‑1 (MQ=255)
                                                                                                  ggTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCg                                                                                                                     >  4:121441/1‑73 (MQ=255)
                                                                                                          aGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTcc                                                                                    >  3:145514/1‑98 (MQ=255)
                                                                                                          aGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTcc                                                                                    <  4:145514/98‑1 (MQ=255)
                                                                                                          aGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGaaaa                                                                        <  1:953067/110‑1 (MQ=255)
                                                                                                          aGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGaaaa                                                                        >  2:953067/1‑110 (MQ=255)
                                                                                                            aaaTTTATGCTGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCt                                                                                                                                       >  2:559317/1‑45 (MQ=255)
                                                                                                            aaaTTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCt                                                                                                                                       <  1:559317/45‑1 (MQ=255)
                                                                                                              aTTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGc                                >  1:868851/1‑146 (MQ=255)
                                                                                                              aTTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTACTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGc                                <  2:868851/146‑1 (MQ=255)
                                                                                                                ttATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGtt                            >  2:472120/1‑148 (MQ=255)
                                                                                                                    gatgGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTa                                                                 <  2:706142/107‑1 (MQ=255)
                                                                                                                    gatgGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTa                                                                 >  1:706142/1‑107 (MQ=255)
                                                                                                                           tCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTgg                                 <  3:406495/132‑1 (MQ=255)
                                                                                                                           tCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTgg                                 >  4:406495/1‑132 (MQ=255)
                                                                                                                            cTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                                                                                                   <  3:898123/65‑1 (MQ=255)
                                                                                                                            cTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                                                                                                   >  4:898123/1‑65 (MQ=255)
                                                                                                                             ttGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtat                         <  2:781044/138‑1 (MQ=255)
                                                                                                                             ttGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAAAAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtat                         >  1:781044/1‑138 (MQ=255)
                                                                                                                              tGCATGACGTGTTGGATCCTTCGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                           <  3:477243/119‑1 (MQ=255)
                                                                                                                              tGCATGACGTGTTGGATCCTTCGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                           <  2:611486/119‑1 (MQ=255)
                                                                                                                              tGCATGACGTGTTGGATCCTTCGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                           >  1:611486/1‑119 (MQ=255)
                                                                                                                              tGCATGACGTGTTGGATCCTTCGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                           >  4:477243/1‑119 (MQ=255)
                                                                                                                                      gtgtTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTtct      >  1:878168/1‑148 (MQ=255)
                                                                                                                                      gtgtTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTtct      <  2:878168/148‑1 (MQ=255)
                                                                                                                                       tgtTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTtc       >  1:297320/1‑146 (MQ=255)
                                                                                                                                       tgtTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTtc       <  2:297320/146‑1 (MQ=255)
                                                                                                                                        gtTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                           >  1:728911/1‑109 (MQ=255)
                                                                                                                                        gtTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGcc                                           <  2:728911/109‑1 (MQ=255)
                                                                                                                                        gtTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTgg                                 >  3:1021999/1‑119 (MQ=255)
                                                                                                                                        gtTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTgg                                 <  4:1021999/119‑1 (MQ=255)
                                                                                                                                          tGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                                                                                                   <  1:1101279/51‑1 (MQ=255)
                                                                                                                                          tGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                                                                                                   >  2:1101279/1‑51 (MQ=255)
                                                                                                                                          tGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCca                                          <  4:1107905/108‑1 (MQ=255)
                                                                                                                                          tGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCca                                          >  3:1107905/1‑108 (MQ=255)
                                                                                                                                            gATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtat                         >  3:400111/1‑123 (MQ=255)
                                                                                                                                            gATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtat                         <  4:400111/123‑1 (MQ=255)
                                                                                                                                               ccTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTc                                                                                     <  4:1081666/60‑1 (MQ=255)
                                                                                                                                               ccTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTc                                                                                     >  3:1081666/1‑60 (MQ=255)
                                                                                                                                               ccTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTg                                  <  2:1058099/111‑1 (MQ=255)
                                                                                                                                               ccTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTg                                  >  1:1058099/1‑111 (MQ=255)
                                                                                                                                                 tttGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGtt                                                                                                     >  3:854969/1‑42 (MQ=255)
                                                                                                                                                 tttGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGtt                                                                                                     <  4:854969/42‑1 (MQ=255)
                                                                                                                                                 tttGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTCTCa    <  4:200685/139‑1 (MQ=255)
                                                                                                                                                 tttGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTCTCa    >  3:200685/1‑139 (MQ=255)
                                                                                                                                                 tttCTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAATAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTCTCATg  <  2:966595/141‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
TTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTCTCATG  >  NC_000913/832787‑833072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: