Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,593,917 G→A A648V (GCG→GTG)  tmcA ← elongator methionine tRNA (ac4C34) acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,593,9170GA96.3% 266.6 / ‑0.7 82A648V (GCG→GTG) tmcAelongator methionine tRNA (ac4C34) acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (27/51);  total (30/52)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GACCATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGTA  >  NC_000913/2593820‑2594013
                                                                                                 |                                                                                                
taCCATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCAcc                                                                                                <  2:104859/99‑1 (MQ=255)
  ccATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTg                                                                                              >  1:1803060/1‑100 (MQ=255)
  ccATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTg                                                                                              >  2:454823/1‑100 (MQ=255)
   cATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTgc                                                                                             >  1:1153558/1‑100 (MQ=255)
     tGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTgcgc                                                                                           >  2:1115219/1‑100 (MQ=255)
     tGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTgcgc                                                                                           >  1:121428/1‑100 (MQ=255)
      gcgaTAAAATTGAAGCAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGcgcg                                                                                          <  1:504807/98‑1 (MQ=255)
      gagaTAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGcgcg                                                                                          <  1:504749/100‑1 (MQ=255)
         aTAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGcc                                                                                       <  1:454821/100‑1 (MQ=255)
         aTAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGcc                                                                                       >  1:1547275/1‑100 (MQ=255)
          tAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCct                                                                                      >  1:1875024/1‑100 (MQ=255)
           aaaaTTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCctc                                                                                     <  1:1120972/100‑1 (MQ=255)
           aaaaTTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCctc                                                                                     <  1:427875/100‑1 (MQ=255)
            aaaTTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCctct                                                                                    >  2:1440716/1‑100 (MQ=255)
               ttGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGCTCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTcc                                                                                 <  1:1366316/100‑1 (MQ=255)
               ttGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTcc                                                                                 <  2:561164/100‑1 (MQ=255)
               ttGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTcc                                                                                 >  2:735485/1‑100 (MQ=255)
                tGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTccc                                                                                <  2:501603/100‑1 (MQ=255)
                     gAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGa                                                                           >  2:1018774/1‑100 (MQ=255)
                     gAGTCAGCGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGa                                                                           <  2:467318/100‑1 (MQ=255)
                       gTCAATGAAACAATTGCCATTGCGTTATGCTATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGAc                                                                          <  2:362179/99‑1 (MQ=255)
                        tCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACgg                                                                        <  2:621962/100‑1 (MQ=255)
                        tCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACgg                                                                        <  2:458291/100‑1 (MQ=255)
                         cAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGa                                                                       >  1:494602/1‑100 (MQ=255)
                          aGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGAc                                                                      >  2:1504165/1‑100 (MQ=255)
                           gTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGAcc                                                                     <  2:1933900/100‑1 (MQ=255)
                           gTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGAcc                                                                     <  1:1147768/100‑1 (MQ=255)
                            tGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCa                                                                    <  2:1444283/100‑1 (MQ=255)
                             gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                   >  1:1048453/1‑100 (MQ=255)
                             gAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAg                                                                   >  1:1944738/1‑100 (MQ=255)
                               aataatTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCCCCTGCGCGGCCTCTTCCCGCTGACGGAcc                                                                     >  1:32924/1‑96 (MQ=255)
                               aataatTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTa                                                                 <  1:735484/100‑1 (MQ=255)
                                ataatTGCCATTGCGTTATGCGATCGCGCAGACGCTGAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTaa                                                                >  2:1371773/1‑100 (MQ=255)
                                 taatTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAAc                                                               <  2:1069333/100‑1 (MQ=255)
                                 taatTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAAc                                                               <  1:1432674/100‑1 (MQ=255)
                                    ttGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATc                                                            <  2:906097/100‑1 (MQ=255)
                                     tGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCt                                                           <  2:1130736/100‑1 (MQ=255)
                                     tGCCATTGCGTTATGCGAACGCGCAGACGCTCAGTGAGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCt                                                           >  1:567824/1‑100 (MQ=255)
                                       ccATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTg                                                         <  2:1323879/100‑1 (MQ=255)
                                       ccATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTg                                                         >  2:1346569/1‑100 (MQ=255)
                                        cATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGc                                                        <  1:784470/100‑1 (MQ=255)
                                        cATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGc                                                        <  1:35902/100‑1 (MQ=255)
                                         aTTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCg                                                       <  2:1137740/100‑1 (MQ=255)
                                          ttGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCgg                                                      <  2:105542/100‑1 (MQ=255)
                                           tGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGc                                                     <  1:1255642/100‑1 (MQ=255)
                                              gTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTg                                                  <  1:1492069/100‑1 (MQ=255)
                                               ttATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGa                                                 >  1:137308/1‑100 (MQ=255)
                                                    cGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGtt                                            >  1:424166/1‑100 (MQ=255)
                                                      aTCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTa                                          <  2:1776086/100‑1 (MQ=255)
                                                      aTCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTa                                          <  2:101508/100‑1 (MQ=255)
                                                      aTCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTa                                          >  1:783725/1‑100 (MQ=255)
                                                       tCGCGGAGACGCACAGTGAGAACATGATTAAGGGCGAACAGAAACTGCGCGGCTTCTTGCCGGTGACGGACAAGTAACATCTTGCGGCCTGAAAGta                                            <  1:265801/97‑2 (MQ=255)
                                                        cgcgCAGACGCTCAGTGCGAACATCATTGAGTGCGAACAGCACCTGCGCGGCGTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAg                                        <  1:781476/100‑1 (MQ=255)
                                                        cgcgCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAg                                        <  1:1098278/100‑1 (MQ=255)
                                                        cgcgCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAg                                        >  1:1422103/1‑100 (MQ=255)
                                                        cgcgCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAg                                        <  2:1875036/100‑1 (MQ=255)
                                                           gcAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAgtgt                                     <  2:322829/100‑1 (MQ=255)
                                                            cAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGtgtg                                    >  1:723535/1‑100 (MQ=255)
                                                              gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGt                                  <  2:272823/100‑1 (MQ=255)
                                                              gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGt                                  <  2:1824381/100‑1 (MQ=255)
                                                              gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGt                                  <  1:1224834/100‑1 (MQ=255)
                                                              gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGt                                  <  1:1207431/100‑1 (MQ=255)
                                                              gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGt                                  <  1:535932/100‑1 (MQ=255)
                                                              gACGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGt                                  >  1:743918/1‑100 (MQ=255)
                                                               aCGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTa                                 <  1:1018773/100‑1 (MQ=255)
                                                                cGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTAc                                >  1:154564/1‑100 (MQ=255)
                                                                cGCTCAGTGCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTAc                                <  2:54459/100‑1 (MQ=255)
                                                                    cAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATaa                            >  2:1530544/1‑100 (MQ=255)
                                                                        gCGAACATCATTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTgc                        >  2:903662/1‑100 (MQ=255)
                                                                           aaCATCATTAAGTGCGAACAGCACCTGCGCGGCATCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTgcgct                     <  1:475700/100‑2 (MQ=255)
                                                                            aCATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCAt                    >  2:1778242/1‑100 (MQ=255)
                                                                               tcatTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAAGATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCAc                 <  1:1079557/100‑1 (MQ=255)
                                                                               tcatTAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCAc                 <  2:1547285/100‑1 (MQ=255)
                                                                                  ttAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTgg              <  2:1153560/100‑1 (MQ=255)
                                                                                   tAAGTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGc             >  2:1131093/1‑100 (MQ=255)
                                                                                      gTGCGAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGtt          <  2:1389097/100‑1 (MQ=255)
                                                                                          gAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTc      <  2:1217028/100‑1 (MQ=255)
                                                                                          gAACAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTc      <  2:1214858/100‑1 (MQ=255)
                                                                                           aaCAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCt     >  2:398181/1‑100 (MQ=255)
                                                                                            aCAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTg    <  1:1923483/100‑1 (MQ=255)
                                                                                            aCAGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTg    <  1:819002/100‑1 (MQ=255)
                                                                                              aGCACCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGTa  <  2:245057/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GACCATGAGATAAAATTGAAGGAGTCAGTGAAATAATTGCCATTGCGTTATGCGATCGCGCAGACGCTCAGTGCGAACATCATTAAGTGCGAACAGCGCCTGCGCGGCCTCTTCCCGCTGACGGACCAGTAACATCTTGCGGCCTGAAAGTTTAAGTGTGGTACATAACTGCGCATCACTGGCGTTTTTCTGTA  >  NC_000913/2593820‑2594013

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: