Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,715,401 A→T 6.8% D166E (GAT→GAA 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_0029473,715,4010AT 6.8% 182.4 / 4.3 59D166E (GAT→GAApaaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (22/33);  new base T (2/2);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GGGATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCA  >  NC_002947/3715252‑3715549
                                                                                                                                                     |                                                                                                                                                    
gggATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaa                                                                                                                                                      >  2:1258441/1‑150 (MQ=255)
gggATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaa                                                                                                                                                      <  1:270259/150‑1 (MQ=255)
 ggATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                                     <  1:886025/150‑1 (MQ=255)
 ggATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                                     <  1:1063845/150‑1 (MQ=255)
 ggATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                                     >  1:453391/1‑150 (MQ=255)
        gACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGcc                                                                                                                                              >  1:173569/1‑150 (MQ=255)
           cGGCACGATTTCCCGCGCCCGGCGGCGCCCGGCCCGAGCCGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTa                                                                                                                                           <  2:209932/150‑1 (MQ=255)
            ggCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAg                                                                                                                                          >  1:316626/1‑150 (MQ=255)
             gCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                             <  1:772418/146‑1 (MQ=255)
             gCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                             >  2:772418/1‑146 (MQ=255)
              cACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCCGGCATTTCCGGCATCGAATCTGGGCCGTTggc                                                                                                                                        >  2:916005/1‑150 (MQ=255)
              cACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAggc                                                                                                                                        >  1:676965/1‑150 (MQ=255)
                 gATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAggcggc                                                                                                                                     >  2:755738/1‑150 (MQ=255)
                  aTTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGcc                                                                                                                                    <  2:1222560/150‑1 (MQ=255)
                   tttCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCt                                                                                                                                   >  2:572233/1‑150 (MQ=255)
                    ttCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCtt                                                                                                                                  <  2:676965/150‑1 (MQ=255)
                       ccGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTtcat                                                                                                                               >  1:691066/1‑150 (MQ=255)
                        agcgcCATGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTtcatc                                                                                                                              >  2:1282314/2‑150 (MQ=255)
                            ccAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                          <  1:318951/150‑1 (MQ=255)
                            ccAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                          >  1:511695/1‑150 (MQ=255)
                                  ggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                          >  2:625424/1‑144 (MQ=255)
                                  ggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                          <  1:625424/144‑1 (MQ=255)
                                  ggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGa                                                                                                                    <  1:4735/150‑1 (MQ=255)
                                   gcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGAc                                                                                                                   <  1:551037/150‑1 (MQ=255)
                                       cgGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGaaa                                                                                                               <  1:572233/150‑1 (MQ=255)
                                        gggCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAAc                                                                                                              <  1:815719/150‑1 (MQ=255)
                                          gCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGc                                                                                                            <  2:173569/150‑1 (MQ=255)
                                            cTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGtt                                                                                                                      <  1:861735/138‑1 (MQ=255)
                                            cTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCAGCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGtt                                                                                                                      >  2:861735/1‑138 (MQ=255)
                                             tGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAg                                                                                                         <  1:1258441/150‑1 (MQ=255)
                                                 ctgctgCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGa                                                                                                     >  1:745328/1‑150 (MQ=255)
                                                  ggcttcGGTTTTGTGCTGGCTGCGCAGGGTAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAt                                                                                                    <  2:511695/149‑1 (MQ=255)
                                                      gcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                                <  1:322154/150‑1 (MQ=255)
                                                              gtgCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGaa                                                                                        <  1:1756/150‑1 (MQ=255)
                                                               tgCTGGCTGCGCAGGCCAAAGGCATCTTGGTGGGCGCGCGGGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCCAATCGATGCCGTAGGCGGCCTTCAGCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAAc                                                                                       <  2:386218/150‑1 (MQ=255)
                                                                     ctgCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCg                                                                                 <  1:261089/150‑1 (MQ=255)
                                                                        cgcAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGcc                                                                              <  1:488490/150‑1 (MQ=255)
                                                                              gCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGaaa                                                                                                               <  2:799078/111‑1 (MQ=255)
                                                                              gCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGaaa                                                                                                               >  1:799078/1‑111 (MQ=255)
                                                                                aaaGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGccc                                                                      >  2:465729/1‑150 (MQ=255)
                                                                                  aGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGa                                                                    >  2:1085916/1‑150 (MQ=255)
                                                                                   ggCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGaa                                                                   <  2:300342/150‑1 (MQ=255)
                                                                                       tCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGcc                                                               <  1:916005/150‑1 (MQ=255)
                                                                                       tCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGcc                                                               <  2:316626/150‑1 (MQ=255)
                                                                                                    gcgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATTCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAc                                                  <  2:1274080/150‑1 (MQ=255)
                                                                                                      gagaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAccc                                                <  2:548099/150‑1 (MQ=255)
                                                                                                         aTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATc                                             <  1:204107/150‑1 (MQ=255)
                                                                                                             cGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTa                                                   <  1:54930/140‑1 (MQ=255)
                                                                                                             cGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTa                                                   >  2:54930/1‑140 (MQ=255)
                                                                                                                 cTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAAAAGCTAGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACgg                                     >  1:461996/1‑150 (MQ=255)
                                                                                                                   gTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTg                                   >  1:463418/1‑150 (MQ=255)
                                                                                                                           aCGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGgtgt                                                                                            >  1:362987/1‑85 (MQ=255)
                                                                                                                           aCGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGgtgt                                                                                            <  2:362987/85‑1 (MQ=255)
                                                                                                                            cGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACa                          <  1:795320/150‑1 (MQ=255)
                                                                                                                                ttCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCCCGCGCCATGGa                      <  2:691066/150‑1 (MQ=255)
                                                                                                                                   aGCCGTTTCCGGCATCGAATCGATGCCGGATGCGGCCTTCATCAGCGGGGTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTc                   >  2:79037/1‑150 (MQ=255)
                                                                                                                                     ccGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAAc                                                                                                              >  2:447735/1‑57 (MQ=255)
                                                                                                                                     ccGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAAc                                                                                                              <  1:447735/57‑1 (MQ=255)
                                                                                                                                     ccGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCca                 <  1:465729/150‑1 (MQ=255)
                                                                                                                                        tttCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCacac              <  1:755738/150‑1 (MQ=255)
                                                                                                                                          tacGGCATCGATTCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtg                                                                                                   >  2:711765/3‑63 (MQ=255)
                                                                                                                                          tacGGCATCGATTCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtg                                                                                                   <  1:711765/61‑1 (MQ=255)
                                                                                                                                             ggCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCg         <  1:131199/150‑1 (MQ=255)
                                                                                                                                              gCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCgg        >  2:352567/1‑150 (MQ=255)
                                                                                                                                                atcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATca                                            <  2:578715/112‑1 (MQ=255)
                                                                                                                                                atcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATca                                            >  1:578715/1‑112 (MQ=255)
                                                                                                                                                atcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGACGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGcc      <  1:1139650/150‑1 (MQ=255)
                                                                                                                                                 tcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCa     <  1:1282314/150‑1 (MQ=255)
                                                                                                                                                   gaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGc   <  1:1085916/150‑1 (MQ=255)
                                                                                                                                                    aatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCa  >  2:1028731/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCA  >  NC_002947/3715252‑3715549

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