Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,674,958 A→T 13.3% intergenic (+11/+14) 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,9580AT13.3% 158.4 / 13.1 52intergenic (+11/+14)tqsA/pntBpheromone AI‑2 transporter/pyridine nucleotide transhydrogenase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base T (3/4);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACG  >  NC_000913/1674812‑1675104
                                                                                                                                                  |                                                                                                                                                  
tttttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCaga                                                                                                                                                  <  2:1034850/149‑1 (MQ=255)
   ttGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGc                                                                                                                                               <  3:1669834/149‑1 (MQ=255)
       ggATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCGTc                                                                                                                                           <  3:1529639/149‑1 (MQ=255)
             ttgttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAgg                                                                                                                                                <  1:486086/138‑1 (MQ=255)
             ttgttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAgg                                                                                                                                                >  2:486086/1‑138 (MQ=255)
              tgttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCg                                                                                                                                             <  3:1692378/140‑1 (MQ=255)
              tgttAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCTGCTGAGGCCg                                                                                                                                             >  4:1692378/1‑140 (MQ=255)
                    gACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGc                                                                                                                              >  4:1430479/1‑149 (MQ=255)
                       cGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAgg                                                                                                                                                >  1:107470/1‑128 (MQ=255)
                       cGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAgg                                                                                                                                                <  2:107470/128‑1 (MQ=255)
                            ggTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGtt                                                                                                                                     >  4:1697657/1‑134 (MQ=255)
                            ggTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGtt                                                                                                                                     <  3:1697657/134‑1 (MQ=255)
                                      tttCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCAc                                                                                                            <  3:1430479/149‑1 (MQ=255)
                                                              tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                                                                                 >  4:587183/1‑88 (MQ=255)
                                                              tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCagag                                                                                                                                                 <  3:587183/88‑1 (MQ=255)
                                                              tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAggg                                                                                                                                       <  4:813709/98‑1 (MQ=255)
                                                              tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAggg                                                                                                                                       >  3:813709/1‑98 (MQ=255)
                                                              tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAggg                                                                                                                                       <  4:361108/98‑1 (MQ=255)
                                                              tCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAggg                                                                                                                                       >  3:361108/1‑98 (MQ=255)
                                                                aaaaTTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCa                                                                                  >  4:936406/1‑149 (MQ=255)
                                                                       cgcTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGgt                                                                           <  1:598850/149‑1 (MQ=255)
                                                                          ttGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAgg                                                                                                                                        <  3:534843/85‑1 (MQ=255)
                                                                          ttGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAgg                                                                                                                                        >  4:534843/1‑85 (MQ=255)
                                                                           tGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGttt                                                                       >  1:1717083/1‑149 (MQ=255)
                                                                             aacaaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTc                                                                     >  3:771566/1‑149 (MQ=255)
                                                                              acaaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTcc                                                                    >  2:777731/1‑149 (MQ=255)
                                                                               caaacaGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCt                                                                   <  1:777731/149‑1 (MQ=255)
                                                                                             tCAAAGCATCGCCGTTCTGTTAAGCGCTCTCAATAAAGAGTGACGGCCTCAGCAGAGGCTGTCAGGCTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTc                                                                     <  4:1227495/133‑1 (MQ=255)
                                                                                             tCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTc                                                                     >  3:1227495/1‑133 (MQ=255)
                                                                                              cAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcagcag                                                                                                                                                   >  3:1234372/1‑54 (MQ=255)
                                                                                              cAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTcagcag                                                                                                                                                   <  4:1234372/54‑1 (MQ=255)
                                                                                                    aTCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTa                                                                                                                                    <  3:1295102/63‑1 (MQ=255)
                                                                                                    aTCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTa                                                                                                                                    >  4:1295102/1‑63 (MQ=255)
                                                                                                       gCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCacac                                           <  3:936406/149‑1 (MQ=255)
                                                                                                       gCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCacac                                           <  4:771566/149‑1 (MQ=255)
                                                                                                          gTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGc                                               >  3:363181/1‑142 (MQ=255)
                                                                                                          gTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGc                                               <  4:363181/142‑1 (MQ=255)
                                                                                                                  aaGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAtgt                                                                               <  2:569830/102‑1 (MQ=255)
                                                                                                                  aaGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCAtgt                                                                               >  1:569830/1‑102 (MQ=255)
                                                                                                                          tcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGtt                           >  4:1181887/1‑146 (MQ=255)
                                                                                                                          tcAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGtt                           <  3:1181887/146‑1 (MQ=255)
                                                                                                                           cAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTg                                                <  4:1337382/124‑1 (MQ=255)
                                                                                                                           cAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTg                                                >  3:1337382/1‑124 (MQ=255)
                                                                                                                                         ggCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGc                                  >  3:1389347/1‑124 (MQ=255)
                                                                                                                                         ggCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGc                                  <  4:1389347/124‑1 (MQ=255)
                                                                                                                                         ggCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaaa            <  2:1550166/146‑1 (MQ=255)
                                                                                                                                         ggCCTCTGCTGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTaaa            >  1:1550166/1‑146 (MQ=255)
                                                                                                                                         ggCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTg                                                >  4:1732470/1‑110 (MQ=255)
                                                                                                                                         ggCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTg                                                <  3:1732470/110‑1 (MQ=255)
                                                                                                                                          gCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGaa                    <  3:327064/137‑1 (MQ=255)
                                                                                                                                          gCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGaa                    >  4:327064/1‑137 (MQ=255)
                                                                                                                                           ccTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATc                       >  2:287127/1‑133 (MQ=255)
                                                                                                                                           ccTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATc                       <  1:287127/133‑1 (MQ=255)
                                                                                                                                           ccTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACaa       >  1:464139/1‑149 (MQ=255)
                                                                                                                                            cTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtt                                                   >  3:961489/1‑104 (MQ=255)
                                                                                                                                            cTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGtt                                                   <  4:961489/104‑1 (MQ=255)
                                                                                                                                             tCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATc     >  1:20405/1‑149 (MQ=255)
                                                                                                                                               agcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTgg                                                                             <  2:490246/75‑1 (MQ=255)
                                                                                                                                               agcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTgg                                                                             >  1:490246/1‑75 (MQ=255)
                                                                                                                                               agcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCAc   >  4:112410/1‑149 (MQ=255)
                                                                                                                                               agcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCAc   >  2:1095825/1‑149 (MQ=255)
                                                                                                                                                gcagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACg  >  1:659432/1‑149 (MQ=255)
                                                                                                                                                 cagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTc                          >  3:714640/1‑124 (MQ=255)
                                                                                                                                                 cagAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTc                          <  4:714640/124‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTTTGGGGATGGTTGTTAGGACCGGTGGGTATGCTGCTTTCCGTGCCGTTGACAATTATTGTCAAAATTGCGCTTGAACAAACAGCGGGAGGTCAAAGCATCGCCGTTCTGTTAAGCGATCTCAATAAAGAGTGACGGCCTCAGCAGAGGCCGTCAGGGTTACAGAGCTTTCAGGATTGCATCCACGCTGGCTTTGGCGTCACCAAACAGCATGTGGGTGTTTTCCTTGAAGAACAGCGGGTTTTGCACACCAGCATAGCCAGTGTTCATCGAACGTTTAAAGACAATCACG  >  NC_000913/1674812‑1675104

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: