Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,135,804 T→G I190L (ATT→CTT)  pitB ← phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,135,8040TG100.0% 114.2 / NA 40I190L (ATT→CTT) pitBphosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/25);  total (15/25)

TACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTAAACCGAAGAACC  >  NC_000913/3135585‑3136037
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
tacgggtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGt                                                                                                                                                                                                                           <  1:266166/238‑1 (MQ=255)
 acgggtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTcc                                                                                                                                                                                                                         <  1:291556/239‑1 (MQ=255)
                      cGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCgcagc                                                                                                                                                                                                     <  1:90949/238‑1 (MQ=255)
                                                     gCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGc                                                                                                                                                                     <  1:822654/238‑1 (MQ=255)
                                                           cccACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCgcagca                                                                                                                                                                                                    <  1:342127/202‑1 (MQ=255)
                                                              accaGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTcc                                                                                                                                                                                                                         <  1:105032/178‑1 (MQ=255)
                                                                     ccAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTTTGGTCTTTCGGAAAGCGGGGAAGGCGGGCCCGCTTTTTTGTCCCGCTCCCGTAGCGTCGCCGCAGGAAAATCAg                                                                                                                                                                                         >  2:342127/1‑202 (MQ=37)
                                                                       aGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTcc                                                                                                                                                                                                                         <  1:111986/169‑1 (MQ=255)
                                                                       aGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTGCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTcc                                                                                                                                                                                                                         <  1:920441/169‑1 (MQ=255)
                                                                        gCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTTCGATCTTCCGGAAAGCGGGGGAGACGGTCACGCTTTTTTTTCCCGCTCCCGTa                                                                                                                                                                                                               >  2:105032/1‑178 (MQ=37)
                                                                        gCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAg                                                                                                                                                                                         <  1:636143/199‑1 (MQ=255)
                                                                                 ggCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTCGTTTGCCTTTTTTTTTTTTTCGATCTTTCCGAATGCGGGGAAGACGGTCACCCCTTTTTCTCCCGCTCCCGtc                                                                                                                                                                                                               >  2:920441/1‑168 (MQ=25)
                                                                                 ggCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGTCTTTTTTTTTTTTGCGATCTTTCGGGATGCGGGGAAGAAGGTCACGCTTTTTCTTTCCCCTcccgaa                                                                                                                                                                                                               >  2:111986/1‑165 (MQ=37)
                                                                                  gCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGCTTTCGGTTGACCTTTTTTTTTTTTCGGTCCTCCGGAATGCGGGGAAGAAGGGCACGCCTTTTTGTTCCGCCCCAGTAGCCGCGCACCAGGGAAAACAGGGCCCCCCGCa                                                                                                                                                                               >  2:636143/1‑200 (MQ=19)
                                                                                      aTCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCACTACGCGTCCAGAATGGCGGGTTCTGTTTGCCCTTTTTTTTTTTGCGATTTTCCGGAATGCGGTGAAGGCGGTCACGCTTTTTTGTTCCGCCCCCCTAGCGTCGCAGCAAGGAGACACGGCCTCCCGCCACGACCACGCCGCCGCGCAGGGGAACAATCACCGGCGAGaaa                                                                                                                                    >  2:1046856/1‑237 (MQ=19)
                                                                                      aTCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTCGTTTGCCTTTTTTCTTTTTTTGATCTTCCGGAATGCGGGGAAGACGGGCACGCTTTTTTGTCCCCCTCCCGGAGCGTCGCCGCAGGAAAAATAGGGCCCCCGCCACGACCCAGCCGGCCAGAGGGGGAAAAAAAAGCGGGGGGaaaaa                                                                                                                                  >  2:290234/1‑238 (MQ=21)
                                                                                      aTCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTGTGATCTTCCGGAATGCGGGGAAGACGGGCACGCTTTTTTGTTCCGCTCCAGTAGCGTCGCAGCAGGCATATAAGGACCCCCGCCATGACCAAGCCGTCGAGAGGGGGAAACAACATGGAGGAGaaaaa                                                                                                                                  >  2:305180/1‑238 (MQ=37)
                                                                                       tCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGATTACGGTTTCCTTTTTTTTTTTTTGGATCTTTCGGAATGCGGGGAAGAAGGGCAAGATTTTTTGTCCCGCTCCAGTAGCGGCGCGCCAGGAACATTAGGCCTCCCGCCAAGAACACGCCGACGATAGGGGGGAAAATAAGCGAGGAGAAAAtt                                                                                                                                 >  2:909491/1‑239 (MQ=25)
                                                                                            ttttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGt                                                                                                                            <  1:1119857/239‑1 (MQ=255)
                                                                                                     tcgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACg                                                                                                                   <  1:410707/239‑1 (MQ=255)
                                                                                                      cgtTCGCGGCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACg                                                                                                                   <  2:679035/238‑1 (MQ=255)
                                                                                                               cgggCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAAc                                                                                                         <  1:558873/236‑1 (MQ=255)
                                                                                                               cGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAAc                                                                                                         <  1:290234/239‑1 (MQ=255)
                                                                                                                    gAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGcatc                                                                                                    <  2:550571/239‑1 (MQ=255)
                                                                                                                                       cGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTcc                                                                                                                                                                                                                         <  2:580990/105‑1 (MQ=255)
                                                                                                                                          aaaCAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTaa                                                                               <  1:305180/238‑1 (MQ=255)
                                                                                                                                                 cAGCGCAATACGCGTCCAGAATGGCGGTTTTCGGTTTCCTTTTTTTTTTTTTCGGTCTTCCGGGATGCGGGGAAGCCGGGCACGCTTTTTCTTCCCGCTCCAGTa                                                                                                                                                                                                               >  1:580990/1‑105 (MQ=255)
                                                                                                                                                 cAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTTCCTTTTTTTTTTTTGCGATCTTCCCGGATGCGGTGAAGACGGTCACGCTTTTTTTTTCCGCCCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCCCAATGACCCGGCCGACGAAAGGGGAAACAATCAGCGAGGGGAAAATTTTGGTCCCTTCACGCAGGGTTTACGCATCCATCACCGAACAGCCGGTTAACAGCGCGt                                                                       >  2:83881/1‑239 (MQ=255)
                                                                                                                                                            gcgTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCg                                                            <  1:909491/239‑1 (MQ=255)
                                                                                                                                                                     gtGGCGGTTTACGTTTGCCTTTTTTCTTTGTGCGATCTTCCGGAATGCGGGGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTAGCAGCATGAATATCAGGCCTCCCGCAATGACCAGGACGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCAGCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGcc                                                       <  2:313607/234‑1 (MQ=255)
                                                                                                                                                                             ttACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCAc                                            <  1:543233/238‑1 (MQ=255)
                                                                                                                                                                                      ccTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTgg                                                                                                                             >  1:175102/1‑148 (MQ=255)
                                                                                                                                                                                      ccTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTgg                                                                                                                             <  2:175102/148‑1 (MQ=255)
                                                                                                                                                                                       cTTTTTTCTTTTTGCGATTTTCCGGAAGGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGgt                                  <  2:902809/238‑1 (MQ=255)
                                                                                                                                                                                           tttCTTTTTGCGATCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCg                             <  2:1110031/239‑1 (MQ=255)
                                                                                                                                                                                                        tCTTCCGGAATGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCGCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGAAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCg                 <  1:909567/238‑1 (MQ=255)
                                                                                                                                                                                                                  tGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGcgcg                                                                        <  2:23902/173‑1 (MQ=255)
                                                                                                                                                                                                                  tGCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGcgcg                                                                        >  1:23902/1‑173 (MQ=255)
                                                                                                                                                                                                                   gCGGTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTAAACCgaaga     >  2:636398/1‑239 (MQ=255)
                                                                                                                                                                                                                      gTGAAGACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTAAACCGAAGAAcc  >  1:407149/1‑239 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
TACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTAAACCGAAGAACC  >  NC_000913/3135585‑3136037

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