Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,421,700 T→A 5.8% intergenic (+136/+4) pepF1_2 → / ← phoU Oligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,421,7000TA 5.8% 300.5 / 8.0 89intergenic (+136/+4)pepF1_2/phoUOligoendopeptidase F, plasmid/Phosphate‑specific transport system accessory protein PhoU
Reads supporting (aligned to +/- strand):  ref base T (28/54);  new base A (2/3);  total (32/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTAAATGTCT  >  USA300TCH1516_ALE/1421613‑1421790
                                                                                       |                                                                                          
tACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAg                                                                                         <  1:2036180/91‑1 (MQ=255)
 aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt                                                                                        <  1:1229476/91‑1 (MQ=255)
 aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt                                                                                        <  1:3802755/91‑1 (MQ=255)
 aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt                                                                                        <  1:359269/91‑1 (MQ=255)
 aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt                                                                                        <  1:303109/91‑1 (MQ=255)
 aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGt                                                                                        <  1:302587/91‑1 (MQ=255)
 aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTAAAGt                                                                                        <  1:1302510/91‑1 (MQ=255)
 aCTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTAAAGt                                                                                        <  1:280861/91‑1 (MQ=255)
    aaTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGATTAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTAAAGTTAt                                                                                     <  1:1234774/91‑1 (MQ=255)
    aaTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAt                                                                                     >  1:393034/1‑91 (MQ=255)
     aTTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTAtt                                                                                    <  1:1382678/91‑1 (MQ=255)
      ttGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTg                                                                                   <  1:417286/91‑1 (MQ=255)
         tgAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTc                                                                                <  1:2186216/91‑1 (MQ=255)
         tgAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTc                                                                                <  1:2186710/91‑1 (MQ=255)
           agaATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGt                                                                              <  1:3028835/91‑1 (MQ=255)
           agaATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGt                                                                              <  1:1714619/91‑1 (MQ=255)
            gaATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTa                                                                             >  1:4122263/1‑91 (MQ=255)
                cTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAAtgt                                                                         <  1:3231001/91‑1 (MQ=255)
                cTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAAtgt                                                                         <  1:3231378/91‑1 (MQ=255)
                 ttGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAAtgtg                                                                        >  1:2180232/1‑91 (MQ=255)
                 ttGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAAtgtg                                                                        >  1:165993/1‑91 (MQ=255)
                 ttGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAAtgtg                                                                        >  1:3824567/1‑91 (MQ=255)
                  tGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATgtgt                                                                       >  1:1952930/1‑91 (MQ=255)
                   ggTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAATGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTa                                                                      >  1:4136578/1‑91 (MQ=255)
                    gTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGAGTAc                                                                     >  1:3461496/1‑91 (MQ=255)
                       tttttGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTg                                                                  <  1:135039/91‑1 (MQ=255)
                        ttttGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGt                                                                 <  1:260213/91‑1 (MQ=255)
                         tttGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGtt                                                                <  1:2672776/91‑1 (MQ=255)
                         tttGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGtt                                                                >  1:2062328/1‑91 (MQ=255)
                          ttGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTta                                                               <  1:681940/91‑1 (MQ=255)
                          ttGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTta                                                               >  1:79203/1‑91 (MQ=255)
                           tGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTtaa                                                              <  1:3279007/91‑1 (MQ=255)
                           tGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTtaa                                                              >  1:3055270/1‑91 (MQ=255)
                            gAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTtaaa                                                             <  1:1852474/91‑1 (MQ=255)
                             aGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTtaaat                                                            <  1:559786/91‑1 (MQ=255)
                               ggTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTtaaataa                                                          >  1:2647308/1‑91 (MQ=255)
                               ggTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTtaaataa                                                          >  1:2646980/1‑91 (MQ=255)
                                  ttATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTaaataaaaa                                                       <  1:2120063/91‑1 (MQ=255)
                                   tATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTaaataaaaat                                                      >  1:1050065/1‑91 (MQ=255)
                                    aTGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTaaataaaaata                                                     <  1:1604186/91‑1 (MQ=255)
                                     tGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTaaataaaaataa                                                    <  1:634046/91‑1 (MQ=255)
                                      gAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAaataaaaataaa                                                   <  1:1904893/91‑1 (MQ=255)
                                      gAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAaataaaaataaa                                                   <  1:3741696/91‑1 (MQ=255)
                                         taataGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAacac                                                <  1:1888465/91‑1 (MQ=255)
                                         taataGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAacac                                                >  1:3939963/1‑91 (MQ=255)
                                         taataGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAacac                                                >  1:1251949/1‑91 (MQ=255)
                                          aataGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACt                                               >  1:2039838/1‑91 (MQ=255)
                                             aGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTc                                            <  1:3108731/91‑1 (MQ=255)
                                                 atatatatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCg                                        <  1:373491/91‑1 (MQ=255)
                                                  tatatatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGa                                       <  1:1503891/91‑1 (MQ=255)
                                                  tatatatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGa                                       >  1:1669263/1‑91 (MQ=255)
                                                   atatatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGAt                                      <  1:201243/91‑1 (MQ=255)
                                                   atatatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGAt                                      <  1:1758925/91‑1 (MQ=255)
                                                   atatatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGAt                                      >  1:2743550/1‑91 (MQ=255)
                                                    tatatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATg                                     >  1:3280180/1‑91 (MQ=255)
                                                     atatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGt                                    >  1:217975/1‑91 (MQ=255)
                                                     atatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGt                                    <  1:1308024/91‑1 (MQ=255)
                                                      tatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGtt                                   >  1:1721899/1‑91 (MQ=255)
                                                      tatTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGtt                                   <  1:2875104/91‑1 (MQ=255)
                                                       atTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTta                                  <  1:2391363/91‑1 (MQ=255)
                                                        tttATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTtaa                                 <  1:1694550/91‑1 (MQ=255)
                                                        tttATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTtaa                                 <  1:1796046/91‑1 (MQ=255)
                                                           aTCGTAAGGGATTTTATAGTAATAGCTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTtaataa                              <  1:73219/91‑1 (MQ=255)
                                                           aTCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTtaataa                              <  1:1527459/91‑1 (MQ=255)
                                                           aTCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTtaataa                              <  1:673379/91‑1 (MQ=255)
                                                           aTCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTtaataa                              <  1:73204/91‑1 (MQ=255)
                                                           aTCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTtaataa                              <  1:3834097/91‑1 (MQ=255)
                                                                  gggATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATc                       <  1:4053015/91‑1 (MQ=255)
                                                                  gggATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATc                       <  1:808834/91‑1 (MQ=255)
                                                                  gggATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATc                       >  1:3203757/1‑91 (MQ=255)
                                                                  gggATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATc                       >  1:1139340/1‑91 (MQ=255)
                                                                  gggATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATc                       <  1:808583/91‑1 (MQ=255)
                                                                    gATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCAc                     <  1:1522095/91‑1 (MQ=255)
                                                                      ttttATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCa                   <  1:1108422/91‑1 (MQ=255)
                                                                      ttttATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCa                   <  1:2424596/91‑1 (MQ=255)
                                                                        ttATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAAt                 <  1:3407572/91‑1 (MQ=255)
                                                                        ttATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAAt                 <  1:3017073/91‑1 (MQ=255)
                                                                        ttATAGTAATAACTTGAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAAt                 >  1:3538937/1‑91 (MQ=255)
                                                                        ttATAGTAATAACTTCAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAAt                 >  1:3538924/1‑91 (MQ=255)
                                                                         tataGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATa                <  1:3365276/91‑1 (MQ=255)
                                                                          ataGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATAc               <  1:3895096/91‑1 (MQ=255)
                                                                          ataGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATAc               <  1:346341/91‑1 (MQ=255)
                                                                              taataaCTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTc           >  1:2123466/1‑91 (MQ=255)
                                                                              taataaCTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTc           >  1:303364/1‑91 (MQ=255)
                                                                                ataaCTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTa         >  1:695601/1‑91 (MQ=255)
                                                                                ataaCTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTa         <  1:2445807/91‑1 (MQ=255)
                                                                                 taaCTTAAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTaa        >  1:2669888/1‑91 (MQ=255)
                                                                                   aCTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTAAAt      >  1:695911/1‑91 (MQ=255)
                                                                                   aCTTAAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTAAAt      >  1:886786/1‑91 (MQ=255)
                                                                                     tttAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTAAATGt    <  1:1529158/91‑1 (MQ=255)
                                                                                      ttAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTAAATGTc   >  1:3318045/1‑91 (MQ=255)
                                                                                       tAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTAAATGTCt  >  1:4032414/1‑91 (MQ=255)
                                                                                       |                                                                                          
TACTAATTGTGAGAATCTTGGTTTTTTTGAGGGTTTATGAATAATAGAAATATATATTTATCGTAAGGGATTTTATAGTAATAACTTTAAGTTATTGTTCGTAATGTGTACCTGTTAAATAAAAATAAACACTTTCAGCGATGTTAATAATATGATCACCAATACGTTCTAAATGTCT  >  USA300TCH1516_ALE/1421613‑1421790

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