Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,436,416 A→G 10.2% L28P (CTT→CCT)  USA300HOU_2299 ← APC family amino acid‑polyamine‑organocation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,436,4160AG10.2% 137.9 / 9.9 61L28P (CTT→CCT) USA300HOU_2299APC family amino acid‑polyamine‑organocation transporter
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base G (3/3);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-01

GATCCTGTTGTAGGGTTTAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACTAAATAAATCCATAAAAAAACTTATTGCCG  >  CP000730/2436281‑2436556
                                                                                                                                       |                                                                                                                                            
gATCCTGTTGTAGGGTTTAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCACACCCATGAAAAGcccc                                                                                                                                         >  2:338683/1‑141 (MQ=255)
 aTCCTGTTGTAGGGTTTAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGccccc                                                                                                                                        >  2:303689/1‑140 (MQ=255)
  tCCTGTTGTAGGGTTTAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCAc                                                                                                                                       <  1:205386/141‑1 (MQ=255)
   ccTGATGTAGGGTTTAAATAAATCATTTCCCACATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCAGTCCATTTAATTGTACTTGTAGCACCCATGAAATGCCCCACg                                                                                                                                      >  2:91824/1‑141 (MQ=255)
     tgttgtAGGGTATAAATAAATCATTTCCCCCATTGCTCTAATGATTAAAAATAAAAAGATAGCCTCAATTAAATATGCAAGGATAACTGATGTGCCTGTCCATTTAATTGTACTTGTCGAACCCATGAAAAGCCCCACGcc                                                                                                                                    >  2:597884/1‑141 (MQ=255)
        tgtAGGGTTTAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAAt                                                                                                                                 <  2:530777/141‑1 (MQ=255)
        tgtAGGGTTTAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAAt                                                                                                                                 <  2:695165/141‑1 (MQ=255)
         gtAGGGTTTAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAACGCCCCGCGCCCATc                                                                                                                                >  1:224109/1‑141 (MQ=255)
               tttAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTaccc                                                                                                                          >  1:422294/1‑140 (MQ=255)
                ttAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTacccc                                                                                                                         >  2:52526/1‑139 (MQ=255)
                ttAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTaccac                                                                                                                         >  2:52471/1‑141 (MQ=255)
                ttAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTaccac                                                                                                                         >  1:390559/1‑141 (MQ=255)
                     taaatCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGt                                                                                                                    <  1:242275/141‑1 (MQ=255)
                      aaatCATTCTCCCCACTGCTTTCATGATTAAAAATAAAAATATCCCCGCAATTAAATCTTCACGCATAACTGAGGCTCCTTCGCATTTCATTGCATTTGTCGCACCCATGAACAGCCCCACGCCAATCGTACCACCAAGTg                                                                                                                   <  2:16609/141‑1 (MQ=255)
                         tCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCaa                                                                                                                <  1:400391/141‑1 (MQ=255)
                             ttCCCCCATTGCTCTCATGATTAAAAATACAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCat                                                                                                            >  2:328290/1‑141 (MQ=255)
                              tCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCCCCAAGTGCAATCatt                                                                                                           >  1:665462/1‑141 (MQ=255)
                               cccccATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTACATTTAATTGTACTTGTCGCACCAATGAAAAGCCCCACGACAATCGTACTACCAAGTGCAATCattt                                                                                                          >  1:50796/1‑141 (MQ=255)
                                     ttGCTCTCAAGAATAAAAAGAAAAAGCTACCCCAAATTAACGATGCAAGGATAACTGATGGGCCTGTCCACTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCatttgaattt                                                                                                    <  1:39118/141‑1 (MQ=255)
                                       gCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGTAAATCCTTGCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCatttgaatttga                                                                                                  <  1:774465/141‑1 (MQ=255)
                                              tGATTAAAAATAAAAAGCTACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCAc                                                                                           <  1:508861/141‑1 (MQ=255)
                                                  taaaaataaaaaGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGCCGGGCACTCaa                                                                                       >  1:557877/1‑141 (MQ=255)
                                                  aaaaaataaaaaGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCAACCATGAAAAGCCCCACGCCAATCGCACCACCAAGTGCAATCAGTTGAATGTTACGGGCACTCaa                                                                                       <  2:545077/140‑1 (MQ=255)
                                                   aaaaataaaaaGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAAt                                                                                      >  1:466591/1‑141 (MQ=255)
                                                   aaaaataaaaaGATACCCGAAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAAt                                                                                      >  1:419695/1‑141 (MQ=255)
                                                        taaaaaGATACCCGCAATTAAATATGCAAGGATAACTGAAGGGCCTGTCCATGTAATTGTAATTGTCGAACCCATGAAAAGCCCCAAGCAAATCGAACCACAAATAGCGAACATTTAAATTTCTAGGACACTCCATcccat                                                                                 >  2:41257/1‑139 (MQ=255)
                                                         aaaaaGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCtt                                                                                <  2:373748/141‑1 (MQ=255)
                                                             aGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTa                                                                            >  1:467972/1‑141 (MQ=255)
                                                                   ccGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCAt                                                                      >  1:344291/1‑141 (MQ=255)
                                                                    cGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAGGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCAtt                                                                     >  2:323884/1‑141 (MQ=255)
                                                                    cGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAGGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCAtt                                                                     >  2:671733/1‑141 (MQ=255)
                                                                    cGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAGGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCAtt                                                                     >  2:323347/1‑141 (MQ=255)
                                                                    cGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCAtt                                                                     >  1:464112/1‑141 (MQ=255)
                                                                     gCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCAttt                                                                    >  1:201050/1‑140 (MQ=255)
                                                                     gCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCAttt                                                                    >  1:623664/1‑140 (MQ=255)
                                                                           aaaTATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCCCGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTcc                                                              <  1:342660/141‑1 (MQ=255)
                                                                             agatGCACGGATCCCTGCTGGGCCTGTCTATTTAATTGGACTTTTCGCACCCATTAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCAt                                                            <  2:743450/139‑1 (MQ=255)
                                                                                    aGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAActtct                                                     <  1:594508/141‑1 (MQ=255)
                                                                                        tAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTAcc                                                 >  1:201407/1‑141 (MQ=255)
                                                                                        tAACGGATGGGCGTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTGAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTAcc                                                 >  1:557308/1‑141 (MQ=255)
                                                                                          aCTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTAcccc                                               >  2:493051/1‑141 (MQ=255)
                                                                                          aCTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTAcccc                                               >  1:389879/1‑141 (MQ=255)
                                                                                           cTGATGCGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCa                                              <  1:564262/141‑1 (MQ=255)
                                                                                             taTGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTCACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTTCCCCAtc                                            <  2:336165/140‑1 (MQ=255)
                                                                                             gaggggCCTGTCCAATTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACTGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCAtc                                            <  2:509712/138‑1 (MQ=255)
                                                                                             gATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCACAtc                                            >  1:751854/1‑141 (MQ=255)
                                                                                               tGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCAtctc                                          >  1:769420/1‑141 (MQ=255)
                                                                                                  gCCTGTCCATTTAATTGTACTTGTCGCACCCATGAACAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATATGATTCATATTTATGCCGCCTTGTAGTCCATTGTTCTCCATACCTTCTTAACCCTTCTCtga                                       <  2:556325/141‑1 (MQ=255)
                                                                                                    cTGTCCATTTAATTGTACTTGTCGCACCCATGAAAGGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCtgatt                                     <  1:323347/141‑1 (MQ=255)
                                                                                                    cTGTCCATTTAATTGTACTTGTCGCACCCATGAAAGGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCtgatt                                     <  1:671733/141‑1 (MQ=255)
                                                                                                    cTGTCCATTTAATTGTACTTGTCGCACCCATGAAAGGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCtgatt                                     <  1:323884/141‑1 (MQ=255)
                                                                                                       tCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCtgattgat                                  >  1:261314/1‑141 (MQ=255)
                                                                                                         cATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATAc                                <  1:303689/141‑1 (MQ=255)
                                                                                                         cATTTAAGTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATAc                                <  1:528780/141‑1 (MQ=255)
                                                                                                           tttAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACta                              >  1:744603/1‑141 (MQ=255)
                                                                                                            ttAATTGTACTTGTCACACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTCACCCCATCTCCGATTGATACtaa                             <  2:383339/141‑1 (MQ=255)
                                                                                                                    aattGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCATGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACTAAATAAATcc                     <  2:523734/139‑1 (MQ=255)
                                                                                                                     cTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACTAAATAAATCCa                    <  1:328290/141‑1 (MQ=255)
                                                                                                                          cGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACTAAATAAATCCATaaaa               >  1:327834/1‑141 (MQ=255)
                                                                                                                             aCCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACTAAATAAATCCATaaaaaaa            <  1:757941/141‑1 (MQ=255)
                                                                                                                               tcATGAAAAGCCCAACGCCAATCGTACCACCAAGTGCAATCATATGAATTTGACGGGCACTCAATACCCTTTGTGGTTCATTGTTCTGCATAACTTATTACCCCAACTCCGATTGCGACGAATTCAATCCATAAAAAAACt          >  2:161788/2‑141 (MQ=255)
                                                                                                                                     aaaGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACTAAATAAATCCATAAAAAAACTTATTGc    >  1:542562/1‑141 (MQ=255)
                                                                                                                                       aGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACTAAATAAATCCATAAAAAAACTTATTGCCg  <  2:654035/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
GATCCTGTTGTAGGGTTTAAATAAATCATTTCCCCCATTGCTCTCATGATTAAAAATAAAAAGATACCCGCAATTAAATATGCAAGGATAACTGATGGGCCTGTCCATTTAATTGTACTTGTCGCACCCATGAAAAGCCCCACGCCAATCGTACCACCAAGTGCAATCATTTGAATTTGACGGGCACTCAATCCCCTTTGTAGTTCATTGTTCTCCATAACTTCTTACCCCATCTCTGATTGATACTAAATAAATCCATAAAAAAACTTATTGCCG  >  CP000730/2436281‑2436556

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