Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,462,029 T→C 11.2% G379G (GGT→GGC) ‡ 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,0290TC11.2% 251.0 / 18.1 89G379G (GGT→GGC) ‡paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (36/41);  new base C (5/5);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

TGGAACCGCGCCTGATGGGGCTTGGTCCGGTGCCTGCAACTCGCCGGGTGCTGGAACGCGCAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGAC  >  NC_000913/1461761‑1462301
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                
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                                                                                                                                                                                                                                                                              <  6:219706/149‑1 (MQ=255)
                                                                                                                             gggTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGtt                                                                                                                                                                                                                                                                              <  7:298626/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)
                                                                                                                                       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                                                                                                                                                                                                                                                     >  7:157365/1‑149 (MQ=255)
                                                                                                                                                      ccTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCgg                                                                                                                                                                                                                                                     >  6:287342/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                                                                                                                                                                                        >  7:169909/1‑149 (MQ=255)
                                                                                                                                                                                                                   gAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTAcc                                                                                                                                                                                        >  5:230599/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                                                                                                                                                                              <  5:81583/110‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aatgCGTAACGACCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGAt                                                                                                                                                                              >  6:81583/4‑113 (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)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                
TGGAACCGCGCCTGATGGGGCTTGGTCCGGTGCCTGCAACTCGCCGGGTGCTGGAACGCGCAGGGCTGAGTATTCACGATATGGACGTGATTGAACTGAACGAAGCGTTCGCGGCCCAGGCGTTGGGTGTACTACGCGAATTGGGGCTGCCTGATGATGCCCCACATGTTAACCCCAACGGAGGCGCTATCGCCTTAGGCCATCCGTTGGGAATGAGTGGTGCCCGCCTGGCACTGGCTGCCAGCCATGAGCTGCATCGGCGTAACGGTCGTTACGCATTGTGCACCATGTGCATCGGTGTCGGTCAGGGCATCGCCATGATTCTGGAGCGTGTTTGAGCATATCAACCTGCGAGTACCCTACAATGATAACCAATACAAAGCTTGACCCGATTGAAACCGCGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGAC  >  NC_000913/1461761‑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: