Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 185,723 T→A 6.9% H139L (CAT→CTT)  PP_0162 ← ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947185,7230TA 6.9% 161.0 / 3.2 58H139L (CAT→CTT) PP_0162ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base A (2/2);  total (27/31)
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.34e-01

TGACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGG  >  NC_002947/185579‑185866
                                                                                                                                                |                                                                                                                                               
tgACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTAAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGAACGGAGATCCCCATCTCTTCGGCGATGGCCTGATGCt                                                                                                                                              >  6:48093/1‑148 (MQ=255)
  aCCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCa                                                                                                                                           >  8:78728/1‑149 (MQ=255)
   ccGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCATCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTAAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCa                                                                                                                                           <  5:48093/148‑1 (MQ=255)
   ccGTGTTTGTATCAGGGCGGGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAg                                                                                                                                          <  4:296186/149‑1 (MQ=255)
   acttttttGTATCAGGGCGTGCGACGGTAATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCa                                                                                                                                           <  8:91005/145‑1 (MQ=255)
         ttGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCt                                                                                                                                     <  7:274354/148‑1 (MQ=255)
                aGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGAGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGCTCCCCATCTGTTCGGCTATGGCCTGATGCTTCAGCCCCTCGCATTGa                                                                                                                             <  5:78449/149‑1 (MQ=255)
                   gCGTGCGACGGTCATGGATTCAGCAGTGTTCTTCATAGCGCAGCACGTAGCAGTTATAGAGCGCATTGGCGATGTATCGCTCAAGGGTGCGCAGGGATATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGc                                                                                                                           <  7:171213/148‑1 (MQ=255)
                    cGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCca                                                                                                                         <  2:183490/149‑1 (MQ=255)
                                 tGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGcgcg                                                                                                            <  5:269803/149‑1 (MQ=255)
                                   gATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGAc                                                                                                          >  6:78802/1‑149 (MQ=255)
                                        agcagTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTAcc                                                                                                     >  6:247877/1‑149 (MQ=255)
                                          cagTCGTCTTCATAGCGCAGCACCTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCtt                                                                                                   >  2:63344/1‑149 (MQ=255)
                                           agTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTg                                                                                                  >  6:46455/1‑149 (MQ=255)
                                           agTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTg                                                                                                  <  1:133492/149‑1 (MQ=255)
                                               gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                              >  4:332920/1‑149 (MQ=255)
                                               gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                              <  4:60362/149‑1 (MQ=255)
                                               gtcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGtt                                                                                              >  1:322943/1‑149 (MQ=255)
                                                tcTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGttt                                                                                             <  2:322943/149‑1 (MQ=255)
                                                  ttCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGAACCTTGGGttt                                                                                             >  7:54499/1‑147 (MQ=255)
                                                     aTAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAgg                                                                                         <  3:332920/148‑1 (MQ=255)
                                                     aTAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAgg                                                                                         <  5:46455/148‑1 (MQ=255)
                                                          gcagcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCACGAATGCGCGACGTACCTAGGGATTGAGGCCGTc                                                                                    >  3:297832/1‑148 (MQ=255)
                                                             gcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                                <  5:78802/149‑1 (MQ=255)
                                                             gcaCGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGc                                                                                <  1:219632/149‑1 (MQ=255)
                                                                    gCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCGCATCTGTTCGGCGCTGGCCTTATGCTTCAGCCCCTCGCATTGAGCCAGCATGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGt                                                                          <  1:158109/148‑1 (MQ=255)
                                                                          aTAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATAc                                                                   >  7:134791/1‑149 (MQ=255)
                                                                           tAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACg                                                                  >  2:257864/1‑149 (MQ=255)
                                                                                    ttGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATCCTTCAGCCCCACGCATTTAGCCAGCAGGAATCCCCGACGTACCTTGGGTTTGAGCCCGACGAGCATGCCTCCAATACGCTCAAGCaa                                                         >  3:277777/1‑149 (MQ=255)
                                                                                       gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTgg                                                                                                 <  2:257045/106‑1 (MQ=255)
                                                                                       gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTc                                                      >  3:178895/1‑149 (MQ=255)
                                                                                       gCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGCCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTgg                                                                                                 >  1:257045/1‑106 (MQ=255)
                                                                                         gATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                    <  4:178895/149‑1 (MQ=255)
                                                                                            gTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                          >  8:176606/1‑108 (MQ=255)
                                                                                            gTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAg                                                                                          <  7:176606/108‑1 (MQ=255)
                                                                                               gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                    >  7:23314/1‑143 (MQ=255)
                                                                                               gcgcTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCca                                                    <  8:23314/143‑1 (MQ=255)
                                                                                                             gcACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCTCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGtctc                                 >  8:139020/1‑148 (MQ=255)
                                                                                                             gcACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGtctc                                 >  7:209207/1‑148 (MQ=255)
                                                                                                              cACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGa                               >  3:152935/1‑149 (MQ=255)
                                                                                                              cACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGa                               <  4:297832/149‑1 (MQ=255)
                                                                                                                       ccccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCCCGACGAACCAGGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGa                      >  8:184928/1‑149 (MQ=255)
                                                                                                                        cccATCTGTTCGGCGATGGCGTGATGCTTCTGCCCCTCGCATTGAGCCAGCAGGAATTCGCGACGTACCTTGTGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTGCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                     <  4:277777/149‑1 (MQ=255)
                                                                                                                         ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAcc                    <  1:63344/149‑1 (MQ=255)
                                                                                                                         ccATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGAc                     >  1:284068/1‑148 (MQ=255)
                                                                                                                          cATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                          >  2:219328/1‑142 (MQ=255)
                                                                                                                          cATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAg                          <  1:219328/142‑1 (MQ=255)
                                                                                                                           aTCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAg                  <  8:209207/149‑1 (MQ=255)
                                                                                                                              ctcaatgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTg                <  7:178506/142‑1 (MQ=255)
                                                                                                                              ctcaatgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTg                >  8:178506/7‑148 (MQ=255)
                                                                                                                                ttCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCATGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGt                 <  1:257864/145‑1 (MQ=255)
                                                                                                                                  atgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCGGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc           >  6:249693/3‑149 (MQ=255)
                                                                                                                                  atgCGAGGGGCTGAAGCATCAGGCCATCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCGGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc           <  5:249693/147‑1 (MQ=255)
                                                                                                                                      gATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                                      <  6:209370/118‑1 (MQ=255)
                                                                                                                                      gATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCgg                                      >  5:209370/1‑118 (MQ=255)
                                                                                                                                           ccTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCTAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATgg  <  7:78728/149‑1 (MQ=255)
                                                                                                                                               aTGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc           <  4:329329/136‑1 (MQ=255)
                                                                                                                                               aTGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTcc           >  3:329329/1‑136 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGACCGTTTTTGTATCAGGGCGTGCGACGGTCATGGATTCAGCAGTCGTCTTCATAGCGCAGCACGTAGCAGTGATAGAGCGCATTGGCGATGTAGCGCTCAACGGTGCGCACGGAGATCCCCATCTGTTCGGCGATGGCCTGATGCTTCAGCCCCTCGCATTGAGCCAGCAGGAATGCGCGACGTACCTTGGGTTTGAGGCCGTCGAGCATGCGTGCAATACGCTCAAGCAATTCCAGCAGCAGCTCGCGGGTCTCGACTGAGGGGACCAGTGCTTCCGGCACATGG  >  NC_002947/185579‑185866

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