Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,955 A→T 11.5% intergenic (+8/+17) tqsA → / ← pntB pheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,674,9550AT11.5% 106.1 / 4.7 49intergenic (+8/+17)tqsA/pntBpheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (27/15);  new base T (2/4);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCAC  >  NC_000913/1674807‑1675103
                                                                                                                                                    |                                                                                                                                                    
ttGATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTca                                                                                                                                                      >  1:175635/1‑149 (MQ=255)
 tGATTTTTTGGGGATTGTTGTTAGGTCCTGTGGTTAGGCGGCTTTACGTGCCGTTGACAATTATTGTCAATATTGCGCTTGAACCAACAGCGGGAGGGCAAAGCATCGCCGTGCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcag                                                                                                                                                     <  1:606357/149‑1 (MQ=255)
    ttttttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCATcagcag                                                                                                                                                  >  1:1021546/1‑149 (MQ=255)
      ttttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGCGTGACGGCCTcagccgcg                                                                                                                                                >  1:350893/1‑145 (MQ=255)
           gggATGGTTGTTAGGACCGGGGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGt                                                                                                                                           <  1:37997/149‑1 (MQ=255)
            ggATGGTTGTTAGGACCGGGGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTGGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAAAAAAGAGTGACGGCCTCTGCTGAGGCCGTc                                                                                                                                          <  1:746779/149‑1 (MQ=255)
              aTGGTTGTGAGGACCGGGGGGTATGCTGCTTTCCGTGCCGTAGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGGCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGTCac                                                                                                                                        <  1:259892/149‑2 (MQ=255)
                   tgttATGACCGGGGGGTATGCTGCTTTCCGGGCCGTTGACAATTATTGGAAAAATTGCGCTTGAACAAACAGCGGGATGTCAAAGCATGGCCGTTCTGTTAATCGATCTCAATAAAGCGTGACGGCCTCAGCAGATGCCGTCATGGTTa                                                                                                                                   <  2:962849/149‑1 (MQ=255)
                       aGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATAGCGCTGGAACACACAGCGGGAGGGCAAAGAATCGCCGTTCTGTTAAGCGATCTCATTAAAGAGTGACGGCCTCAGACGAGGCCCTCAGGGTTACaga                                                                                                                               >  2:582923/1‑149 (MQ=255)
                         gACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGc                                                                                                                             >  1:498126/1‑149 (MQ=255)
                           ccGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCtt                                                                                                                           >  2:497640/1‑149 (MQ=255)
                                       ctgctTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCGGAGGCCGTCAGGGTTACAGCGCTTTCAGGATTGCAt                                                                                                               >  1:822363/1‑149 (MQ=255)
                                           tttCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCAc                                                                                                           >  1:363296/1‑149 (MQ=255)
                                             tCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGc                                                                                                         <  2:175635/149‑1 (MQ=255)
                                                gTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGGTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCGTTCAGGATTGCATCCACGCTgg                                                                                                      >  2:29859/1‑149 (MQ=255)
                                                 tGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAAGATTGCATCCACGCTGGc                                                                                                     >  2:907009/1‑149 (MQ=255)
                                                                   tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAg                                                                                   >  1:67282/1‑149 (MQ=255)
                                                                      aaaTTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAt                                                                                <  1:497640/149‑1 (MQ=255)
                                                                       aaTTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATAGCCGTTCAGATAAGCGATCTCAATAAAGAATGAAGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGAATCCATCCAAGCAGGCTTTGGCGTCACCAAAcatcatc                                                                               >  1:732163/1‑148 (MQ=255)
                                                                             gcTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACTTCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGgtg                                                                         >  1:473556/1‑149 (MQ=255)
                                                                                      aacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTg                                                                >  1:800558/1‑149 (MQ=255)
                                                                                       acaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGGTTTCCTTga                                                               >  1:543350/1‑149 (MQ=255)
                                                                                        caGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTAACATATCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTAACCAAAAAGCATGTGCGTTTTTTCCTTgaa                                                              >  1:605915/1‑149 (MQ=255)
                                                                                         aGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTgaag                                                             >  2:730455/1‑149 (MQ=255)
                                                                                           cGGGAGGTCAAAGCATCGCCGTTCTGTTAATCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCATGGTTACAGAGCTTTCAGGATTTCCTCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGTTGTGTTCCTTgaagaa                                                           >  2:480621/1‑149 (MQ=255)
                                                                                                      aGCATGGCCGGTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtttt                                                <  1:131147/149‑1 (MQ=255)
                                                                                                      aGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCAGCAAACAGCATGTGGGTGTTTTCCTTGAAGATCAGCGGGtttt                                                >  2:535392/1‑149 (MQ=255)
                                                                                                        cATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGc                                              <  1:827616/149‑1 (MQ=255)
                                                                                                        cATCGCCGTGCTGTTAAGCGATAGCAATAATGCGTGACTGCCTCAGCAGAGGCCGTCAGGGTTCCAGCGCTTTCAGGATTGCATCCACGCGGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGc                                              <  1:738042/149‑1 (MQ=255)
                                                                                                         agcGCCGGTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCCGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGca                                             <  2:350893/147‑1 (MQ=255)
                                                                                                         aTCGCCGTTCTGTTAAGCGATCTCAATAAAGAGGGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGca                                             <  1:398276/149‑1 (MQ=255)
                                                                                                          tCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGcac                                            >  1:6050/1‑149 (MQ=255)
                                                                                                                    gTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGGTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAACAACAGCGGGTTTTGCTCACCCGCagag                                  >  1:51497/1‑149 (MQ=255)
                                                                                                                     ttAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGGCGTCAGGGGTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGc                                 <  2:435363/149‑1 (MQ=255)
                                                                                                                       aaGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCCCCAACCAGCATGTGGGTGTTTTCCTTCAAGAACAGCGGGTTTTGCACACCCGCATAGCCa                               >  1:143138/1‑149 (MQ=255)
                                                                                                                        aGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAg                              <  1:730455/149‑1 (MQ=255)
                                                                                                                        aGCGAGCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAg                              <  2:67282/149‑1 (MQ=255)
                                                                                                                          cGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCTTTGAAGCACAGCGGGTTTTGCACACCAGCATAGCCAgtg                            >  2:248799/1‑149 (MQ=255)
                                                                                                                            aTCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGtt                          <  2:363296/149‑1 (MQ=255)
                                                                                                                               tcAATAAAGCGGGTCGGCCTCTTCAGAGGCCGTCAGGGTTTCAGAGCTTTCATGATTGCATCCACGCTGGCTTTGGCGTCACCAATCAGCATTTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCAt                       <  2:732163/149‑1 (MQ=255)
                                                                                                                                cAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATc                      >  1:4118/1‑149 (MQ=255)
                                                                                                                                 aaTAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGTCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTAATCg                     >  2:356707/1‑149 (MQ=255)
                                                                                                                                   tAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGaa                   >  1:340817/1‑149 (MQ=255)
                                                                                                                                       gagTGACGGTCTCAGCAGAGTCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCAGTCTTGCTTTGGCGTCACCAAACAGCATGTTTGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGtt               <  1:258017/149‑1 (MQ=255)
                                                                                                                                        cctgccGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTGTGGCGTCGCCCACCAGCATGTGGGTGTTTTCCTTGTAGATAAGCGGGTTTTGCTCACCAGCCTCGCCAGTGTTCATCGAACGttt              >  1:691088/6‑149 (MQ=255)
                                                                                                                                         ctGACGGCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTa             >  1:170655/2‑149 (MQ=255)
                                                                                                                                           gACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaaa           <  1:385673/149‑1 (MQ=255)
                                                                                                                                            aCGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaa            >  2:578351/1‑147 (MQ=255)
                                                                                                                                            aCGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaa            <  1:578351/147‑1 (MQ=255)
                                                                                                                                            aCGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTGGAAGAACAGCGGGGTTTGCACACCCGCATAGCCAGTGTTGATCGCACGTTTAAAg          >  2:180616/1‑149 (MQ=255)
                                                                                                                                             cGGCCGCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCGGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTGTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGa         <  1:726073/149‑1 (MQ=255)
                                                                                                                                                ccTCAGCAGAGGCCGTCAGCGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGGTCCTCGCACGTTTACAGACaa      >  2:613895/1‑149 (MQ=255)
                                                                                                                                                 cTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTACATCCACGCTGGCTTTGGCGTCACCAAACAGCTTGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAAt     >  1:309658/1‑149 (MQ=255)
                                                                                                                                                   cagcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGTTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCa   <  2:1021546/149‑1 (MQ=255)
                                                                                                                                                   cagcagAGGCCGTCAGGGCTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCa   >  1:1007459/1‑149 (MQ=255)
                                                                                                                                                    agcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCAc  >  2:384798/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGATTTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCAC  >  NC_000913/1674807‑1675103

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: