Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000731 25,669 (A)7→6 coding (322/333 nt) USA300HOU_pUSA300HOUMR0031 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073125,6630A.96.7% 267.8 / ‑0.2 61coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (28/31);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAAT  >  CP000731/25527‑25799
                                                                                                                                         |                                                                                                                                        
tatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaa                                                                                                                                      <  1:35487/141‑1 (MQ=255)
tatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaa                                                                                                                                      <  1:155821/141‑1 (MQ=255)
 atGTTAGTTACACTAGTGACTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTTATTATAATTGTAGAACTATTTCTTATGATGTTTAATTATTATACTACACATaaaaaa                                                                                                                                     <  2:191329/141‑1 (MQ=255)
            acTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTaa                                                                                                                         >  1:26288/1‑141 (MQ=255)
               aGTGCGTAAATTTACAGGTTTGTCACCAAGTAGTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCAATATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAgg                                                                                                                      <  2:141991/141‑1 (MQ=255)
                gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                                                                     <  1:137600/141‑1 (MQ=255)
                gTGCTTAAAGTTAAAGGTTTGTTACCAAGTATTGAAGGATCTTGGTTTTCTAAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGt                                                                                                                     <  2:274468/141‑1 (MQ=255)
                    ttAAATTTACAGGTTTTTCAACAAGTATTGAAGAATCTTCGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTaaa                                                                                                                 <  1:222639/141‑1 (MQ=255)
                      aaaTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGTTTTTTAATTACTGTACTACAAATAAAAAATATACATCTAAAGGTTAAAtt                                                                                                               >  1:65673/1‑141 (MQ=255)
                       aaTTTACAGGTTTGTCACCAAGTATTGCAGTATACTGATTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTa                                                                                                              >  2:130982/1‑141 (MQ=255)
                       aaTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTa                                                                                                              >  2:14356/1‑141 (MQ=255)
                         tttACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAAt                                                                                                            >  1:88875/1‑141 (MQ=255)
                         tttACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATCGTTCTCTGGTGACAAAAGCTGTTTGCAACTTTTACTTCATATAAAAAATATACATCTAAAGGTTAAATTAAt                                                                                                            <  1:135926/141‑1 (MQ=255)
                            aCAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaa                                                                                                         >  2:16236/1‑141 (MQ=255)
                             cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGaaa                                                                                                        <  1:242350/141‑1 (MQ=255)
                                gTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                                                                     <  1:105189/141‑1 (MQ=255)
                                     tCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAt                                                                                                >  2:224583/1‑141 (MQ=255)
                                     tCACAAAATATTGAAGGATCTTGGTTTTCTCAATATACAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAt                                                                                                <  2:93793/141‑1 (MQ=255)
                                         cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagg                                                                                            >  1:36678/1‑140 (MQ=255)
                                                tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAATAATAGAATTAGAGAAtata                                                                                     >  1:159667/1‑141 (MQ=255)
                                                  aaGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGa                                                                                   <  2:146157/141‑1 (MQ=255)
                                                     gATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAACAATAGAATTGGAGAATATAGaaaa                                                                                >  1:78041/1‑141 (MQ=255)
                                                       tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                                                              <  1:16236/141‑1 (MQ=255)
                                                         ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                                                                            <  2:209578/141‑1 (MQ=255)
                                                         ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                                                                            >  2:274800/1‑141 (MQ=255)
                                                         ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                                                                            >  1:77776/1‑141 (MQ=255)
                                                          tGGTGTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                                                                           <  1:37171/141‑1 (MQ=255)
                                                               ttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGaa                                                                      >  2:88670/1‑141 (MQ=255)
                                                                tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAt                                                                     >  1:29904/1‑141 (MQ=255)
                                                                  tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAttt                                                                   >  1:159175/1‑141 (MQ=255)
                                                                  tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAttt                                                                   <  1:232886/141‑1 (MQ=255)
                                                                   cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATTCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTATAGAATATAGAAAACAAAAAGAAATTtc                                                                  >  2:223187/1‑141 (MQ=255)
                                                                      tatCCAGAAATTTTCTTTATAATTGTATCAATATTTCTTATGCTGTTTAATTACTATACTAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTTGAGAATATAGAAAACAAAAAGGAATTTCTCa                                                               <  2:194737/141‑1 (MQ=255)
                                                                      tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTACAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                                                               <  1:130982/141‑1 (MQ=255)
                                                                      tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGCGCATATAGAAAACATAAAGGAATTTCTCa                                                               <  1:124457/141‑1 (MQ=255)
                                                                      tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                                                               <  1:14356/141‑1 (MQ=255)
                                                                      tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                                                               >  1:199860/1‑141 (MQ=255)
                                                                      tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                                                               >  2:52016/1‑141 (MQ=255)
                                                                            aaaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAc                                                         <  2:233149/140‑1 (MQ=255)
                                                                              aatttaatttATAATAGTAGCAAAAATTATTATTATGTTTAATTAATAAAATAAAAAAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt                                                       <  2:137086/141‑1 (MQ=255)
                                                                                    cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGTATATGGAACAGAAAAAGGAATTTTTCAAGAAACTACAAGTa                                                 >  1:93051/1‑141 (MQ=255)
                                                                                    cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTa                                                 >  2:232022/1‑141 (MQ=255)
                                                                                     tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATAGACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTaa                                                <  1:76221/141‑1 (MQ=255)
                                                                                           aTTGTAGCACAATTTCTTATGCTGGTTAATAACTATACTAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAAGAAATTCTCAAGAAACTTTAAGTAAAGAAtt                                          <  1:272244/141‑1 (MQ=255)
                                                                                           aTTGTAACACTATTTCTTATGCTGTTTAATTACTATACTAAAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAAtt                                          <  2:123862/141‑1 (MQ=255)
                                                                                            ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTa                                         <  1:251482/141‑1 (MQ=255)
                                                                                            ttGGAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTa                                         <  1:99074/141‑1 (MQ=255)
                                                                                               tAGCACTATTACTTATGCTGTTTAATTACTATAATACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaaaa                                      <  2:193528/141‑1 (MQ=255)
                                                                                                        tttATTATGCTGTTTAATTAATATAATAAAAAAAAAAAAAATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAg                              <  2:136577/141‑1 (MQ=255)
                                                                                                         ttCATATGCTGTTTAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAgg                             <  2:65673/141‑1 (MQ=255)
                                                                                                                gCTGTTTAATTACTATACTACAAATAAAAAATATACATCTACATGTTCAAATTATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTa                      <  1:238269/141‑1 (MQ=255)
                                                                                                                     ttAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAAt                 >  2:204745/1‑141 (MQ=255)
                                                                                                                      tAATTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCACACTATTACTg                >  1:59116/1‑141 (MQ=255)
                                                                                                                       aaTTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGc               <  1:112267/141‑1 (MQ=255)
                                                                                                                       aaTTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGc               <  1:88670/141‑1 (MQ=255)
                                                                                                                       aaTTACTACAATACAAATAAAAAATATACAACTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGc               <  1:82906/141‑1 (MQ=255)
                                                                                                                        aTTACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCa              >  1:4019/1‑141 (MQ=255)
                                                                                                                         ttACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCaa             >  1:145049/1‑141 (MQ=255)
                                                                                                                         ttACTATACTACAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCaa             >  1:145184/1‑141 (MQ=255)
                                                                                                                          tactatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAAAACCATTAATTGAGAAAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAAt            <  1:283824/141‑1 (MQ=255)
                                                                                                                           aataaaatacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAAAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATa           <  2:62083/134‑1 (MQ=255)
                                                                                                                             tatactaCAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATaaa         >  2:57904/1‑139 (MQ=255)
                                                                                                                                 ctacAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGaaaaa     >  1:215770/1‑140 (MQ=255)
                                                                                                                                    cAAATAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACaat  >  1:37599/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TATGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAAT  >  CP000731/25527‑25799

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: