Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,430,691:1 +A 8.2% intergenic (‑115/‑270) suhB_2 ← / → birA_2 Inositol‑1‑monophosphatase/Bifunctional ligase/repressor BirA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,430,6911.A 8.2% 241.1 / 15.1 61intergenic (‑115/‑270)suhB_2/birA_2Inositol‑1‑monophosphatase/Bifunctional ligase/repressor BirA
Reads supporting (aligned to +/- strand):  ref base . (26/30);  new base A (2/3);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

TATTTGCTGTAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAATATCATAGCA  >  USA300TCH1516_ALE/2430550‑2430829
                                                                                                                                              |                                                                                                                                          
tATTTGCTGTAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTg                                                                                                                                             <  1:544299/141‑1 (MQ=255)
  tttGCTGTAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGc                                                                                                                                           <  1:388079/141‑1 (MQ=255)
    tGCTGTAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTa                                                                                                                                         >  1:630372/1‑141 (MQ=255)
       tGTAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGcc                                                                                                                                      <  1:694447/141‑1 (MQ=255)
        gTAATGTTTTATCTGTCATTTCTTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCg                                                                                                                                     <  1:764980/141‑1 (MQ=255)
        gTAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCg                                                                                                                                     <  1:93007/141‑1 (MQ=255)
         tAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTGACTAGCCGa                                                                                                                                    <  1:290454/141‑1 (MQ=255)
         tAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTGACTAGCCGa                                                                                                                                    <  1:528775/141‑1 (MQ=255)
          aaTGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGCACGTGACTAGCCGat                                                                                                                                   <  2:147208/141‑2 (MQ=255)
             gTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAATGATCGTCTAGTCACGTGCTAGCCGACaaa                                                                                                                                <  2:75900/141‑1 (MQ=255)
              ttttATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAAt                                                                                                                               <  1:307339/141‑1 (MQ=255)
              ttttATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAAt                                                                                                                               >  1:151840/1‑141 (MQ=255)
               tttATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATa                                                                                                                              <  2:44455/141‑1 (MQ=255)
               tttATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATa                                                                                                                              >  1:343426/1‑141 (MQ=255)
               tttATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACACTGTTTTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATa                                                                                                                              <  1:445988/141‑1 (MQ=255)
                   tCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaa                                                                                                                          <  1:516021/141‑1 (MQ=255)
                      gTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAACAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaagga                                                                                                                       >  2:316509/1‑141 (MQ=255)
                      gTCATTTCGTCCACCACTTCTCATATAATTATCATTTTAATATGACTCTATAGTAAAAGAATCAACAATACAACTGATGACTTCTTCATTTTATGCATAAAAAAAACGGCTAGTCAAGTGCTAGCCGACAAATAgaaagga                                                                                                                       <  1:705089/141‑1 (MQ=255)
                       tCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaaggaa                                                                                                                      >  1:774787/1‑141 (MQ=255)
                       tCATTTCGTCCACCACTTCTCATATCATTATCATTTTATGATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAgaaaggaa                                                                                                                      >  1:21637/1‑141 (MQ=255)
                          tttCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAaagt                                                                                                                   <  1:625666/141‑1 (MQ=255)
                            tCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAaagtaa                                                                                                                 >  1:220615/1‑141 (MQ=255)
                             cGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAaagtaag                                                                                                                >  1:297423/1‑141 (MQ=255)
                                 caccacTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGtaat                                                                                                            <  1:34296/141‑1 (MQ=255)
                                        tctcATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATAtt                                                                                                     <  2:261233/141‑1 (MQ=255)
                                           cATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGACAGTAAGTAATAAATATTGaa                                                                                                  >  1:32856/1‑141 (MQ=255)
                                           cATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGaa                                                                                                  <  1:603887/141‑1 (MQ=255)
                                            atatCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAg                                                                                                 <  1:492549/141‑1 (MQ=255)
                                                 attatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAtgtt                                                                                            <  2:313051/141‑1 (MQ=255)
                                                  ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAtgttg                                                                                           >  2:238882/1‑141 (MQ=255)
                                                  ttatcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGAtgttg                                                                                           <  2:281474/141‑1 (MQ=255)
                                                     tcattTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGa                                                                                        <  1:311254/141‑1 (MQ=255)
                                                        ttttATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGt                                                                                     >  2:204618/1‑141 (MQ=255)
                                                          ttattattaCCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTaa                                                                                   >  1:632683/1‑141 (MQ=255)
                                                                    ccTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGAt                                                                         <  1:554749/141‑1 (MQ=255)
                                                                          ttAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGGAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAg                                                                   <  1:378782/141‑1 (MQ=255)
                                                                            aaaaGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCt                                                                 <  2:205885/141‑1 (MQ=255)
                                                                              aaGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTAt                                                               <  1:731238/141‑1 (MQ=255)
                                                                              aaGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTAt                                                               >  1:539070/1‑141 (MQ=255)
                                                                                   tCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGtt                                                          <  1:70919/141‑1 (MQ=255)
                                                                                      acaaTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtata                                                       <  1:75721/141‑1 (MQ=255)
                                                                                       caaTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGGGAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtatat                                                      >  1:274174/1‑141 (MQ=255)
                                                                                       caaTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtatat                                                      >  1:607673/1‑141 (MQ=255)
                                                                                       caaTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtatat                                                      >  1:374077/1‑141 (MQ=255)
                                                                                       caaTACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTtatat                                                      >  2:238123/1‑141 (MQ=255)
                                                                                            cAACTGAAGACTTCTTCATTTGATGCGCAAAAAAAGCGGCGAGTCAGTTGCGATCCAGCAGAGAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGctc                                                 <  1:271185/141‑1 (MQ=255)
                                                                                                           tcATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTaa                                  >  1:134559/1‑141 (MQ=255)
                                                                                                             aTTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATc                                <  1:259992/141‑1 (MQ=255)
                                                                                                                ttATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAATCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCt                             <  2:333574/141‑1 (MQ=255)
                                                                                                                ttATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCt                             <  2:25854/141‑1 (MQ=255)
                                                                                                                  aTGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTcc                           >  2:254736/1‑141 (MQ=255)
                                                                                                                   tGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTccc                          >  1:717273/1‑141 (MQ=255)
                                                                                                                      aTAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCt                       <  1:183442/141‑1 (MQ=255)
                                                                                                                               tCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGt              <  1:231115/141‑1 (MQ=255)
                                                                                                                                cGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTa             >  2:207123/1‑141 (MQ=255)
                                                                                                                                cGGCTAGCACGTGACTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTa             >  2:315681/1‑141 (MQ=255)
                                                                                                                                  gCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAat           <  1:524995/141‑1 (MQ=255)
                                                                                                                                   cTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAata          >  2:91897/1‑141 (MQ=255)
                                                                                                                                     tagcACGTGACTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAatat         >  2:89439/4‑141 (MQ=255)
                                                                                                                                           gTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAATATCATAGca  >  2:145128/1‑141 (MQ=255)
                                                                                                                                           gTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAATATCATAGca  >  2:145124/1‑141 (MQ=255)
                                                                                                                                           gTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAATATCATAGca  >  2:117527/1‑141 (MQ=255)
                                                                                                                                           gTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAATATCATAGca  >  1:782656/1‑141 (MQ=255)
                                                                                                                                              |                                                                                                                                          
TATTTGCTGTAATGTTTTATCTGTCATTTCGTCCACCACTTCTCATATCATTATCATTTTATTATTACCCTATATTAAAAGAATCAACAATACAACTGAAGACTTCTTCATTTTATGCATAAAAAAATCGGCTAGTCACGTGCTAGCCGACAAATAGAAAGGAAAGTAAGTAATAAATATTGAAGATGTTGTGATGTAACTTGAACGATTAAAAGCTATCTGTTATATAGCTCTACCCCTTTGTTTAATCGCTCCCCCTGTTACAAGTAATATCATAGCA  >  USA300TCH1516_ALE/2430550‑2430829

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