Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,421,950 C→T 6.5% A75T (GCG→ACG)  rluB ← 23S rRNA pseudouridine 2633 (=2605 standard) pseudouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,421,9500CT 6.5% 191.1 / 3.4 61A75T (GCG→ACG) rluB23S rRNA pseudouridine 2633 (=2605 standard) pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base C (22/35);  new base T (2/2);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GCTTATTGGCGAATTCGCCGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACACCTAGCTCTGTTACTACTT  >  BS168_reference_genome/2421812‑2422075
                                                                                                                                          |                                                                                                                             
gCTTATTGGCGAATTCGCCGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTgcgc                                                                                                                               <  2:100498/139‑1 (MQ=255)
 cTTATTGGCGAATTCGCCGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGCCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCt                                                                                                                              >  1:6787/1‑139 (MQ=255)
  ttATTGGCGAATTCGCCGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTg                                                                                                                             <  8:94702/139‑1 (MQ=255)
           aaTTCGCCGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGCCGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCCCTTTGCGTCCTTTATCATCTTTCGCTGCGGAAATaa                                                                                                                    >  7:277223/1‑139 (MQ=255)
           aaTTCGCCGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGAGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACATTGCGTCCTTTATCATCTTGCGCTGCGGaaaaaa                                                                                                                    >  5:29813/1‑139 (MQ=255)
               cGCCGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGcc                                                                                                                >  6:103539/1‑139 (MQ=255)
                 ccGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCctc                                                                                                              <  8:80334/139‑1 (MQ=255)
                 ccGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCctc                                                                                                              <  8:172266/139‑1 (MQ=255)
                  cGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCctct                                                                                                             <  3:217168/139‑1 (MQ=255)
                      aTTTGTTAAACCCAAAAGACCGCTCGTATAATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGt                                                                                                         <  5:84695/139‑1 (MQ=255)
                      aTTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGt                                                                                                         <  1:244550/139‑1 (MQ=255)
                      aTTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGt                                                                                                         <  8:154909/139‑1 (MQ=255)
                      aTTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGt                                                                                                         <  6:120576/139‑1 (MQ=255)
                       ttGGTTAAAAGCAAAAGCCCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGtt                                                                                                        >  2:100409/1‑139 (MQ=255)
                          gTTAAAAGGAAAAGACCGCTCGTATCCTAGTCCAGCCGCCCAATAGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTACTTTATCATCTTGCGCTGCGGAAATAACGCCGCATGGTTta                                                                                                      >  6:257170/1‑138 (MQ=255)
                            aaaaaCCAAAAGCAGGCGCGTCTCTTAGTCCTCACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGGGTCCTTTATCATCTTGCGCTGCGGAATTAACGCGTCTTGGTGtata                                                                                                    <  6:45180/137‑1 (MQ=255)
                               aaGCAAACGACCGCGCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAag                                                                                                <  8:277223/139‑1 (MQ=255)
                                aGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAaga                                                                                               <  2:49490/139‑1 (MQ=255)
                                    aaaGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAATTGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAgaga                                                                                             >  1:38951/1‑137 (MQ=255)
                                    aaaGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAgaga                                                                                             >  1:38945/1‑137 (MQ=255)
                                           gCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCg                                                                                                                           <  1:163441/100‑1 (MQ=255)
                                           gCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCg                                                                                                                           >  2:163441/1‑100 (MQ=255)
                                            cTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGTTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCg                                                                                   >  5:265937/1‑139 (MQ=255)
                                            cTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCg                                                                                   >  1:234733/1‑139 (MQ=255)
                                            cTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCg                                                                                   >  5:203884/1‑139 (MQ=255)
                                            cTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACAACTTTGCGTCCTTTATCATCTTGCGTTGAGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAAcag                                                                                   >  5:166351/1‑137 (MQ=255)
                                                    aTAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATAATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTcg                                                                           >  3:252345/1‑139 (MQ=255)
                                                        tCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGtt                                                                       <  3:65052/139‑1 (MQ=255)
                                                         ccAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTc                                                                      <  5:103539/139‑1 (MQ=255)
                                                          cAGACGCCCAATCGTATAAATGCGGTGCGTAATATTTATAAAAATGTCCGTCACATCTTTGCTTCCTTTATCAACTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCTCGTTc                                                                      <  2:226583/138‑1 (MQ=255)
                                                          cAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTc                                                                      <  6:116824/138‑1 (MQ=255)
                                                          cAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTc                                                                      <  5:279670/138‑1 (MQ=255)
                                                          cAGACGCCCAATCGGATAAAATCGCTGCGGAATATATTTAAAAAAGTCCGTCACCACTTTTCGTCCTATATCATCTTGCGCTGCGGAAATAAAGCCTCTAGGTTTATATAAGAGAAAGTAAACCTGCTCTTCACGtac                                                                      >  1:226583/1‑136 (MQ=255)
                                                             aCGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGt                                                                                                         <  7:61749/100‑1 (MQ=255)
                                                             aCGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGt                                                                                                         >  8:61749/1‑100 (MQ=255)
                                                             aCGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAAcc                                                                  >  1:48861/1‑139 (MQ=255)
                                                             aCGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAAcc                                                                  >  6:207527/1‑139 (MQ=255)
                                                                cccAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCtt                                                                >  2:159938/1‑138 (MQ=255)
                                                                   aaTCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCTCTGAGTAAATAACGCGTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAAc                                                            <  3:131755/139‑1 (MQ=255)
                                                                   aaTCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAAc                                                            <  2:174072/139‑1 (MQ=255)
                                                                       ggATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGt                                                        <  5:207527/139‑1 (MQ=255)
                                                                           aaaTGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTAc                                                    >  6:253409/1‑139 (MQ=255)
                                                                            aaTGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTAcc                                                   <  6:203884/139‑1 (MQ=255)
                                                                             aTGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCt                                                  >  6:8072/1‑139 (MQ=255)
                                                                                 gcTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCaa                                               >  2:265364/1‑138 (MQ=255)
                                                                                 gcTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAAt                                              <  7:226434/139‑1 (MQ=255)
                                                                                   tGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCGTTAACCCGTTTCCCTCAATCt                                            >  2:269619/1‑139 (MQ=255)
                                                                                      ggAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGAt                                         <  1:100409/139‑1 (MQ=255)
                                                                                                 aaaaacGTCCGTCACCACTTTGCGTCGTTTATCATGTTGCTTTGCGTAAATAACGCCGATTTGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGt                               <  6:166351/138‑1 (MQ=255)
                                                                                                 aaaaaaGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGTTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGt                               <  6:265937/138‑1 (MQ=255)
                                                                                                 aaaaaaGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGt                               <  5:8072/138‑1 (MQ=255)
                                                                                                 aaaaaaGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGt                               <  5:253409/138‑1 (MQ=255)
                                                                                                 aaaaaaGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGt                               <  2:6787/138‑1 (MQ=255)
                                                                                                       gTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTa                        <  4:252345/139‑1 (MQ=255)
                                                                                                       gTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTa                        <  2:3588/139‑1 (MQ=255)
                                                                                                            tCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACAcc                   <  1:265364/139‑1 (MQ=255)
                                                                                                               ccacTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACACCTAg                <  2:234733/139‑1 (MQ=255)
                                                                                                                cacTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACACCTAGc               <  7:206251/139‑1 (MQ=255)
                                                                                                                cacTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACACCTAGc               <  5:257170/139‑1 (MQ=255)
                                                                                                                             ttATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACACCTAGCTCTGTTACTACtt  >  7:23576/1‑139 (MQ=255)
                                                                                                                              tATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACACCTAGCTCTGTTACTACtt  <  2:38951/138‑1 (MQ=255)
                                                                                                                              tATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGGTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACACCTAGCTCTGTTACTACtt  <  2:38945/138‑1 (MQ=255)
                                                                                                                                          |                                                                                                                             
GCTTATTGGCGAATTCGCCGTCATTGGTTAAAAGCAAAAGACCGCTCGTATCATAGTCCAGACGCCCAATCGGATAAATGCGCTGCGGAATATTTTTAAAAAAGTCCGTCACCACTTTGCGTCCTTTATCATCTTGCGCTGCGGAAATAACGCCTCTTGGTTTATATAAGAGAAAGTAAACCGGCTCTTCGCGTTCAACCTTTAACCCGTTTACCTCAATCTGATCAGAACCGGTGACCTTAACACCTAGCTCTGTTACTACTT  >  BS168_reference_genome/2421812‑2422075

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