Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,462,028 G→A 12.2% G379D (GGT→GAT) ‡ paaJ → 3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,462,0280GA12.2% 236.4 / 19.2 88G379D (GGT→GAT) ‡paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (36/41);  new base A (5/6);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

CCGGGGTGGAACCGCGCCTGATGGGGCTTGGTCCGGTGCCTGCAACTCGCCGGGTGCTGGAACGCGCAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGAC  >  NC_000913/1461755‑1462301
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                 
caggggTGGAACCGCGCCTGACGGGGGTTGGACCGGTGCCTGCAGCTGGCCGGGTGCTGGAACGCGCAGGGCGGAGTATTCACGATATGGACGTGATTGAACTGAACGAGGCGTTGGCGGCCCAGGCGTTGGGTGTACTACGCGAATGGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACgg                                                                                                                                                                                                                                                                                   <  4:231733/272‑1 (MQ=255)
      tGGAACCGCGCCTGATGGGGCTTGGTCCGGTGCCTGCAACTCGCCGGGTGCTGGAACGCGCAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTAc                                                                                                                                                                                                                                                                            <  1:63950/275‑1 (MQ=255)
          acctggcCTCATGGGGCGTGGTCGGGTGCCTGCAACTCGCCGGGTGATGGAACGCGCAGGGGTGAGTATTCAAGATATGGACGTGATTGACGTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCAGCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCAt                                                                                                                                                                                                                                                                        <  4:372669/270‑1 (MQ=255)
              cgcCTGATGGGGCTTGGCCCGGTGCCTGCAACTCGCCGGGTGCTGGAACGCGCAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCCGCCATGAGCTGCATCGGCGTACCGGTCGTTACGCAtct                                                                                                                                                                                                                                                                      >  2:231737/1‑271 (MQ=255)
                   gATGGGGCTTGGTCCGGTGCCGGCAACTCGCCGGGGGCTGGAACGCGTAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGGTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCAcc                                                                                                                                                                                                                                                                <  2:228673/274‑1 (MQ=255)
                                                                cgcAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTTCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCCTCGGTGTCGGTCCGGGCATCGCCATGATTCCGGAGCgt                                                                                                                                                                                                                   >  2:276832/1‑274 (MQ=255)
                                                                  cAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTAc                                                                                                                                                                                                                                                                            >  4:354975/1‑215 (MQ=255)
                                                                  cAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCgtgt                                                                                                                                                                                                                 <  2:476976/274‑1 (MQ=255)
                                                                     gagtataagagagcGATATGGGAGTGATTGAGCTGAACGAAGCGTTCGCGGCCCAGGCGTGGGGGGTACTACGCGAATTGGGGGTGCCTGGTGATGCCCCACATGTTAGCCCCAACGGAGGCGCTATCGCCTTAGGCAATCCGTTGGGAATGGGTGGTGCCCGCCTGGCGCTGGCTGCCAGCCATGGGCTGCATCGGCGTGACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTg                                                                                                                                                                                                              <  4:121829/261‑1 (MQ=255)
                                                                      gCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCAt                                                                                                                                                                                                                                                                        <  2:354980/215‑1 (MQ=255)
                                                                               tCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATTTCaa                                                                                                                                                                                                    >  2:45733/1‑274 (MQ=255)
                                                                               tCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGACATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGGGCATAGGTGTAGGTCAGGGAATCGCCATGCTTCTGGAGCGTGTTTGAGCATATCAAc                                                                                                                                                                                                   >  2:347971/1‑275 (MQ=255)
                                                                                    attgggcagttaTTGAACTGAACGAAGCGTGCGCGGCCCGGGGGTTGGGTGGACTACGCGAATTGGGGCGGCCGGAGGATGCCCCACAGGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGc                                                                                                                                                                                               <  4:461562/263‑1 (MQ=255)
                                                                                       tGGACGTGTGTGGACTGGACTAAGCGTTCGCGCCCCAGGCGTGGGGTGTCCCACGCGAATTGGGGCTGCCTGAGTATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGGGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAg                                                                                                                                                                                            <  4:226037/274‑1 (MQ=255)
                                                                                         tgcGTGATGTGAAATGAACGAGGCGGTCGGGGCTCAGGGGTTGGGGGACCGACGGCGATGGGGGCTGGTGGATGATGGCCCAAATGTTAACCCCAACGGAGGCGCTATCGCCTAGGGCCATCCGTGGGGAATGGGTGGTGGCCGCCTGGGACGGGATGCCGGCCATGAGTTGGATCGGCGTAATGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGGCATGAGTCTGGAGCGTGTTTGAGGATATCAACCTGCGAGt                                                                                                                                                                                           <  4:491237/272‑1 (MQ=255)
                                                                                               aTTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGCGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTTNGTACCATGGGNATc                                                                                                                                                                                                                                                       >  4:528833/1‑207 (MQ=255)
                                                                                               aTTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATGGCCATGATTc                                                                                                                                                                                                                           >  4:49430/1‑235 (MQ=255)
                                                                                               aTTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGTTTCTGGAGCgtg                                                                                                                                                                                                                  >  2:121830/1‑244 (MQ=255)
                                                                                                    aactgaacGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGgtg                                                                                                                                                                                                                                                   <  2:528838/207‑1 (MQ=255)
                                                                                                    aactgaacGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGa                                                                                                                                                                                                                       <  2:49430/235‑1 (MQ=255)
                                                                                                         aacGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCCTTGTTCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGCTTCTTGAGTGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAAcc                                                                                                                                                                          >  1:64723/1‑275 (MQ=255)
                                                                                                                          ccAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATAACAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGAcc                                                                                                                                                          >  2:324819/1‑274 (MQ=255)
                                                                                                                                   gggTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGtt                                                                                                                                                                                                                                                                              <  7:298626/149‑1 (MQ=255)
                                                                                                                                   gggTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGtt                                                                                                                                                                                                                                                                              <  6:219706/149‑1 (MQ=255)
                                                                                                                                       gtACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGACCGTTACGc                                                                                                                                                                                                                                                                          <  7:322645/149‑1 (MQ=255)
                                                                                                                                           tacGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATtg                                                                                                                                                                                                                                                                      <  6:274938/149‑1 (MQ=255)
                                                                                                                                           tacGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGACCGTTACGccgat                                                                                                                                                                                                                                                                      <  5:229045/149‑5 (MQ=255)
                                                                                                                                             cgcgAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACgg                                                                                                                                                                                                                                                                                   <  4:505021/134‑1 (MQ=255)
                                                                                                                                             cgcgAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATtgtg                                                                                                                                                                                                                                                                    >  5:380203/1‑149 (MQ=255)
                                                                                                                                             cgcgAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCt                                                                                                                                                                                                                          >  4:75693/1‑191 (MQ=255)
                                                                                                                                               cgAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATc                                                                                                                                                                                                                                     >  4:407380/1‑178 (MQ=255)
                                                                                                                                                 aaTTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAg                                                                                                                                                                                                                      <  2:75694/191‑1 (MQ=255)
                                                                                                                                                  aTTGTGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGAATCGCCATGATTCTGGAGCGTGCTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCTCGTCTGTtgatga                                                                                                                                 >  2:357015/1‑275 (MQ=255)
                                                                                                                                                   ttGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCa                                                                                                                                                                                                                                 <  2:407385/178‑1 (MQ=255)
                                                                                                                                                       ggCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGc                                                                                                                                                                                                                                                          >  6:58176/1‑149 (MQ=255)
                                                                                                                                                          tGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATc                                                                                                                                                                                                                                                       <  6:467129/149‑1 (MQ=255)
                                                                                                                                                            ccTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCgg                                                                                                                                                                                                                                                     >  6:287342/1‑149 (MQ=255)
                                                                                                                                                            ccTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCgg                                                                                                                                                                                                                                                     >  7:157365/1‑149 (MQ=255)
                                                                                                                                                              tgatgaTGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGgtg                                                                                                                                                                                                                                                   >  7:406661/1‑149 (MQ=255)
                                                                                                                                                                          acaTGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCa                                                                                                                                                                                                                                       >  6:3738/1‑149 (MQ=255)
                                                                                                                                                                           actGTTACCCCCAACCGAGGGGCTATCGCCTTAGGCATTCCTTTGGGAATGAGTGGTGCCCGCCTGGCACGGGCTGCCAGCCATGCGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCg                                                                                                         <  4:324815/272‑1 (MQ=255)
                                                                                                                                                                                   ccccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTTAGGGCATCGCCATgg                                                                                                                                                                                                                              >  6:182824/1‑148 (MQ=255)
                                                                                                                                                                                    cccAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGAt                                                                                                                                                                                                                             >  6:74446/1‑149 (MQ=255)
                                                                                                                                                                                            aGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGc                                                                                                                                                                                                                     <  5:58176/149‑1 (MQ=255)
                                                                                                                                                                                                  tatctgctgaGGCCATGCGTTGGGAATGAGGGGTGACCGACTGGCACTGGCTGACAGCCATGAGCGGCAGCGGCGGAACGGTGGGTACGCATTGTGCACAAGGTGCATCGGGGTCGGTCAGGGCAGCGCCATGATGCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACAcgc                                                                                 <  1:154470/266‑1 (MQ=255)
                                                                                                                                                                                                  tATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGttt                                                                                                                                                                                                               <  5:3738/149‑1 (MQ=255)
                                                                                                                                                                                                    tCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCAtgt                                                                                                                                                                                                                                                            >  4:285038/1‑102 (MQ=255)
                                                                                                                                                                                                    tCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTGAGGGCATCGCCATGATTCTGGAGCGTGTTTGa                                                                                                                                                                                                             <  5:182824/149‑1 (MQ=255)
                                                                                                                                                                                                    tCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGa                                                                                                                                                                                                             <  6:203608/149‑1 (MQ=255)
                                                                                                                                                                                                       ccTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCa                                                                                                                                                                                                          >  5:442830/1‑149 (MQ=255)
                                                                                                                                                                                                       ccTTAGGCCATACGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCa                                                                                                                                                                                                          >  5:250345/1‑149 (MQ=255)
                                                                                                                                                                                                        cTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCAt                                                                                                                                                                                                                                                        <  2:285042/102‑1 (MQ=255)
                                                                                                                                                                                                           aGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATc                                                                                                                                                                                                      <  7:214499/149‑1 (MQ=255)
                                                                                                                                                                                                            ggCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCTTCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCa                                                                                                                                                                                                     >  7:85455/1‑149 (MQ=255)
                                                                                                                                                                                                              ccATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAAc                                                                                                                                                                                                   <  6:380203/149‑1 (MQ=255)
                                                                                                                                                                                                                aTCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCACATCAACCt                                                                                                                                                                                                 <  7:53231/149‑1 (MQ=255)
                                                                                                                                                                                                                 tCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTTGCTTCCAGCCATGAGCTGCACCGGCGTAACGGTCGTTACGCATAGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATTATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATCCGAAGCTTGACCCGATTGAAACCGCAGCTGTTGATTAGTTACAGCAGTTGCAACCCCAGCGTCGCAAGAGTAAGCCAAAACACCCGTATGgaactcgtcc                                                                  >  4:233418/1‑268 (MQ=255)
                                                                                                                                                                                                                  ccGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGc                                                                                                                                                                                               <  8:157365/149‑1 (MQ=255)
                                                                                                                                                                                                                       gggAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCg                                                                                                                                                                                              >  2:309838/1‑145 (MQ=255)
                                                                                                                                                                                                                         gAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAcc                                                                                                                                                                                        >  5:230599/1‑149 (MQ=255)
                                                                                                                                                                                                                         gAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAcc                                                                                                                                                                                        >  7:169909/1‑149 (MQ=255)
                                                                                                                                                                                                                           aTGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTa                                                                                                                                                                                          <  4:309834/145‑1 (MQ=255)
                                                                                                                                                                                                                               gtggtgCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCCTCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTAAAATGATAACCAATACAAAGCTTGACCCCATTGAAACCGCGTCCGTTGATGAGTTCAAGGCGTTGCAAACACAGCGTCCGGAATGGACGCCCAAACACGCGTATGAAAATGTGCCGATGTATAgggcca                                                    >  4:189164/1‑271 (MQ=255)
                                                                                                                                                                                                                                      ccGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAAc                                                                                                                                                                           <  6:250345/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ccTGGCACTGACTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATCTCAACCTGCGAGTACCCTACAATGATAACCaa                                                                                                                                                                        <  8:406661/149‑1 (MQ=255)
                                                                                                                                                                                                                                             gCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAGCCTACAATGATAACCAATAc                                                                                                                                                                     >  8:344704/1‑148 (MQ=255)
                                                                                                                                                                                                                                             gCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATAc                                                                                                                                                                     <  7:344704/148‑1 (MQ=255)
                                                                                                                                                                                                                                              cACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACgcagca                                      <  2:189166/274‑1 (MQ=255)
                                                                                                                                                                                                                                                     tGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGa                                                                                                                                                            <  5:74446/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     atccagcccTGTGATGGATCGGCGTAAAGAGCGTTACGGAATGTGGACCGTGGTAATCGGTGTCGGACAGGGCATCAGCATGATGCTGGAGCGTGTTTGAACATATCAACCTGCGAGTACCCTACAATGAGAACCAGTAAAAAGCTTGACCCAATTGAAACCGCGTCTGTTGATGAGTTACAGGCGGTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACa                              <  4:276828/273‑1 (MQ=255)
                                                                                                                                                                                                                                                        cagccaTGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGAccc                                                                                                                                                         >  5:393778/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          gccaTGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCGATACAAAGCTTGACCCGa                                                                                                                                                       <  5:361416/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCt                                                                                                                                                               >  2:227194/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                      aTCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATTATAACCAATACAAAGCTTGACACGATTGAAACCGCGt                                                                                                                                           >  6:33095/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      aTCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGa                                                                                                                                                       <  8:147931/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                      aTCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGa                                                                                                                                                       >  7:147931/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                       aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGAt                                                                                                                                                                              >  6:81583/4‑113 (MQ=255)
                                                                                                                                                                                                                                                                       aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGAt                                                                                                                                                                              <  5:81583/110‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTg                                                                                                                                                             >  5:24490/4‑130 (MQ=255)
                                                                                                                                                                                                                                                                       aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTg                                                                                                                                                             <  6:24490/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGt                                                                                                                                           <  8:458680/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGt                                                                                                                                           >  7:458680/4‑148 (MQ=255)
                                                                                                                                                                                                                                                                       aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTc                                                                                                                                          >  7:162742/4‑149 (MQ=255)
                                                                                                                                                                                                                                                                         tgCGTAACGACCGTTACGCAGTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTg                                                                                                                                        >  6:121053/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                         ggCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGAc                                                                                                                                                           <  4:227190/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ggCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTg                                                                                                                                        <  6:393778/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gCGTAACGGTGTTTCCGCGGTGTGCACCATGTGCATCGGTGTCGGTCAGGGGATCGTCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGACGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGa            <  1:135910/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gCGTAACGGTCGGTACGCATTGTGCACCAGGTGCATAGGGGTCGGTTAGGGAATCGCCATGATGCTGGAGCGTGTTTGAGCATAGCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGa            <  1:3543/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gCGTAACGGGCGTTACGCATTGGGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGa            <  2:297969/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                              aaCGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGGTTGAAACCGCGTCTGTtgat                                                                                                                                   >  5:245634/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                  gTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGTCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGAc  <  2:498573/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                 
CCGGGGTGGAACCGCGCCTGATGGGGCTTGGTCCGGTGCCTGCAACTCGCCGGGTGCTGGAACGCGCAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGAC  >  NC_000913/1461755‑1462301

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