Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,364,760 A→G 6.2% S111S (AGT→AGC qcrA ← menaquinol:cytochrome c oxidoreductase (iron‑sulfur subunit)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,364,7600AG 6.2% 188.5 / 3.9 65S111S (AGT→AGCqcrAmenaquinol:cytochrome c oxidoreductase (iron‑sulfur subunit)
Reads supporting (aligned to +/- strand):  ref base A (25/35);  new base G (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GAAGCCGTCTTTTACTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAAT  >  BS168_reference_genome/2364628‑2364896
                                                                                                                                    |                                                                                                                                        
gAAGCCGTCTTTTACTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGGCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTATGATACCTATTc                                                                                                                                    >  6:1294/1‑139 (MQ=255)
  aGCCGTCTTTTACTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCa                                                                                                                                  <  5:109826/139‑1 (MQ=255)
  aGCCGTCTTTTACTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCa                                                                                                                                  <  3:330445/139‑1 (MQ=255)
     cGTCTTTTACATCTTGCTCATAATGATCAAGCGGCGCACGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGtt                                                                                                                               <  5:58195/139‑1 (MQ=255)
      gTCTTTTACTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGATTTTATTCGGGTTTTTAGGATCACTATTCCAGTTc                                                                                                                              <  5:350880/139‑1 (MQ=255)
            gaCTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTa                                                                                                                        <  7:226290/138‑1 (MQ=255)
             aCTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTa                                                                                                                        <  7:200131/138‑1 (MQ=255)
             aCTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTa                                                                                                                        <  5:75802/138‑1 (MQ=255)
             aCTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTa                                                                                                                        <  3:135008/138‑1 (MQ=255)
              cttcttGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTac                                                                                                                       <  3:151270/138‑1 (MQ=255)
                        aTAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGt                                                                                                                                >  8:128028/1‑119 (MQ=255)
                        aTAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGt                                                                                                                                <  7:128028/119‑1 (MQ=255)
                        aTAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACAATTACAGTTAACCGTACACCCTAAATGc                                                                                                            >  8:226290/1‑139 (MQ=255)
                               tCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCa                                                                                                                             >  3:177935/1‑115 (MQ=255)
                               tCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCa                                                                                                                             <  4:177935/115‑1 (MQ=255)
                                    agGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTGCGTAATGGCAGGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCAAGTTCACCGTACACACTAAATGCTTACAAATTGGc                                                                                                >  8:52815/1‑139 (MQ=255)
                                    agGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGc                                                                                                >  8:101836/1‑139 (MQ=255)
                                    agGCGAAAGAGGAGGCGTGCCAGGGACATTGGTGCAGTCCTTTTAATAGATTCCGTAATTGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTGACCGTACACCCTAAATGCATACAAATTGGc                                                                                                >  8:33815/1‑139 (MQ=255)
                                         cAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAg                                                                                           <  8:287538/139‑1 (MQ=255)
                                                   cgTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTAAAAATTGGCGAAAGTGCTACAAtt                                                                                 <  4:18456/139‑1 (MQ=255)
                                                               aTTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTa                                                                                      >  1:260055/1‑122 (MQ=255)
                                                               aTTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTa                                                                                      <  2:260055/122‑1 (MQ=255)
                                                               aTTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCGTACAAATTGGCGAAAGTGCTACACTTTCATACCCGgtt                                                                     >  2:302595/1‑139 (MQ=255)
                                                                   gTGCCGTCCTTTTCCTAGAGTCCGGAATGGCCTGGACAAAAGAATTTATTCTGGTTTTTAGGCTCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTg                                                                  <  3:240722/138‑1 (MQ=255)
                                                                   gTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTg                                                                  <  7:199976/138‑1 (MQ=255)
                                                                   gTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTg                                                                  <  8:109719/138‑1 (MQ=255)
                                                                         tCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTa                                                                                                                        >  6:24601/1‑78 (MQ=255)
                                                                         tCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTa                                                                                                                        <  5:24601/78‑1 (MQ=255)
                                                                         tCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATTCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGAcc                                                            >  8:286746/1‑138 (MQ=255)
                                                                         tCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGAccc                                                           >  8:24095/1‑139 (MQ=255)
                                                                           cTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAgg                                                        <  1:279107/139‑1 (MQ=255)
                                                                                 aTAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGAt                                                   <  3:309845/139‑1 (MQ=255)
                                                                                   agagTCCATAATGGCATGGACAAAAGAATTTATTCGGATTTTTAGGATCGCTATTCCAGTTCACCGTCCACCCTAAATGCTTACAAATTGGCGAAAGCGCTACAATTTCATCCCCGCTTTTCTAGACCCAGGCTGAGCt                                                 >  2:346258/1‑139 (MQ=255)
                                                                                     agTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCtt                                                >  3:33435/1‑138 (MQ=255)
                                                                                      gTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTGACGGTACAGCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGa                                              <  7:33815/139‑1 (MQ=255)
                                                                                      gTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAACTGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGa                                              <  2:303810/139‑1 (MQ=255)
                                                                                      gTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGa                                              <  1:279786/139‑1 (MQ=255)
                                                                                      gTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGa                                              <  7:101836/139‑1 (MQ=255)
                                                                                      gTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATATCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGa                                              <  2:299353/139‑1 (MQ=255)
                                                                                         cGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTa                                                                                                                        <  8:111407/62‑1 (MQ=255)
                                                                                         cGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTa                                                                                                                        >  7:111407/1‑62 (MQ=255)
                                                                                           tAATGTCATGGACAAAAGTATTTATTCGGGTTTGTAGGATCACTATTCCAGTTCACCGTACACACTATATGCTTACAAATTGTCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGAGtctt                                         <  4:115993/139‑1 (MQ=255)
                                                                                           tAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGActct                                          >  6:51943/1‑138 (MQ=255)
                                                                                           tAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACtctt                                         <  1:281905/139‑1 (MQ=255)
                                                                                                 cATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCt                                                      >  4:165547/1‑120 (MQ=255)
                                                                                                 cATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCt                                                      <  3:165547/120‑1 (MQ=255)
                                                                                                 cATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACt                                   >  3:75396/1‑139 (MQ=255)
                                                                                                 cATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACt                                   >  5:46953/1‑139 (MQ=255)
                                                                                                 cATGGACAAAAGAATTTATTCGGATTTTTAGGATCGCTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGCGCTACAATTTCATCCCCGCTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGAc                                    >  7:283585/1‑138 (MQ=255)
                                                                                                   tGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATCCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTTACTCTTCTGACTc                                  >  6:39587/1‑138 (MQ=255)
                                                                                                     gACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATa                               >  7:38436/1‑139 (MQ=255)
                                                                                                       cAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATAcc                             <  7:84105/139‑1 (MQ=255)
                                                                                                       cAAAAGAATTTATTCGGATTTTTAGGATCGCTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGCGCTACAATTTCATCCCCGCTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATAcc                             <  8:283585/139‑1 (MQ=255)
                                                                                                          aaGAATTTATTCGGGTTTTTAGGCTCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGATAGTGCTACAATTTCATCCCCGTTGTTGAAGACCCAGGCTGATCTTGACTCTTCTGGCTCGTACCAgg                          >  8:332792/1‑139 (MQ=255)
                                                                                                          aaGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAACTGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAgg                          >  4:83432/1‑139 (MQ=255)
                                                                                                          aaGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAgg                          >  5:149618/1‑139 (MQ=255)
                                                                                                            gAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCa                        <  5:1294/139‑1 (MQ=255)
                                                                                                               tttATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATc                      >  4:326909/1‑138 (MQ=255)
                                                                                                                  aTTCGGGTTTGTAGGATCACTATTACAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGAc                   <  7:118160/138‑1 (MQ=255)
                                                                                                                       ggTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGAtt             >  2:177706/1‑139 (MQ=255)
                                                                                                                              aGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAAtttt       <  6:122345/138‑1 (MQ=255)
                                                                                                                               ggATCGCTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGCGCTACAATTTCATCCCCGCTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGGTTAAttttg      <  8:203352/138‑2 (MQ=255)
                                                                                                                                  tCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGa                                              <  3:187742/95‑1 (MQ=255)
                                                                                                                                  tCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGa                                              >  4:187742/1‑95 (MQ=255)
                                                                                                                                  tCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAAt  <  4:75396/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                        
GAAGCCGTCTTTTACTTCTTGCTCATAATGATCAAGAGGCGCAAGCGGCGGCGTGCCAGGGACATTGGTGCCGTCCTTTTCATAGAGTCCGTAATGGCATGGACAAAAGAATTTATTCGGGTTTTTAGGATCACTATTCCAGTTCACCGTACACCCTAAATGCTTACAAATTGGCGAAAGTGCTACAATTTCATCCCCGTTTTTGAAGACCCAGGCTGATCTTGACTCTTCTGACTCATACCAGGCATCGACTTGATTAATTTTAAAAT  >  BS168_reference_genome/2364628‑2364896

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: