Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,610,079:1 (T)7→8 100% coding (447/795 nt) USA300HOU_2473 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,610,0721.T57.4% 108.4 / 88.4 54coding (454/795 nt)USA300HOU_2473lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (10/13);  new base T (12/19);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

ACATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATAC  >  CP000730/2609933‑2610205
                                                                                                                                             |                                                                                                                                       
aCATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTACTATCTTTAAAGGTGCCATAGTGTACAAAGAACTTAAAGTTTTCAATctct                                                                                                                                        >  1:340045/1‑141 (MQ=255)
aCATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATctct                                                                                                                                        >  1:704951/1‑141 (MQ=255)
 cATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCtt                                                                                                                                       <  2:578802/141‑1 (MQ=255)
 cATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCtt                                                                                                                                       <  1:681185/141‑1 (MQ=255)
     ataATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCtttttt                                                                                                                                  >  2:98608/1‑141 (MQ=255)
       aaTCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCtttttttt                                                                                                                               >  1:684931/1‑140 (MQ=255)
        atcatcATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCtttttttta                                                                                                                              >  1:747815/1‑139 (MQ=255)
          catcatTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCttttttttaat                                                                                                                            >  2:135609/1‑137 (MQ=255)
             catTATTCAATTGATACTTTGCCGAATAACTTGGTACATTACGATTATATGAAATATCACCATCTTTATAATCCTTAATATCTCTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCttttttttcacttg                                                                                                                         <  2:341887/141‑8 (MQ=14)
              attattCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTc                                                                                                                        <  1:21242/141‑1 (MQ=255)
                  ttCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACACAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCCtt                                                                                                                     >  1:708034/1‑141 (MQ=255)
                   tCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTGTAACTTGTCATTc                                                                                                                    >  1:690855/1‑141 (MQ=255)
                   tCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCg                                                                                                                   >  2:65610/1‑141 (MQ=255)
                   tCAATTGATACTTCGCCGAATAACTTGGTACTTTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTc                                                                                                                    <  1:592367/141‑1 (MQ=17)
                      aTTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTa                                                                                                                <  2:335259/141‑1 (MQ=255)
                        tGAAATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAg                                                                                                               <  2:309517/141‑1 (MQ=255)
                         gATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGt                                                                                                              <  1:597805/141‑1 (MQ=255)
                          atatTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTg                                                                                                             <  1:418870/141‑1 (MQ=255)
                              tttGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTc                                                                                                        <  2:347475/141‑1 (MQ=255)
                               ttGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAATGCACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTc                                                                                                        <  2:284649/141‑1 (MQ=255)
                                    gAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGt                                                                                                  >  1:612961/1‑141 (MQ=255)
                                     aaTAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTa                                                                                                 <  2:287876/141‑1 (MQ=255)
                                      aTAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATGCGGTGGTGGCTTCGTTGGttt                                                                                                >  1:83161/1‑139 (MQ=255)
                                       tAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGAATg                                                                                               >  2:222122/1‑141 (MQ=255)
                                              gTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGATATGAttttt                                                                                         >  1:762843/1‑140 (MQ=255)
                                                  aTTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTATTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATg                                                                                    >  1:257479/1‑141 (MQ=255)
                                                         ttATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCAt                                                                              <  2:240502/141‑1 (MQ=255)
                                                           atatGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCAttt                                                                            <  1:654913/141‑1 (MQ=255)
                                                             atGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCAttttt                                                                          <  2:504651/141‑1 (MQ=255)
                                                                aaaTATCACCATCTTTATGATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTAcc                                                                       <  1:161842/141‑1 (MQ=255)
                                                                      cACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCACTCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATAc                                                                 <  1:682955/141‑1 (MQ=255)
                                                                       aCCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATAcc                                                                >  2:312122/1‑141 (MQ=255)
                                                                               tataATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTtca                                                        <  1:173028/141‑1 (MQ=255)
                                                                                ataATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTtctt                                                       >  2:278060/1‑139 (MQ=255)
                                                                                     taTTAATATCGTTAAAGTTGCCATTTTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATcc                                                  <  2:230068/140‑1 (MQ=255)
                                                                                        tAATATCTTTAAAGTTGCCATTTTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTg                                              <  2:699927/141‑1 (MQ=255)
                                                                                        tAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTg                                              <  2:700324/141‑1 (MQ=255)
                                                                                            atCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTTTTCTTTCTTTGTTCt                                          >  2:737715/1‑141 (MQ=25)
                                                                                            atCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCt                                          >  2:737707/1‑141 (MQ=255)
                                                                                            atCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCt                                          >  2:87205/1‑141 (MQ=255)
                                                                                               tttAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTTTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGc                                        >  2:663236/1‑141 (MQ=255)
                                                                                                 tcaaGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGt                                     <  2:294636/139‑1 (MQ=255)
                                                                                                  aaaGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGt                                     <  1:488764/141‑1 (MQ=255)
                                                                                                    aGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTa                                   >  1:752097/1‑141 (MQ=255)
                                                                                                    aGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTa                                   <  1:242631/141‑1 (MQ=255)
                                                                                                               tGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTGTTCCGTTCCATCTCTACCTTAGACCTTTTCTCATCATCCTTTGGTCTGCCGTTActgtcatctgt                        <  2:166298/141‑1 (MQ=25)
                                                                                                               tGTAAAAAGAACTTAAAGTTTTCAATCTTTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActgtcatctgt                        <  2:106713/141‑1 (MQ=255)
                                                                                                                       gAACTTAAATTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCAc                <  1:98608/141‑1 (MQ=255)
                                                                                                                       gAACTTAAAGTTTTCAATCTCTTTTTTTAACCTCTCAATCCGTAGTTGCTGCCTTGGTATGATTTTATTATGTGCCATTTGTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACt               <  2:542478/141‑1 (MQ=255)
                                                                                                                       aaaCTTAAAGTTTTCAATCTCTTTTTTTAACTTTTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTActgtcatctgccctctctcg               <  2:398918/140‑11 (MQ=255)
                                                                                                                         aCTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTTTTTTTTTTTTGGTCTGCCATTTCTGTCTTCTGTCTTTTCTCtt              >  2:546373/1‑141 (MQ=25)
                                                                                                                         aCTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACtt              <  1:657069/141‑1 (MQ=255)
                                                                                                                         aCTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTAGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCCTTTCACtt              >  2:180571/1‑141 (MQ=255)
                                                                                                                         aCTTAAAGTTTTCAATCTCTTTTTCGTAACTTTTCAATTGGTCGACCAATCCTTGTTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACtt              <  1:569764/141‑1 (MQ=255)
                                                                                                                          cTTAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATAATCCTTTGGTCTGCCGTTACTGTCATCTGTCTTTTCACttt             >  1:657128/1‑140 (MQ=255)
                                                                                                                           ttAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTtat            >  1:270704/1‑141 (MQ=255)
                                                                                                                            tAAAGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTtata           >  1:518591/1‑141 (MQ=255)
                                                                                                                             taaGTTTTCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaa          <  2:67978/140‑1 (MQ=255)
                                                                                                                             aaaGTTTTCAATCTTTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaa          <  2:658932/141‑1 (MQ=255)
                                                                                                                             aaaGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCCTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaaa         >  1:422993/1‑141 (MQ=255)
                                                                                                                              accttttCAACTCCTTTTTTTACCTTGTCATTCGGTAGTGGTCTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCAGCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATaaaa        <  1:222122/138‑1 (MQ=255)
                                                                                                                                   ttCAATCTCTTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCTTTTGATCTGCCGCTACTGTCACCGGTCATTCCACGTCTAAaatata    >  2:55719/1‑139 (MQ=255)
                                                                                                                                    tCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATAc  <  2:612961/141‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                       
ACATTATAATCATCATTATTCAATTGATATTTTGCCGAATAACTTGGTACATTAGGATTATATGAAATATCACCATCTTTATAATCATTAATATCTTTAAAGTTGCCATATTGTACAAAGAACTTAAAGTTTTCAATCTCTTTTTTTAACTTGTCATTCGGTAGTGGCTTCGTTGGTATGATTTTATTATGTTCCATTTTTACCGGATACCTTTTTTCATCATCCTTTGGTCTGCCGTTACTGTCATCTGTCATTTCACTTATAAAATAATAC  >  CP000730/2609933‑2610205

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: