Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 282,578 A→T 8.3% N10I (AAT→ATT)  betT‑II → choline/carnitine/betaine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947282,5780AT 8.3% 175.2 / 5.4 60N10I (AAT→ATT) betT‑IIcholine/carnitine/betaine transporter
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base T (2/3);  total (29/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 = 8.97e-01

TAACTTGAACGCTTGTTCAATGAAAAACGCTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACGGTCGGCTGGTAC  >  NC_002947/282433‑282726
                                                                                                                                                  |                                                                                                                                                     
tAACTTGAACGCTTGTTCAATGAAAAACGCTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATg                                                                                                                                                     <  7:297991/148‑1 (MQ=255)
tAACTTGAACGCTTGTTCAATGAAAAACGCTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATg                                                                                                                                                     <  7:24491/148‑1 (MQ=255)
tAACTTGAACGCTTGTTCAATGAAAAACGCTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATg                                                                                                                                                     <  6:36808/148‑1 (MQ=255)
   cTTGAACGCTTGTTCAATGAAAAACGCTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGACTGCACCGAACCCTTCCCTTGCCAATGGCaa                                                                                                                                                 <  1:28219/149‑1 (MQ=255)
    ttGAACGCTTGTTCAATGAAAAACGCTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGc                                                                                                                                                   >  6:238334/1‑146 (MQ=255)
      gAACGCTGGTTCAATGAAAAACGCTGGTCCGCTTCCGGTGTTCAACAGGTGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCATTGGCAAg                                                                                                                                                <  2:95975/147‑1 (MQ=255)
                  aaTGAAAAACGCTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAAccc                                                                                                                                  >  5:8921/1‑149 (MQ=255)
                             cTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCt                                                                                                                      <  6:8921/149‑1 (MQ=255)
                              tGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCATTGGCAAgggaa                                                                                                                                            <  1:268268/127‑5 (MQ=255)
                              tGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCATTGGCAAgggaa                                                                                                                                            >  2:268268/1‑123 (MQ=255)
                                gTCCGCTGCCTGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACAGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTAc                                                                                                                    <  2:238515/148‑1 (MQ=255)
                                    gCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAAcccc                                                                                                                                 <  8:190078/132‑1 (MQ=255)
                                         cGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTgcggc                                                                                                            >  4:220447/1‑147 (MQ=255)
                                            tgtTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGtt                                                                                                       <  4:84748/149‑1 (MQ=255)
                                                 aaCAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATc                                                                                                  <  3:220447/149‑1 (MQ=255)
                                                    aggaggATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGgtgt                                                                                                                         >  4:48051/1‑123 (MQ=255)
                                                    aggaggATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGgtgt                                                                                                                         <  3:48051/123‑1 (MQ=255)
                                                    aggaggATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCt                                                                                                >  2:157059/1‑148 (MQ=255)
                                                     ggaggaTTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCa                                                                                              <  5:238334/149‑1 (MQ=255)
                                                                 tGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTg                                                                                   <  4:146918/148‑1 (MQ=255)
                                                                 tGCGGAACCCGGTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTg                                                                                   <  4:297057/148‑1 (MQ=255)
                                                                  gCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTgg                                                                                  <  1:157059/148‑1 (MQ=255)
                                                                     gAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGgt                                                                                                                           <  5:310027/104‑1 (MQ=255)
                                                                     gAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGgt                                                                                                                           >  6:310027/1‑104 (MQ=255)
                                                                        cccGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTgcg                                                                                                              >  4:116903/1‑114 (MQ=255)
                                                                        cccGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTgcg                                                                                                              >  6:192582/1‑114 (MQ=255)
                                                                         ccGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCAtc                                                                          >  1:171991/1‑149 (MQ=255)
                                                                         ccGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACACCGGGGGGTTcatgcg                                                                                                                  >  6:329485/1‑104 (MQ=255)
                                                                         ccGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACACCGGGGGGTTcatgcg                                                                                                                  <  5:329485/109‑6 (MQ=255)
                                                                           gTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCAtctc                                                                        <  2:171991/149‑1 (MQ=255)
                                                                                 tCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACTGCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAAcc                                                                   >  1:261723/1‑149 (MQ=255)
                                                                                   cGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTATACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCg                                                                 >  6:43615/1‑148 (MQ=255)
                                                                                            cGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTATACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAg                                                       <  5:43615/149‑1 (MQ=255)
                                                                                             gTCATTCGACGAGGTGCACGCCCTATGAGTGCCCCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCAGCAAGCCGCAgg                                                      >  2:234734/1‑149 (MQ=255)
                                                                                             gTCATTCGACGAGGTGCACGCCCAATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAgg                                                      >  6:66911/1‑149 (MQ=255)
                                                                                                    gacgaGGTGCACGCCCAATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTg                                                <  5:66911/148‑1 (MQ=255)
                                                                                                          gTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCt                                          >  8:321394/1‑148 (MQ=255)
                                                                                                          gTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCt                                          <  7:321394/148‑1 (MQ=255)
                                                                                                            gCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGt                                                                                                        <  5:312124/84‑1 (MQ=255)
                                                                                                            gCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGt                                                                                                        >  6:312124/1‑84 (MQ=255)
                                                                                                            gCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGt                                                     >  3:255510/1‑135 (MQ=255)
                                                                                                            gCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGt                                                     <  4:255510/135‑1 (MQ=255)
                                                                                                                  ccTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGACCCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAgcg                                   >  3:18157/1‑147 (MQ=255)
                                                                                                                  ccTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCACCCCGCAAGCCGCAGGTGACTGGTTGCTGGCAgcg                                   >  2:331428/1‑147 (MQ=255)
                                                                                                                   cTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTg                                         >  7:165670/1‑140 (MQ=255)
                                                                                                                   cTATGAGTGCACCGAACACTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTg                                         <  8:165670/140‑1 (MQ=255)
                                                                                                                             aCCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGggcgg                      >  3:317871/1‑149 (MQ=255)
                                                                                                                                      cttcccttGCCAATGGCAAGATCCGCATGAACCCCGCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTg                                            >  1:53005/1‑120 (MQ=255)
                                                                                                                                      cttcccttGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTg                                            <  2:53005/119‑1 (MQ=255)
                                                                                                                                       ttcccttGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCaa                                                              <  1:322583/100‑1 (MQ=255)
                                                                                                                                       ttcccttGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCaa                                                              >  2:322583/1‑100 (MQ=255)
                                                                                                                                       ttcccttGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGtta                                                   <  7:192873/111‑3 (MQ=255)
                                                                                                                                       ttcccttGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGtta                                                   >  8:192873/1‑109 (MQ=255)
                                                                                                                                       ggatcttGCCATTGGCAAGGTCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACg              <  6:335899/144‑1 (MQ=255)
                                                                                                                                       ggatcttGCCATTGGCAAGGTCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACg              >  5:335899/5‑148 (MQ=255)
                                                                                                                                          ccttGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACGGTc           <  1:234734/148‑1 (MQ=255)
                                                                                                                                           cttGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACGGTCgg         >  5:355725/1‑149 (MQ=255)
                                                                                                                                           cttGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACGGTCgg         >  7:59885/1‑149 (MQ=255)
                                                                                                                                              gCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACGGTCGGc        <  4:18157/147‑1 (MQ=255)
                                                                                                                                              gCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACGGTCGGCTg      >  5:30954/1‑149 (MQ=255)
                                                                                                                                                  aTGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACGGTCGGCTGGtac  <  6:355725/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TAACTTGAACGCTTGTTCAATGAAAAACGCTGGTCCGCTGCCGGTGTTCAACAGGAGGATTCCTTTGCTGAACCCGTTCATTCCGAACCTCACGTCATTCGACGAGGTGCACGCCCTATGAGTGCACCGAACCCTTCCCTTGCCAATGGCAAGATCCGCATGAACCCCCCGGTGTTCTACTTTGCGGCGAGTTTCATCCTCATCTTTGGCCTGGTAGTCATCTCCAACCCGCAAGCCGCAGGTGACTGGTTGCTGGCAGCGCAGAACTGGGCGGCCAACACGGTCGGCTGGTAC  >  NC_002947/282433‑282726

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