Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,715,401 A→T 8.1% 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 8.1% 347.2 / 16.1 112D166E (GAT→GAA) paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (46/56);  new base T (2/7);  total (49/63)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 123 total aligned reads displayed.

GGATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGC  >  NC_002947/3715253‑3715548
                                                                                                                                                    |                                                                                                                                                   
ggATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTACGCTCGTGTGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                                    <  1:1978284/150‑1 (MQ=255)
 gATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaatc                                                                                                                                                   <  2:767323/150‑1 (MQ=255)
 gATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaatc                                                                                                                                                   <  1:398499/150‑1 (MQ=255)
    ttCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGAt                                                                                                                                                <  2:247224/150‑1 (MQ=255)
     tCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATg                                                                                                                                               <  1:805818/150‑1 (MQ=255)
      cGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGc                                                                                                                                              >  1:1095055/1‑150 (MQ=255)
        aCCGGCACGCTTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                            <  2:441156/150‑1 (MQ=255)
        aCCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                            <  2:105883/150‑1 (MQ=255)
        aCCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCg                                                                                                                                            <  1:2133052/150‑1 (MQ=255)
          cGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTa                                                                                                                                          >  1:1048742/1‑150 (MQ=255)
          cGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTa                                                                                                                                          >  2:332360/1‑150 (MQ=255)
           ggCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCggaa                                                                                                                                         <  2:2364291/150‑4 (MQ=255)
            gCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAgg                                                                                                                                        >  1:2074457/1‑150 (MQ=255)
             cACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCggaaac                                                                                                                                       <  1:1815045/150‑6 (MQ=255)
              aCGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAg                                                                                                                                         <  1:764665/147‑1 (MQ=255)
              aCGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAg                                                                                                                                         >  2:764665/1‑147 (MQ=255)
               cGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCggaaacgg                                                                                                                                     <  1:2305712/150‑8 (MQ=255)
               cGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAggcgg                                                                                                                                     <  1:1517923/150‑1 (MQ=255)
                  tttCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCt                                                                                                                                  <  1:209866/150‑1 (MQ=255)
                     cccGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTtca                                                                                                                               <  2:607506/150‑1 (MQ=255)
                         cgcCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAg                                                                                                                           <  2:2367733/150‑1 (MQ=255)
                          gcCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGc                                                                                                                          <  1:987996/150‑1 (MQ=255)
                           ccAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                         <  2:100954/150‑1 (MQ=255)
                           ccAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                         <  1:1879730/150‑1 (MQ=255)
                           ccAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                         <  1:2313883/150‑1 (MQ=255)
                             aGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCCATGCCCTACGCAGCCACCACCAGCgcg                                                                                                                       >  2:138409/1‑148 (MQ=255)
                             aGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTAGGCCACTTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCggg                                                                                                                       <  1:297349/150‑1 (MQ=255)
                               gcggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGtt                                                                                                                     <  1:616668/150‑1 (MQ=255)
                                  gcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGAc                                                                                                                  <  1:332360/150‑1 (MQ=255)
                                      cgGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGTTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGaaa                                                                                                              <  2:1272045/150‑1 (MQ=255)
                                           cataGCTGCTGGCGCTGTCTTATGGCAGCGCAGGGGAAAGGAATTCGGGTCGGTGCGCGAGATGCCGAAGTATTCTGCCACGTTTTCCGCCCTTTTTGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGTGGGTTTACGAAACGCCa                                                                                                         <  2:337590/147‑1 (MQ=255)
                                           cTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCa                                                                                                         <  2:2074457/150‑1 (MQ=255)
                                                ctgctgCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                         <  2:1577825/129‑1 (MQ=255)
                                                ctgctgCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCg                                                                                                                         >  1:1577825/1‑129 (MQ=255)
                                                   ctgcGGATTTGTGCTGGCTGCGCAGGGCAACAGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtgg                                                                                                 >  1:467361/1‑150 (MQ=255)
                                                    tgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                            >  2:109199/1‑106 (MQ=255)
                                                    tgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                            <  1:109199/106‑1 (MQ=255)
                                                     gcGGCTTTGTGCTGGCTGAGAAGGGTCAAGGCATCCTGTTCGGCGCGCGAGATGTCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                               <  2:1832630/150‑1 (MQ=255)
                                                         cTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGgtgt                                                                                           >  1:1683733/1‑150 (MQ=255)
                                                          tttGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATggtggt                                                                                             >  2:851252/1‑147 (MQ=255)
                                                          tttGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATggtggt                                                                                             <  1:851252/147‑1 (MQ=255)
                                                          tttGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTc                                                                                          <  2:1576967/150‑1 (MQ=255)
                                                             gtgCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGATGATGTcata                                                                                       <  2:601791/150‑4 (MQ=255)
                                                               gctggctgCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACa                                                                                     <  2:570636/150‑1 (MQ=255)
                                                                ctggctgCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                                    <  2:948044/86‑1 (MQ=255)
                                                                ctggctgCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaat                                                                                                                                                    >  1:948044/1‑86 (MQ=255)
                                                                 tggctgCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGc                                                                                   >  2:1695026/1‑150 (MQ=255)
                                                                  ggctgCGCAGGGCAACGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCt                                                                                  <  2:1683733/150‑1 (MQ=255)
                                                                          aGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCa                                                                          >  1:1048960/1‑150 (MQ=255)
                                                                               aaaGGCATCCTGGTCGGCGCGCGGGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGccc                                                                     >  2:761720/1‑150 (MQ=255)
                                                                               aaaGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGccc                                                                     >  2:2264414/1‑150 (MQ=255)
                                                                                 aGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGa                                                                   >  2:2127416/1‑150 (MQ=255)
                                                                                       ccTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGAACGCCCGAATGTCt                                                             <  2:2280024/150‑1 (MQ=255)
                                                                                           gTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTc                                                         >  2:1708757/1‑150 (MQ=255)
                                                                                            tCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTcc                                                        >  2:1795889/1‑150 (MQ=255)
                                                                                              ggCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGa                                                                                        <  2:1371707/116‑1 (MQ=255)
                                                                                                 gcgcgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACtt                                                   >  1:1047544/1‑150 (MQ=255)
                                                                                                   gcgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCgg                                                                                                                        <  1:793663/79‑1 (MQ=255)
                                                                                                       gaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGc                                                               <  1:913570/132‑1 (MQ=255)
                                                                                                       gaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGc                                                               >  2:913570/1‑132 (MQ=255)
                                                                                                       gaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTc                                                         <  2:2040500/138‑1 (MQ=255)
                                                                                                        aTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATc                                            <  1:140068/150‑1 (MQ=255)
                                                                                                         ggCCGAACAGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATca                                           >  1:696997/2‑150 (MQ=255)
                                                                                                          gCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTc                                                         >  1:1236710/1‑135 (MQ=255)
                                                                                                            cGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGaa                                                                  >  2:254215/1‑124 (MQ=255)
                                                                                                            cGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGaa                                                                  <  1:254215/124‑1 (MQ=255)
                                                                                                                 tGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGt                                   >  1:1184746/1‑150 (MQ=255)
                                                                                                                 tGTTCGGCCACGTTTTCAGCAGTTACCAAAATCGAATGTATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGt                                   <  2:633184/150‑1 (MQ=255)
                                                                                                                   ttCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                               <  2:1245032/88‑1 (MQ=255)
                                                                                                                   ttCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAtggtg                                                                                               >  1:1245032/1‑88 (MQ=255)
                                                                                                                   ttCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACgg                                    <  1:1793012/147‑1 (MQ=255)
                                                                                                                   ttCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACgg                                    >  2:1793012/1‑147 (MQ=255)
                                                                                                                     cGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAc                               <  1:2264414/150‑1 (MQ=255)
                                                                                                                       gCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgc                             >  2:153118/1‑150 (MQ=255)
                                                                                                                        ccACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGc                                                                                                       <  2:513843/75‑1 (MQ=255)
                                                                                                                        ccACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGc                                                                                                       >  1:513843/1‑75 (MQ=255)
                                                                                                                        ccACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgcg                            >  2:40516/1‑150 (MQ=255)
                                                                                                                        ccACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgcg                            >  2:1097010/1‑150 (MQ=255)
                                                                                                                         cACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAACCGCAAGCCGATGGGGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGa                           <  2:805246/150‑1 (MQ=255)
                                                                                                                         cACGTTGTCCGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGa                           <  1:1361134/150‑1 (MQ=255)
                                                                                                                          aCGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGAc                          >  2:334254/1‑150 (MQ=255)
                                                                                                                             ttttCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATg                       <  2:1184746/150‑1 (MQ=255)
                                                                                                                             ttttCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATg                       <  1:1708757/150‑1 (MQ=255)
                                                                                                                             ttttCAGCCGTTTCCGGCATCCAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAATCCAACGAAGAGGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATg                       <  1:2222535/150‑1 (MQ=255)
                                                                                                                                  aGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGa                                                                                                    <  2:2188663/68‑1 (MQ=255)
                                                                                                                                  aGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGa                                                                                                    >  1:2188663/1‑68 (MQ=255)
                                                                                                                                  aGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTc                  >  2:256944/1‑150 (MQ=255)
                                                                                                                                   gCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTcc                 >  2:343457/1‑150 (MQ=255)
                                                                                                                                   gCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTcc                 >  2:20841/1‑150 (MQ=255)
                                                                                                                                      gTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAAc                                      >  1:2120947/1‑126 (MQ=255)
                                                                                                                                       tttCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCacac             >  2:456960/1‑150 (MQ=255)
                                                                                                                                         tacGGCATCGATTCGATGTCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGa                                                      <  1:1932539/105‑1 (MQ=255)
                                                                                                                                           cGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAc                                                 >  1:1995107/1‑110 (MQ=255)
                                                                                                                                           cGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAc                                                 <  2:1995107/110‑1 (MQ=255)
                                                                                                                                           cGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACaccacc         >  2:326204/1‑150 (MQ=255)
                                                                                                                                               atcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGcc     <  2:1095055/150‑1 (MQ=255)
                                                                                                                                                tcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTa                                                  <  1:1687057/104‑1 (MQ=255)
                                                                                                                                                tcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCa    <  1:1795889/150‑1 (MQ=255)
                                                                                                                                                tcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCa    <  1:20841/150‑1 (MQ=255)
                                                                                                                                                  gaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGTACTCCACACCACCGGCCAGc  >  2:508411/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGC  >  NC_002947/3715253‑3715548

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